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Preface 


The  primary  characteristic  of  wireless  communications  is  variability.  Any  change  in  the  posi¬ 
tion  of  the  transmitter  or  receiver,  or  in  fact,  of  any  object  in  the  vicinity  of  the  wireless  system 
can  cause  the  quality  of  the  transmission  to  change  This  effect  is  referred  to  as  fading.  Fad 
ing  is  caused  by  the  transmitted  signal  making  its  way  to  the  receiver  via  multiple  paths,  and 
these  multipath  components  can  then  destructively  interfere  to  significantly  reduce  the  received 
power.  Another  important  characteristic  of  wireless  communications  is  that  users  must  share  a 
common  resource,  bandwidth. 

Consideration  of  these  problems  and  limitations  led  to  the  development  of  multiple-input 
multiple-output  (M1MO)  wireless  communication  systems  in  the  mid  1990s  and  this  technique 
has  emerged  as  the  most  promising  technique  for  improving  the  performance  of  wireless  dig¬ 
ital  transmission  systems  [1,2].  This  improvement  is  two-fold.  First,  the  limited  resources  of 
a  wireless  communication  system,  such  as  spectrum  and  power,  can  be  efficiently  used  with 
multiple  antennas  to  provide  good  quality  and  large  capacity  to  a  wide  range  of  applications 
requiring  high  data  rates.  In  short,  M1MO  techniques  allow  higher  data  rates  to  be  obtained 
for  a  given  bandwidth.  The  second  advantage  obtained  through  the  use  of  Ml  MO  systems  is 
robustness  to  fading  caused  by  multipath  propagation  of  the  radio  transmissions  As  a  conse¬ 
quence  of  these  benefits,  MIMO  has  gone  through  the  adoption  curve  for  commercial  wireless- 
systems  to  the  point  that  today  (2009),  all  high  throughput  commercial  standards  (i.e.  WiMax, 
Wi-Fi  etc.)  have  adopted  MIMO  as  part  of  the  optional,  if  not  mandatory,  portions  of  their 
standards  [3] 

A  consequence  of  the  MIMO  wireless  communications  revolution  is  that  multipath  scatter¬ 
ing  of  the  transmitted  signal  has  moved  from  being  a  problem  to  being  a  valuable  resource. 
This  has  happened  to  such  an  extent  that  MIMO  systems  have  come  to  rely  on  this  resource. 
However,  away  from  dense  urban  and  indoor  environments,  for  example,  in  Australian  rural 
and  remote  areas  with  wide  open  spaces  and  flat  or  smooth  undulating  terrain  such  rich  scat¬ 
tering  can  be  hard  to  find. 

The  research  reported  in  this  work  seeks  to  develop  MIMO  schemes  which  will  give  robust 
and  reliable  performance  in  environments  that  can  change  rapidly  from  a  rich  scattering  en¬ 
vironment  to  clear  line  of  sight  between  transmitter  and  receiver  Wireless  communications 
systems  which  arc  robust  to  changing  environmental  conditions  arc  a  particularly  important 
factor  in  military  communication  systems  as  well  as  for  civilian  emergency  services. 

Central  to  this  research  program  is  the  introduction  of  a  new  degree  of  freedom  at  the  phys¬ 
ical  layer  of  the  wireless  network,  called  polarization.  The  polarization  of  the  electromagnetic 
wave  carrying  the  information  from  the  transmitter  to  the  receiver  is  defined  by  the  orientation 
of  the  electric  field  vector  in  the  plane  perpendicular  to  the  direction  of  propagation.  Tradition¬ 
ally,  wireless  communication  systems  have  used  a  single  transmit  antenna  and  a  single  receive 
antenna  with  linear  vertical  polarization  Linear  horizontal  and  circular  polarization  are  also 
possible  and  have  found  application  in  wireless  point-to-point  (e.g.  IEEE  802.1 1  WiFi)  and 
satellite  communication  systems  respectively.  By  introducing  either  dual-polarized  antennas 
and/or  multiple  antennas  with  different  polarization,  the  dimension  of  polarization  is  added  to 
those  of  time,  frequency  and  space  for  wireless  transmission  at  the  physical  layer  [4], 

Some  existing  wireless  systems  (e.g.  many  WiFi  access  points)  exploit  polarization  in  an 
implicit  way.  These  implicit  techniques  are  analogous  to  traditional  selection  diversity  [4]. 
The  signal  is  launched  with  a  single  polarization,  the  propagation  medium  then  couples  some 
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energy  into  the  cross-polarization  plane,  and  the  signal  is  received  with  a  cross-polarized  an¬ 
tenna  [5]  Exploitation  of  implicit  polarization  diversity  is  suited  to  rich  scattering  (e.g.  dense 
urhan  and  indoor)  env  ironments  [6].  In  this  work,  polarization  is  proposed  as  a  source  of  ex¬ 
plicit  diversity  in  wireless  communications.  Explicit  polarization  diversity  is  particularly  suited 
to  the  Australian  rural  and  remote  environment  where  wide  open  spaces  and  flat  or  smooth  un¬ 
dulating  terrain  give  rise  to  line  of  sight  (LOS)  conditions  between  the  transmitter  and  the 
receiver  The  polarization  diversities  of  transmitted  signals  are  mostly  preserved  by  the  LOS 
environment,  presenting  an  opportunity  for  transmitter  diversity  as  well  as  receiver  diversity 
techniques. 

Another  aspect  of  the  research  reported  here  is  to  develop  fast,  fixed  complexity,  decod¬ 
ing  algorithms  for  space-time  codes  that  are  robust  to  the  changing  environment.  The  major 
drawback  of  spacc-timc  codes  is  the  high  complexity  of  their  decoding  algorithms.  This  is 
particularly  true  for  a  numher  of  codes  which  are  best  suited  to  the  exploitation  of  polarization 
diversity.  The  existing  “fast"  decoding  algorithm  for  this  problem,  the  so-called  sphere  de¬ 
coder ,  has  performance  w  hich  depends  crucially  on  the  channel  conditions.  When  the  channel 
is  close  to  singular,  that  is,  when  the  channel  between  the  hase  station  and  terminal  is  close  to 
pure  line  of  sight,  the  sphere  decoder  defaults  to  an  exhaustive  search.  If  such  conditions  per¬ 
sist  the  communication  system  could  be  in  outage  purely  due  to  computational  overload.  The 
fast  decoding  algorithms  developed  here  have  fixed  complexity  across  all  channel  conditions. 

The  report  is  essentially  divided  into  three  parts.  The  first  part,  Chapters  1  and  2,  give  the 
necessary  background  on  MIMO  communications  and  space-time  block  codes  and  goes  on  to 
develop  some  new  techniques  which  we  later  use  to  analyse  polarimetric  MIMO  schemes. 

The  second  part  consists  of  Chapters  3  and  4  It  reports  on  our  development  of  fast  de¬ 
coding  algorithms  for  space-time  codes.  We  develop  fast  decoding  algorithms  for  a  number 
of  important  space-time  codes,  that  have  fixed  complexity  across  channel  conditions.  Our 
fast  decoding  algorithms  are  based  on  a  technique  called  conditional  optimization ,  which  is 
widely  used  in  statistical  estimation  and  signal  processing  to  reduce  the  search  space  of  an 
optimization  prohlem  by  taking  advantage  of  the  possibility  of  analytically  optimizing  over 
some  subset  of  the  parameters,  conditioned  on  the  remaining  parameters.  This  technique  has 
hardly  hecn  exploited  in  discrete  optimization  problems  associated  with  decoding  in  wireless 
communications.  We  show  how  to  apply  this  method  to  a  large  class  of  full  rate,  full-diversity 
space  time  block  codes,  which  includes  multiplexed  orthogonal  designs,  to  give  fast  maximum 
likelihood  (ML)  detection  with  low  computational  complexity  The  technique  is  also  applied  to 
the  Golden  Code  and  perfect  space-time  hlock  codes,  in  an  approximate  form,  to  obtain  essen¬ 
tially  maximum  likelihood  performance  with  greatly  reduced  complexity.  Chapter  4  analyses, 
in  a  general  way,  the  structure  of  space-time  codes  that  allow  fast  maximum  likelihood  decod¬ 
ing.  We  show  that  the  perfect  space-time  codes  of  Oggier  et  al.  [7]  are  in  fact,  multiplexed 
quasi-orthogonal  designs  and  so  have  fast,  fixed  complexity,  exact  ML  decoding  algorithms. 

The  third  part  of  this  report  explores  the  benefit  of  polarization  diversity  for  MIMO  sys¬ 
tems.  In  Chapter  2  we  analyse  interference  cancellation  for  multiple  Alamouti  schemes  using 
a  Bayesian  approach  which  provides  a  unified  framework  with  which  to  understand  the  rela¬ 
tionship  between  various  signal  detection  techniques.  The  approach  leads  to  a  new  parameter 
which  can  be  used  to  predict  the  performance  of  systems  hased  on  Alamouti  coding,  with¬ 
out  needing  to  resort  to  simulations.  This  parameter  is  used  to  analyse  various  transmission 
schemes  involving  polarization  diversity  in  Chapter  5.  We  show  that  a  certain  scheme  which 
multiplexes  space  time  codes  across  polarization  gives  significant  henefits  in  stability  of  per- 
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formance  for  the  Ml  MO  system  in  the  presence  of  even  a  small  LOS  component.  In  particular, 
this  scheme  is  shown  to  give  inherent  stability  with  respect  to  the  relative  orientation  of  a  pair 
of  spatially  separated  dual  polarized  transmit  antenna.  This  provides  an  important  indieation 
of  how  polarization  can  be  used  to  develop  more  robust  and  reliable  wireless  communication 
systems.  We  go  on  to  investigate  systems  using  dual  polarized  transmit  and  dual  polarized 
receive  antenna  in  various  conditions  from  pure  LOS,  Kieean  and  pure  scattering.  It  is  shown 
that  the  performance  of  the  Golden  Code  can  be  made  stable  across  propagation  conditions  by 
the  use  of  dual  polarized  transmit  and  dual  polarized  receive  antennas.  However,  this  is  true 
only  if  the  transmitter  and/or  receiver  arc  fixed  in  position  or  are  allowed  to  rotate  in  a  plane  of 
alignment. 

In  mobile  communications,  the  transmitter  and/or  receiver  rotate  with  respect  to  each  other, 
throwing  the  dual  polarimctne  antennas  out  of  alignment  and  degrading  the  performance  sig¬ 
nificantly.  To  improve  the  stability  in  transmission  and/or  reception  under  the  rotation  of  the 
transmitter  and/or  receiver,  wc  investigate  the  use  of  a  triad  antenna  at  the  transmitter  and  the 
receiver  in  Chapter  6.  Wc  analytically  show  that  the  capacity  of  a  system  using  triad  antennas 
is  preserved  under  the  relative  rotation  of  the  transmitter  and  receiver.  This  is  not  the  case  for 
the  dual  polarized  system.  We  introduce  a  3  x  3  space -time  code  suitable  for  the  triad  system 
which  gives  full  rate,  through  which  the  gains  achievable  with  triad  antennas  are  demonstrated 
This  code  also  has  the  advantage  of  being  fast  to  decode  due  to  our  fast  decoding  algorithm 
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Notation 


E  Expectation 
IFl  /i  x  7i  Identity  matrix 
Q(  )  Quantiser  (Slicer) 
det  Determinant  of  matrix 
Z  Integer 

Z[t]  Gaussian  integer 
C  Complex  number 
R  Real  number 
Q  Rational  number 

C  A  xA/  N  x  M  matrices  of  complex  number 
!Rn  77.  dimensional  vectors  of  real  numbers 
IV  Trace  of  matrix 
diag  Diagonal  of  a  matrix 
*  complex  conjugate 
t  Hermitian  conjugate 
bR(  )  real  part 


Abbreviations 


AWGN 

Additive  White  Gaussian  Noise 

BLAST 

Bell  Lab  Layered  Space-Time  Architecture 

CDMA 

Code  Division  Multiple  Access 

CSI 

Channel  State  Information 

FDGE 

Enhanced  Data  for  Global  Evolution 

GSM 

Global  System  for  Mobile 

1C 

Interference  Cancellation 

LDC 

Linear  Dispersion  Code 

LOS 

Linc-of-Sight 

MAP 

Maximum  a  Posteriori 

MIMO 

Multiple-Input  Multiple-Output 

ML 

Maximum  Likelihood 

OD 

Orthogonal  Design 

OSTBC 

Orthogonal  Space-Time  Block  Code 

PCS 

Personal  Communication  System 

QOD 

Quasi-Orthogonal  Design 

QOSTBC 

Quasi-Orthogonal  Space-Time  Block  Codes 

SINR 

Signal  to  Interference  and  Noise  Ratio 

SM 

Spatial  Multiplexing 

SNR 

Signal  to  Noise  Ratio 

STBC 

Space-Time  Block  Code 

UMTS 

Universal  Mobile  Telecommunications  System 

WiFi 

Wireless  Fidelity 

WiMAX 

Worldwide  Interoperability  for  Microwave  Access 

XPD 

Cross  Polarization  Discrimination 

ZF 

Zero  Forcing 

Chapter  1 :  MIMO  Communications 


1.1  Introduction 

Recent  advances  in  wireless  communication  systems  have  provided  high  speed,  high  through¬ 
put  over  wireless  channels,  providing  high  quality  information  exchange  between  portable  de 
vices  located  anywhere  in  the  world  Wired  communication  has  more  stability,  better  perfor¬ 
mance  and  higher  reliability,  but  it  greatly  constrains  the  environments  or  locations  in  which 
the  system  can  be  used.  Wireless  communication  provides  the  benefits  of  portability,  mobility 
and  accessibility.  However,  this  represents  both  freedom  to  the  end  user  and  also  a  number  of 
challenges  for  the  system  designer  [8]. 

Wireless  communication  is  the  transfer  of  information  over  a  distance  without  the  use  of 
electrical  conductors  or  wires  Instead,  the  information  rides  on  electromagnetic  waves,  with 
the  consequence  that  the  information  undergoes  attenuation  (fading)  due  to  the  interaction  of 
the  electromagnetic  waves  with  the  physical  environment.  The  signal  attenuation  by  wireless 
propagation  may  be  due  to  the  distance  between  communicating  nodes,  referred  to  as  path  loss , 
or  due  to  shadowing  from  obstacles  such  as  buildings,  or  due  to  constructive  and  destructive 
interference  of  multiple  reflections  of  the  electromagnetic  wave,  referred  to  as  multipath  prop¬ 
agation.  Figure  1.1  illustrates  a  typical  outdoor  wireless  propagation  environment  where  the 
mobile  node  is  communicating  with  a  wireless  access  point  (base  station).  The  signal  trans¬ 
mitted  from  the  mobile  may  reach  the  base  station  directly  (line  of  sight)  or  through  multiple 
reflections  on  local  scatterers  such  as  buildings,  trees,  hills  etc.  The  received  signal  is  a  combi¬ 
nation  of  multiple  delayed  copies  of  the  transmitted  signal.  This  multipath  propagation  makes 
wireless  a  challenging  communication  environment. 


Figure  LI:  Radio  ptvpagalioti  environment 


Fortunately,  this  challenge  can  be  turned  into  an  opportunity  through  the  use  of  multiple 
antennas.  The  multiple  antennas  may  be  in  spatially  separate  locations  and/or  have  different 
polarizations  ( polarization  diversity ).  Due  to  different  local  scattering  environments,  suffi¬ 
ciently  separated  antenna  elements  provide  almost  independent  fading  channels.  How  much 
separation  is  appropriate  depends  on  the  environment.  In  rural  macro  cells,  many  wavelengths 
of  separation  may  be  required  to  de-correlate  antennas,  whereas  in  indoor  environments  a  half 
wavelength  separation  may  be  sufficient  [9].  With  sufficient  spacing,  multiple  antennas  at 
both  ends  of  a  wireless  link  can  be  used  to  exploit  the  statistical  independence  of  the  multi¬ 
ple  channels  connecting  transmitter  to  receiver  to  increase  spectral  efficiency  thought  the  uses 
of  spatial  multiplexing.  The  existence  of  multiple  independent  channels  can  also  be  used  to 
improve  link  reliability  through  the  use  of  transmit  diversity  techniques  A  transmit  diversity 
technique  which  spreads  the  transmitted  information  across  the  multiple  independent  paths, 
and  which  is  particularly  attractive  because  it  does  not  require  channel  state  information  (CSI) 
at  the  transmitter,  is  space- time  coding  [10]. 


1.2  Spatial  Multiplexing 

Spatial  multiplexing  (SM)  is  a  technique  which  uses  multiple  antennas  at  both  the  transmitter 
and  receiver  to  increase  the  transmission  rate  by  exploiting  multipath.  That  is,  it  uses  the 
scattering  characteristics  of  the  propagation  environment  to  enhance  the  transmission  rate  by 
treating  the  multiplicity  of  scattering  paths  as  separate  parallel  sub-channels.  Bell  Labs  was 
the  first  to  demonstrate  a  laboratory  prototype  in  1998,  known  as  BLAST  (Bell  Labs  Layered 
Space-Time  Architecture,  sec  Figure  1 .2),  where  SM  is  the  principal  technique  used  to  improve 
the  performance  of  a  wireless  communication  system.  Bell  Labs  accomplished  this  by  splitting 
a  single  user’s  data  stream  into  multiple  sub-streams  and  using  an  array  of  transmitter  antennas 
to  simultaneously  launch  the  parallel  sub-streams.  All  the  sub-streams  are  transmitted  in  the 
same  frequency  band,  so  spectrum  is  used  very  efficiently.  Since  the  user’s  data  is  being  sent  in 
parallel  over  multiple  antennas,  the  effective  transmission  rate  increases  roughly  linearly  with 
the  number  of  antennas. 


Transmit  Antennas  Receive  Antennas 


Figure  1.2:  MI  MO  system 


2 


At  the  receiver,  if  the  multipath  scattering  is  sufficient,  i.e.  the  channel  matrix  has  full  rank, 
then  the  received  signal  at  eaeh  receive  antenna  is  a  non-zero  linear  combination  of  indepen 
dent  transmitted  sub-streams.  Using  sophisticated  signal  processing,  these  slight  differences  in 
scattering  allow  the  sub-streams  to  be  separated  and  recovered  In  its  most  basic  form  known 
as  vertical  (V)-BLAST,  the  detection  process  uses  a  combination  of  nulling  and  symbol  can¬ 
cellation  to  successively  compute  the  decision  statistics  which  are  then  used  to  form  estimates 
of  the  underlying  data  symbols. 

Spatial  multiplexing  relies  on  the  scattering  richness  in  the  propagation  environment.  The 
lack  of  such  conditions  reduces  the  achievable  rates  of  spatial  multiplexing  techniques  from 
their  theoretical  projections  under  ideal  assumptions  [11,  12]. 


1.3  Space-Time  Block  Codes 


Spaee-time  block  codes  (STBC)  introduced  by  Tarokh  et  al.  [10],  improve  the  reliability  of 
communication  over  fading  channels  by  correlating  signals  across  different  transmit  antennas. 
The  transmitted  data  stream  is  encoded  in  bloeks  which  are  distributed  across  space  and  time. 

Consider  a  communication  system  with  nt  transmit  and  nT  receive  antennas.  Assume  that 
the  transmitted  symbols  from  the  nt  transmit  antennas  at  time  slot  /  is 

x,  =  (x},x? . *;")  i  <  t  <  i  (i.i) 


where  xj  is  the  transmitted  symbol  from  antenna  i  at  time  slot  and  all  transmitted  symbols 
have  the  same  duration.  A  nt  x  r  space-time  codeword  matrix  can  be  defined  as 


X  =  (xi ,  x2, - xT) 


xj 

T2 

Xi 


x\ 

xl 


(12) 


The  ith  row  x1  =  (xi,x2, . . .  ,  x^)  is  the  data  sequence  transmitted  from  the  iih  transmit 

antenna  and  the  jth  column  x3  =  (xj,x^, - a:"')  is  the  space-time  symbol  transmitted  at 

time  j,  1  <  j  <  7  .  In  a  general  form,  an  STBC  can  be  seen  as  a  mapping  of  T  complex 
sy  mbols  (x\ ,  x2,  ,  xt)  onto  a  matrix  X  of  dimension  nt  x  7 


(xi ,  X2 . xt)  — »  X 


(1.3) 


1.3.1  Diversity 

Error  probability  on  an  additive  white  Gaussian  noise  (AWGN)  channel  decays  exponentially 
with  SNR,  and  the  challenge  of  communication  over  Rayleigh  fading  channels  is  that  error 
probability  deeays  only  inversely  with  SNR.  The  received  symbols  y  are  given  by 

y  —  hX  -f  ?i  (1.4) 

where  X  is  the  transmitted  codeword,  h  is  the  vector  of  fading  (or  channel)  coefficients,  and  w 
is  Gaussian  noise.  We  assume  a  quasi-statie  model  where  the  fading  coefficients  are  constant 
over  a  frame  and  ehange  independently  from  one  frame  to  the  next  If  we  normalize  X  so  that 
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the  average  energy  per  complex  symbol  is  1  (SNR  =  l/Aro),  then  the  probability  of  error  given 
perfect  channel  state  information  is 


Pr  (Xi 


V  I 


and  averaging  over  the  thannel  statistics  gives 


Pr(X, 


X;) 


y/SNR/i(X,  X,)(X, 


(1.5) 


(1.6) 


where 

g<1)-  is  F”*  n,) 

The  matrix  (X,  X;)(Xi  —  XJ)t  is  Hermitian,  so  there  is  a  unitary  matrix  U  for  which 


l/(X,  -  X,)(X,  -  Xj)'U'  =diag(A?,-.- ,A*,)  (1.8) 

where  the  entries  A,  are  the  singular  values  of  X,  X;.  The  change  of  basis  h  —>  Uh 
preserves  the  Rayleigh  distribution  of  h  and  if  Ai ,  •  •  •  ,  An,  are  nonzero,  then 


Pr(X,  -*Xj)  = 


(SNR/2)-"' 
A?  •A*, 


(1.9) 


and  the  matrix  Xi  X,  is  nonsingular  and  maximal  diversity  is  achieved.  We  can  associate 
the  following  quantities  with  a  spaec-timc  code: 


•  The  diversity  gain  is  the  minimum  of  the  rank  of  X!  X;  over  all  codewords  in  the 
code. 


•  The  coding  gain  is  the  minimum  value  of  the  product  of  the  non-zero  singular  values  of 
X,  -  X;  over  all  codewords. 


There  is  a  fundamental  trade-off  between  using  the  multipath  to  improve  throughput  by  mul¬ 
tiplexing  and  improving  reliability  through  using  a  space-time  code.  This  can  be  quantified 
using  the  follow  ing  definitions  of  diversity  order  and  multiplexing  rate. 

A  scheme  which  has  an  average  error  probability  Pe(SNR)  as  a  function  of  SNR  that 
satisfies 


lun 


SNR—  oc 


log 


7Jr(SNR) 
log  SNR 


-d 


(1.10) 


is  said  to  have  a  diversity  order  of  d.  In  other  words,  a  scheme  with  diversity  order  d  has  an 
error  probability  at  high  SNR  that  behaves  as  Pe(SNR)  ~  SNR  "d. 

A  scheme  which  has  a  transmission  rate  of  R( SNR)  as  a  function  of  SNR  is  said  to  have  a 
multiplexing  rate  r  if 


Urn 

SNR— oo  log  SNR 


(1.11) 


In  other  words,  the  system  has  a  rate  of  r  log  SNR  at  high  SNR.  There  is  a  trade-off  between 
achievable  diversity  and  multiplexing  rate  which  is  expressed  in  terms  of  the  supremum  d°p,(r) 
of  the  diversity  gain  achievable  by  any  scheme  with  multiplexing  rate  r.  The  following  theorem 
is  due  to  Zheng  and  Tsc  [13], 
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Theorem  1  Let  K  —  nr)  where  lit  and  nr  are  the  number  of  transmit  and  rec  eive 

antennas  and  suppose  that  the  code  block  length  N  >  nt  +  nr  —  1  Then  the  optimal  trade - 
off  curve  d>p,(r)  is  given  by  the  piecewise  linear  function  connecting  points  in  (fc,  cTp,(k)), 
k  =  0,  ■  •  •  ,  K  where 

d>pt(k)  =  (lit  -  k)(nr  -  k).  (1  12) 

It  implies  that  high  multiplexing  rate  comes  at  the  price  of  decreased  diversity  gain  and  high 
diversity  gain  comes  at  the  price  of  low  multiplexing  rate. 

We  now  give  a  description  of  some  of  the  most  important  space*time  codes. 

1.3.2  AlamoutiCode 

The  most  famous  space-time  block  code  was  discovered  by  Alamouti  [14].  and  is  described 
by  a  2  x  2  matrix  where  the  columns  represent  different  time  slots,  the  rows  represent  dif¬ 
ferent  antennas,  and  the  entries  are  the  symbols  to  be  transmitted.  The  Alamouti  codeword  is 
represented  as 

*-(-3  3) 

The  two  rows  (or  columns)  of  X  are  orthogonal  to  each  other, 

XX*  =  (|x,|2  +  |ar2|2)I2. 

Thus,  given  two  codewords  X  and  Y  we  see  that 

det(X  Y)(X-  Y)»  =  (|x,  -  y,|*  +  |*a  -  y2\2)2  (1.15) 

This  means  that  a  subcode  of  the  Alamouti  code  corresponding  to  symbol  pairs  x3  =  (:r,  x ,  xJ2)  6 
C2,  j  —  1 ,  •  •  •  ,  N,  has 

Diversity  gain-  2 

Coding  gain:  min^gf,  ..  ,N)  ||*j  -  •r».-||2 

At  the  receiver,  the  signals  n,  7*2  received  over  two  consecutive  time  slots  are  given  by 

(-3)  -  {-%  !:■;)  (8)  ♦  (8) 

where  /r  1-/12  are  the  path  gains  from  the  two  transmit  antennas  to  the  receiver,  and  the  noise 
samples  m,  112  are  independent  samples  of  a  /ero-mcan  complex  Gaussian  random  variable 
with  variance  2rr2  per  complex  dimension.  Equation  (1.16)  can  be  represented  as 

r  —  7ix  +  n  (1.17) 

where  the  equivalent  matrix  channel  (induced  channel)  Ti  is  orthogonal  i  c  TC^TC  =  I2. 
Throughout  this  work  we  consider  the  case  of  coherent  detection  i.e.  the  receiver  is  assumed 
to  have  perfect  knowledge  of  the  Channel  State  Information  (CS1).  In  practice  CS1  can  be 
obtained  by  introducing  some  pilot  transmissions  that  enable  accurate  channel  estimation.  With 
Ti  known,  the  following  processing  can  be  done  at  the  receiver 

H'r  =  (\hl\2  +  \h2\2)x  +  nl  (1.18) 


(1  13) 

(U4) 
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The  new  noise  term  n  is  still  white,  so  the  ML,  detection  of  x\  X2  is  decoupled. 

The  Alamouti  code  provides  full  diversity  at  full  transmission  rate  for  any  real  or  complex 
signal  constellation,  it  does  not  require  CS1  at  the  transmitter,  and  the  ML  decoding  requires 
only  linear  complexity  processing.  For  these  reasons,  the  Alamouti  code  has  been  adopted 
in  several  wireless  standards  such  as  WCDMA  [15]  and  CDMA2000  [16].  The  Alamouti 
code  also  facilitates  higher  data  rates  through  multiplexing  of  parallel  data  streams  and  the 
addition  of  a  second  antenna  at  the  receiver  that  performs  interference  cancellation.  Data  rates 
of  4  bits/s/Hz  have  been  demonstrated  for  several  wireless  channels  including  UTMS,  GSM, 
FDGE,  IEEE  802.1  In  and  IEEE  802.16  [17].  This  scheme  is  investigated  in  the  next  chapter 
with  a  view  to  extending  its  use  to  a  polarization  diverse  MI  MO  systems. 


1.3.3  Orthogonal  Space-Time  Block  Codes 

The  Alamouti  code  is  an  example  of  an  Orthogonal  Space-Time  Block  Code  (OSTBC)  for 
two  transmit  antennas  and  its  success  provided  an  impetus  for  investigating  the  existence  of 
OSTBCs  for  more  than  two  transmit  antennas.  OSTBCs  are  also  referred  to  as  orthogonal 
designs  (OD)  and  have  the  property  that  ML  decoding  is  linear.  At  the  same  time  these  codes 
achieve  full  diversity.  In  order  to  describe  OSTBCs,  first  consider  an  STBC,  written  in  the 
form  of  a  linear  dispersion  code  fLDC)[l8]  (note  that  LDC  are  discussed  in  some  detail  in 
Section  4.1) 

2T 

X  =  J2x*A'  (1.19) 

t  1 

where  At  is  a  nt  x  T  complex  matrix  and  is  a  set  of  real  scalar  symbols. 

Definition  1  A  STBC  (1.19)  is  an  OSTBC  if  it  satisfies  the  following 

2  r 

XX1  =  (1.20) 

i-i 


where  each  P,  is  a  diagonal  matrix.  It  follows  that 

AiA\  —  T>i,i  =  1 , 27’ 

A,  A]  +  AjA\  =  26ijVi ,  1  <  i  <  j  <  27’. 


(1.21) 


Tarokh  ct  al.  [19]  showed  that,  in  general,  full-rate  or  orthogonal  designs  exist  for  all  real 
constellations  but  only  for  two,  four,  or  eight  transmit  antennas,  while  they  exist  for  all  complex 
constellations  only  for  two  transmit  antennas  (the  Alamouti  scheme).  However,  for  particular 
constellations  it  might  be  possible  to  construct  orthogonal  designs  for  other  cases.  Moreover, 
if  a  rate  loss  is  acceptable,  orthogonal  designs  exist  for  an  arbitrary  number  of  transmit  an¬ 
tennas  [19].  It  is  also  shown  by  Wang  and  Xia  [20]  that  the  rate  of  a  generalized  complex 
orthogonal  design  cannot  exceed  3/4  for  more  than  two  antennas.  An  example  of  the  R  =  3/4 
rate  code  for  four  transmit  antennas  is  the  following: 


Xi 

^2 

X3 

0\ 

-X2 

*1 

0 

X3 

-*3 

0 

x\ 

-X2 

0 

X*2 

XlJ 

(1.22) 
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1.3.4  Quasi-Orthogonal  Space-Time  Block  Codes 

As  mentioned  above,  Tarokh  ct  al.  [19]  showed  that  full-rate  orthogonal  designs  with  complex 
elements  in  its  transmission  matrix  are  impossible  for  more  than  two  transmit  antennas.  This 
has  motivated  Jafarkhani  [21]  to  develop  a  full-rate  complex  design  for  four  transmit  antennas 
for  which  pairs  of  symbols  could  be  decoded  independently.  This  is  only  possible  if  one  accepts 
a  loss  of  diversity  compared  to  a  true  orthogonal  design  Jafarkhani  called  this  class  of  eodes 
quasi-orthogonal  space  time  block  codes  (QOSTBC)  and  introduced  the  4  x  4  eode 


(1.23) 


where  the  structure  of  the  2x‘2  bloeks  mirrors  that  of  the  Alamouti  code.  The  columns 
Xx, ?'  =  1  2,3,4  of  the  matrix  X  divide  into  two  groups  {Xi,X.i}  and  {X2.X3}  with 
the  columns  from  different  groups  being  orthogonal.  Jafarkhani  showed  that  the  maximum 
likelihood  decision  metric  is  a  sum  f(x  1,2*4)  +  g(x 2,2:3)  where  /  is  independent  of  X2,%3 
and  g  is  independent  of  xu  ^4.  The  decoder  hnds  the  pair  (x\,Xa)  that  minimizes  /(j'i,X4) 
and  (independently)  the  pair  (2*2, X3)  that  minimizes  g(x 2,2:3).  Thus  decoding  complexity  is 
quadratic  in  the  size  of  the  signal  constellation. 

The  main  issue  with  QOSTBCs  is  that  they  do  not  have  full  diversity.  However,  as  we 
will  discuss  in  Chapter  4  a  number  of  full-rate,  full-diversity  eodes  can  be  constructed  by 
multiplexing  quasi-orthogonal  STBCs.  We  show  that  codes  with  this  type  of  structure  admit 
fast  ML  decoding  algorithms  based  on  conditional  optimization.  Next  we  will  consider  other 
methods  for  design  of  STBCs. 

1.4  Design  Criteria  for  Space-Time  Block  Codes 

Under  the  assumption  that  perfect  channel  state  information  is  known  at  the  receiver,  Tarokh 
et  al  [10]  developed  the  following  two  design  criteria  for  the  high  SNR  regime. 

•  Rank  Criterion:  Maximize  the  diversity  gain,  i.e  the  minimum  rank  r  of  the  codewords 
difference  Xt  —  X3  overall  distinct  pairs  of  space  time  codewords  Xt,  X3. 

•  Determinant  Criterion:  For  a  given  diversity  r,  maximize  the  coding  gain,  i.c.  the 
minimum  product  of  the  nonzero  singular  values  of  the  difference  Xx  X;  over  all 
distinct  pairs  of  space-time  codewords  X*  X,. 

If  r  =  nt  for  all  pairs  (Xz,  Xj),  we  say  that  the  code  is  full  rank.  If  X  is  full  rank,  we  have 


clet(A,j)  ^  0  for  all  Ai3 


(1.24) 


where  A,;  =  (Xt  X;)(X7  X; )  * ,  and  we  say  that  the  code  has  full  diversity.  This  means 

that  we  can  exploit  all  the  nrnt  independent  channels  available  in  the  MIMO  system 

1.5  Algebraic  Space-Time  Block  Codes 

In  order  to  increase  reliability,  STB  eodes  are  designed  to  have  full  diversity  with  a  large 
minimum  determinant.  Such  eodes  can  be  constructed  by  looking  at  matrix  representations  of 
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algebraic  objects  such  as  Galois  helds  and  cyclic  division  algebras.  In  this  section,  we  give 
descriptions  of  this  class  of  codes  which  have  become  known  as  perfect  space-time  codes. 

An  interesting  result  of  our  research  is  that  these  perfect  codes  can  be  seen  as  multiplexed 
quasi-orthogonal  designs,  (see  Chapter  3  and  Chapter  4),  which  admit  fast  ML  decoding  algo¬ 
rithms.  The  most  important  of  these  codes  is  the  2  x  2  MIMO  system  known  as  the  Golden 
Code. 


1.5.1  The  Golden  Code 


The  Golden  Code  [22]  is  a  2  x  2  block  space-time  eode  that  encodes  four  complex  symbols 
over  two  time  slots  and  achieves  full  diversity.  It  is  known  to  be  the  best  2x2  STBC  in 
having  the  largest  coding  gain  It  achieves  a  trade  off  between  rate  and  reliability  for  both  one 
and  two  receive  antennas  that  is  the  best  possible  in  terms  of  the  diversity-multiplexing  bound 
derived  by  Zheng  and  Tse  [13]  In  fact,  the  Golden  Code  is  incorporated  into  the  IEEE  802.16c 
(WiMAX)  standard  [23]. 

Codewords  in  the  Golden  Code  take  the  form 


X  —  ( a  (  Xl  +  TX2  X3  +  TX4\ 
A  “  rfs  yo  aj  y(x3  4  pxi)  Xi  4-  fix 2 J 


(1.25) 


where  {a:}?  t  €  C  C  Z [i]  are  the  transmitted  symbols,  and  C  is  a  signal  constellation  taken  to 
be  2m -QAM.  The  parameters  r  =  (1  +  and  p  —  (l-\/5)/ 2  are  the  Golden  Ratio.  The 
diagonal  matrix  diagja.  a],  where  a  =  (1  +  ip),  a  =  (1  4-  ir ),  serves  to  equalize  transmitted 
signal  power  across  the  two  transmit  antennas.  The  entries  of  Golden  Code  codewords  are 
drawn  from  Z[i][\/5]  C  Q(t\  v  5).  Note  that  Z[f][\/5]  is  the  ring  of  elements  of  the  form 
(n i  4-  m2)  4  («3  4  i7u)\/5>ni  €  Z  and  Q(i,  \/5)  is  the  field  of  elements  of  the  form 
(aj  4  ia2)  4  (<23  4  i(i4 ) ,  at  €  Q.  Following  [24]  we  rewrite  (1.25)  as 


X 


l  fa  OX  T/ xi  x3\  fT  0\  f  x2  xA 

\A  v°  °y  iv13  Xl  J  +  v°  \lX4  X2J 


The  set  of  integer  matrices  of  the  form 


(,y  x)  <  *i»€Z[i]. 


(1.26) 


(1.27) 


forms  a  matrix  representation  or  the  cyclotomic  ring  Z[C8].  where  Cs  is  a  primitive  8th  root  of 
unity.  This  can  be  seen  explicitly  by  making  the  identification 


0.28) 

and  for  simplicity,  we  will  drop  the  subscript  and  write  £  = 

The  cyclotomic  ring  Z[C]  is  the  nng  of  algehraic  integers  in  the  cyclotomic  field  Q(0  and 
any  element  of  this  field  can  be  written  in  the  form 

a  =  Go  4  01 C  4  O2C2  T  a3^3,  (1.29) 


with  ao,  01 , 02,  03  G  Q.  The  identification  (1 .28)  extends  to  the  matrix  representation 


a  —  Aa 


(  a0  4  ia2  a  1  4  m3 
yz(ai  4  ia3)  Oo  4  ia2 


(1.30) 
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Q(/)  Q(n/2)  Q(?V2) 

Q 

Figure  1.3:  The  sub  fie  Id  structure  of  Q(C) 


The  field  Q(C)  is  a  Galois  extension  of  degree  4  over  Q  with  Gal(Q(0/Q)  —  the 
dihedral  group  of  order  4  There  are  four  field  automorphisms,  ,  j=  1 , 3,  5,  7,  given  by 

<MC)  =  CJ.  (1-31) 

In  fact,  Gal(Q(C)/Q)  =  {03,0$).  The  subfield  structure  of  Q(C)  is  shown  in  Hgure  1.3, 
where  we  note  that  Q2  =  i,  C  -  =  x/2  and  Q  +  C3  =  i\/ 2.  Given  any  element  a  6  Q(C), 

the  field  norm  is  defined  by 

Wo(C>/Q(a)  =  II  ^(a) 

tr6<»al(Q(0/Q)  (1.32) 

—  (flo  4-  «2 )2  +  (af  -ffl])2  4-  4aifl3(ao  —  ^2)  4"  4oo«2(«3  —  fli)* 

In  terms  of  the  matrix  representation  ( 1.30)  we  see  that 


A^Q(C)/q(°)  -  |<let(/la)|2. 

(1.33) 

In  fact,  in  terms  of  the  relative  norm  Aq(o/q(,)*  w  hich  is  defined  by. 

XlO/QoX)  =  o-l(o)<TS(o), 

(1.34) 

we  have 

^Q(C)/Q (i)(a)  ~  det(i4ft). 

(1.35) 

Another  norm  on  Q(C)  we  will  find  useful  is  defined  by 

ll«l|2  =  \  Oi(a)<77(tv)  +  (T3(rt)<75(a)) 

(1.36) 

=  a0  +  n2  +  fl!  +  (,3 

(1.37) 

=  {WaWl, 

(1.38) 

where  ||*||F  denotes  the  matrix  Frobenius  norm. 
The  Golden  Code  ean  be  written  in  the  form 


where 


B>  =  (o  S)  and  B2  =  -i(o  °)  (1.40) 

The  malriees  B\  and  B2  satisfy 

'IV  (SiBj)  =  5^,  (1.41) 

for  i,j  =  1,2.  A  consequence  of  the  structure  of  the  Golden  Code  is  that  for  any  two  code 
words  Xl  and  X;  their  difference  Xt  X_,  is  rank  2  and 

det(X,  -  Xj)(X,  -  X,)f  >  16/25.  (1.42) 

An  observation  that  may  potentially  lead  to  fast  decoding  algorithms  for  the  Golden  Code 
is  that  the  cycloiomic  ring  Z[C]  C  Q(C)  a  Euclidean  domain  with  respect  to  the  field  norm 
Nq(C)/Q*  We  re^er  l^e  rea(^cr  10  l^e  survey  of  cyelolomie  Euclidean  number  fields  by  Akhtar 
[25]  for  more  details.  The  observation  that  Q(C)  is  Euclidean  goes  back  to  Eisenslein  in  1850. 
Thus,  for  any  two  elements  of  a ,  b  6  Z[C]  with  b  ^  0,  there  exist  q ,  r  6  Z[£]  such  that 

a  =  qb  -f  r,  (1 .43) 

with  7  —0  or  Nq(Q/q(r)  <  (b).  In  fact,  we  have 

q  =  [ a/b ]  and  r  —  a  -  [a/b]  6,  (1.44) 

where  for  any  a  6  Q(C),  we  define  [a]  6  Z[C]  by 

l«l  =  L°ol  +  Lcti  1C  +  l«2lC2  +  LnslC3,  (145) 

where  for  x  6  R,  [x]  denotes  the  nearest  integer  to  x  The  Euclidean  algonthm  in  Z[C]  enables 
the  efficient  computation  of  a  greatest  common  divisor  for  any  pair  of  cyelolomie  integers  a, 
b  e  Z[£].  This  provides  the  basis  for  possible  fast  decoding  schemes  based  on  the  extension 
of  an  inhomogeneous  Diophantine  approximation  schemes  such  as  CasscTs  algorithm  [26]. 
Research  in  this  area  is  ongoing. 

Finally  it  is  interesting  to  note  that  the  IEEE  802. 16e  specifications  include  three  M1MO 
profiles  having  2  transmit  antennas,  from  which  two  are  defined  as  M1MO  schemes  on  both  the 
down-link  and  up-link  transmission  of  mobile  WiMAX  systems.  The  first  one  is  the  Alamouti 
code,  which  is  described  in  Section  1  3.2,  namely.  Matrix  A,  which  achieves  full  transmit 
diversity  at  the  expense  of  rate  loss.  The  second  one  is  spatial  multiplexing,  namely  Matrix  B, 
which  achieves  full  rate  at  the  expense  of  diversity  loss.  On  the  other  hand,  the  third  MIMO 
profile  with  two  transmit  antennas  in  IEEE  802. 16e  specifications,  known  as  Matrix  C,  is  both 
a  full-rale  and  full  diversity  code.  More  specifically.  Matrix  C  is  a  variant  of  the  Golden  Code. 

1.5.2  The  Pefect  Space-Time  Block  Codes 

One  of  our  aims  is  the  development  of  polarizalion-limc  codes  for  fully  polarimclric  MIMO 
systems.  For  this  application  we  need  codes  which  have  similar  properties  to  the  Golden  Code 
but  which  are  appropriate  for  more  than  2  transmit  antennas.  The  Perfect  space  lime  block 
codes  will  turn  out  to  be  very  good  candidates. 

Perfect  STBCs  were  first  introduced  by  Oggiere  cl  al  [7]  to  have  full-rale,  full-diversity, 
non-vanishing  determinant  and  to  be  information  lossless.  These  codes  are  constructed  for 
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2  x  2,3  x  3,4  x  4  and  6x6  M1MO  systems.  An  example  of  a  2  x  2  perfect  STBC  is  the 
Golden  Code.  The  following  are  descriptions  of  3  x  3  and  1  x  4  perfect  STBCs.  These  are 
written  in  a  novel  form  that  we  have  introduced  in  order  to  assist  in  the  development  of  our  fast 
decoding  algorithms  which  are  discussed  in  Chapter  3. 

3x3  Perfect  STBC:  The  perfect  3x3  STBC  transmits  nine  complex  (usually  from  an 
N-HEX  constellation)  information  symbols  {xt}?  i  over  three  time  slots  from  three  transmit 
antennas.  The  transmit  codewords  of  the  3  x  3  perfect  STBC  can  be  expressed  as 

2  /#3i  +  l  X3i+2  J"3tf3\ 

X  =  V>+1  j*3<+3  X3i+l  *3»-t2  I  (146) 

i^0  \jx 3i  +  2  jx 3t4-3  X3i+1/ 

where  the  diagonal  matrices  /?,  are 

Bi  =  (1  +  j)I3  +  0 

B2  =  (-1  -2j)I3  +  J02  (1.47) 

B3  =  (-1  -  2j)I3  +  (1  +i)0  +  (1  +  j)€>2 

with  0  -  <liag((?i,  62,03), Oi  —  2ros(2'7r/7).  j  =  e2l"/2.  The  B,  satisfy 

Vr(BjBm)  =  78 jm  (1.48) 

4  x  4  Perfect  STBC:  The  4  x  4  perfect  STBC  transmits  16  complex  (AT— QAM  constellation) 
information  symbols  {xi  }t16  L  over  four  time  slots  from  four  transmit  antennas.  The  codewords 
can  be  expressed  as 


x  =  I>+‘ 

i  0 


(X4tfl  ^4t^2 
iX 4i  f  4  X4t  4- 1 
lX4«+3  2X4 i 4.4 
iX.i,f2  ix 4i  f 3 


J4i  +  3  X4i  +  4\ 
X4i  f2  X4i  +  3  \ 
X41  fl  X41  +  2  I 

/X4t  +  4  X4,fl/ 


where 


—  (1  —  3i)I4  +  i&2 
Bi  =  ( 1  -  3i)0  +  i03 
Z?3  =  -il4  +  (-3  +  4t)©+  (1  -i)03 
/i4  =  (-1  +  i)I,  -  30  +  02  +  03 

with  0  =  <liag;(0i ,  62,  0$,  04 ),  0t  -  2cos(2*7r/15).  The  Bt  satisfy 

’l\(BjBl)  =  15  6jm 


(1.49) 


(1.50) 


(1.51) 


1.5.3  Multiplexed  Alamouti  Blocks 

One  approach  to  high  data  rates  is  to  multiplex  Alamouti  blocks  where  a  high  rate  space-time 
code  might  be  constructed  as  a  "linear"  combination  of  Alamouti  blocks.  These  codes  are  2x2 
STBC  with  full-rate,  full-diversity  and  have  the  structure  that  allows  fast  maximum  likelihood 
decoding  which  will  be  discussed  in  detail  in  Chapter  3. 
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Hottinen-Tirkkonen-Kashacv  [27,  28]:  These  authors  considered  the  2  x  2  full-rate,  full- 
diversity  STBC  which  has  the  form 


where  {xt}?  j  €  Z [*]. 


and 


f/ 


h  - 

*)+r 

r  ■ 

\X2 

x') 

V22 

C*  ■? 

)•  (a: 

)  =  t/( 

1  / 

'1  +  2  - 

1  +  2tn 

=  V7  ( 

1  +  2 1 

!  -i 

(1.52) 


(1.53) 


(1.54) 


which  is  a  unitary  matrix  This  code  is  also  know  n  as  the  Silver  Code  [29]. 

Sezginer-Sari  [30]:  Se/giner  and  Sari  proposed  a  full-rate,  full-diversity  2  x  2  STBC  which 
can  be  expressed  as 

^  ax i  T  bx 3  —  CX2  —  dx 4 
t  ax 2  +  bx 4  exf  +  dxj  . 


=  (n 

Vc 


(1.55) 


where  a,  6,  c  and  d  are  complex-valued  design  parameters  that  are  chosen  to  optimize  coding 
gain.  In  terms  of  the  transmitted  power,  Sezginer  and  Sari  have  expressed  the  desired  eondi 
tions  as 


M2  +  |6|2=  l  =  |r|2  +  |rf|2 

M2  +  |c|2  =  1  =  \b\2  +  \d\2 

The  first  condition  ensures  an  equal  transmit  power  at  each  symbol  time,  while  the  second 
condition  ensures  that  equal  total  power  is  transmitted  for  each  symbol.  These  conditions 
imply  that  all  the  design  parameters  should  have  the  same  magnitude,  i.e. 


|o|  =  Ml  =  kl  =  MI  =  i/v/2 


(1.57) 


Rabiei-AI-I)hahir  [31]:  These  authors  have  considered  a  full-rate,  full-diversity  STBC  for  2 
transmit  antennas  which  takes  the  form 

x  =  Allan  -  0iX2  013-3  +  oia-4\  (1.58) 

\a-2*3  -  023\i  02X,  +  023-2 J 

where  a\,a2<0i  and  02  arc  real  design  parameters  chosen  to  guarantee  that  the  code  is  an 
information  lossless  STBC  and  to  maximize  its  coding  gain  For  this  purpose,  the  following 
conditions  must  be  satisfied 


r  a2, + 0i  =1 

3(«i,/3i, 02.02)  6  R  such  that  <  o2  +  0?  =1 

I  Q1«2  -  0102  =  0 


(1.59) 
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1.6  Multiple-Antenna  Channel  Capacity 

One  important  aspect  of  the  use  of  STBCs  is  their  effect  on  the  capacity  of  the  communication 
channel  [32,  19,  33,  34]  That  is,  if  we  insist  on  using  a  particular  STBC  do  we  degrade 
the  channel’s  capacity  to  carry  information.  In  this  section  we  briefly  discussed  the  mutual 
information  and  capacity  of  a  M1MO  system  which  will  be  useful  for  later  discussion  on  system 
performance,  particularly  for  polarimetrie  channels. 

Consider  a  Ml  MO  Gaussian  channel  characterized  by  a  fixed  nr  x  nt  complex  matrix  H 
[kij  ],  where  htJ  is  the  eomplex  number  representing  path  gain  between  transmit  antenna  i  and 
receive  antenna  j.  At  eaeh  symbol  time,  the  nr -dimensional  receive  signal  veetor  is  represented 
by 

y  =  Hi  4  n  (1 .60) 

where  x  is  the  nr  dimensional  transmitted  veetor  and  the  noise  n  is  a  Gaussian  veetor  with 
nr  i.i.d.  components.  Assuming  that  the  vector  x  has  circularly  complex  Gaussian  distributed 
components  and  H  is  deterministic,  Telctar  [2]  showed  that  the  mutual  information  ean  be 
expressed  as 

/(x.y |H)  =  log2  (let.  (l,„  +  ^HQH')  (1.61) 

where  In,  is  the  nr  x  nr  identity  matrix,  a2  is  the  variance  per  real  dimension  of  the  noise  n 
and  Q  is  the  eovarianee  matrix  of  x. 


Q  =  E[  xx*]. 


(1.62) 


When  H  is  known  perfectly  at  the  transmitter,  then  the  mutual  information  can  be  optimized 
with  water-filling  [2,  35]  to  give  the  capacity 


C  = 


max 

Q.Tr{Q)<r* 


I(x.y) 


(1.63) 


where  Px  is  the  maximum  power  available  at  the  transmitter. 

Consider  the  case  where  the  channel  state  information,  H  is  perfectly  known  only  at  the 
receiver.  Teletar  [2]  conjectures  that  when  the  channel  matrix  is  random,  non-ergodie,  then  in 
the  high  SNR  region,  the  optimal  covariance  matrix  for  the  source  is  given  by 


Qupi  —  (Px/fit )I«i  • 

The  corresponding  mutual  information  is 

I{x  y)  =  log,  det  ( I„,  +  —  HHf 
V  »t 

where  p  =  Px/2<j 2  The  channel  capacity  is  then  given  by 


C  =  Eh  1og2 det  (l,lf  +  ^  HH'j 


(1.64) 


(1.65) 


(1.66) 


where  the  expectation  E  is  taken  over  the  possible  channel  matrices  H 
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1.6.1  Space-Time  Coding  and  Channel  Capacity 

Consider  a  MIMO  system  with  tit  transmit  and  nr  receive  antennas.  A  spaee-time  code  X 
represents  a  STBC  codeword  to  be  transmitted  from  nt  transmit  antennas,  extending  over  1 
time  slots.  The  ease  T  =  1  is  just  spatial  multiplexing.  The  n7  x  T  received  signal  matrix  is 
given  by 

Y  =  ,/^-HX  +  W  (1.67) 

V  Tit 

where  W  €  Cnr  is  the  additive  white  Gaussian  noise  matrix  with  entries  distributed  as 
CAT{ 0,  2a2). 

In  order  to  understand  the  impact  of  the  STBC  on  channel  capacity,  a  new  capacity  C  can 
be  derived  with  the  induced  (equivalent)  channel  matrix  To  this  end  we  first  rewrite  (1 .67)  as 


V  = 


,  -Fix  +  n 

V  nt 


(1.68) 


where  y  is  the  Tnr  x  1  received  signal  vector  and  x  is  a  vector  of  information  symbols  and  Ft 
is  an  equivalent  channel  matrix  Then  the  capacity  of  the  new  equivalent  channel  Ft  is  given 
by 


E 


n 


—  log2  clef  (  Im 


—mi' 

TU 


(1.69) 


If  C  =  C  in  equation  ( 1 .66)  then  the  STBC  X  is  described  as  information  lossless.  An  example 
of  an  information  lossless  code  is  the  Alamouti  code  for  a  2  x  1  MIMO  system. 
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Chapter  2:  Bayesian  Approach  to  Interference  Cancellation 


The  famous  space-time  block  code  discovered  by  Alamouti  [14]  is  a  remarkable  signaling 
seheme  as  it  allows  simple  maximum  likelihood  (ML)  decoding.  The  Alamouti  code  also 
facilitates  higher  data  rates  through  multiplexing  of  parallel  data  streams  and  the  addition  of 
a  second  antenna  at  the  receiver  that  performs  interference  cancellation  In  this  chapter,  we 
consider  interference  cancellation  for  multiple  Alamouti  sehemes  using  a  Bayesian  approach. 

Interference  cancellation  for  two  Alamouti  users  has  been  studied  and  analysed  previously 
in  [36,  32,  37].  Here  we  want  to  be  able  to  extend  this  two  Alamouti  seheme  to  more  complex 
and  varied  situations  and  in  particular  to  use  it  with  polarization  diverse  transmit  and  receive 
antennas.  Our  Bayesian  approach  provides  a  unified  framework  with  which  to  understand  the 
relationship  between  various  signal  detection  techniques  and  leads  to  a  new  parameter  which 
can  be  used  to  prediet  the  performance  of  the  Alamouti  seheme  system  without  needing  to 
resort  to  simulations. 

2.1  System  Model 

The  encoding  rule  for  the  Alamouti  STBC  is  described  by  a  2  x  2  matrix 


(2.1) 


where  the  rows  represent  different  antennas,  the  columns  represent  different  time  slots,  and  the 
entries  are  the  symbols  to  be  transmitted. 

Consider  two  eo-ehannel  users,  eaeh  using  the  Alamouti  code  Let  c  =  (jq,  X2)  and  s  = 
(x3,xa)  be  the  codewords  transmitted  by  the  first  and  second  users  respectively;  rtJ  is  the 
signal  received  at  the  antenna  i  at  time  slot  j ,  and  htJ  denotes  the  channel  gain  from  transmit 
antenna  i  to  the  receive  antenna  j.  The  received  signals  at  the  two  receive  antennas  can  be 
represented  as. 


(2.2) 


where  the  noise  samples  ntJ  are  independent  samples  of  a  complex  Gaussian  random  variable 
with  zero  mean  and  variance  2a2 .  hquation  (2.2)  can  be  rewritten  as 


r i  =  c.H i  +  sQx  4-  n, 
r2  =  C.y.2  +  S&2  +  712 


(2.3) 


where  rq  =  (rn,  -7q2),  7'2  =  (7*21 1  -r'22)  and 


(2.4) 
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Rewrite  (2.3)  as 

r  =  cH  +  sQ  +  n,  (2.5) 

where 

r  =  (n  r2),H  =  (H\  U2)  and  Q  =  (Qi  g2).  (2.6) 

The  problem  here  is  to  dctcet  c  and  s ,  The  challenge  is  that  the  solution  involves  a  great  deal 
more  computation  than  the  simple  Alamouti  scheme,  particularly  for  large  symbol  constella¬ 
tions,  due  to  cross  correlation  terms  in  the  likelihood  function. 

The  Alamouti  code  is  in  fact  a  matrix  representation  of  the  Quaternions.  Thus  the  analysis 
of  multiple  Alamouti  codes  is  facilitated  by  an  understanding  of  the  properties  of  quaternions 
and  quaternion  linear  algebra.  We  turn  to  this  now 


2.2  Properties  of  the  Quaternions 


The  quaternions  M  are  a  non-commutative  extension  of  the  complex  numbers.  A  general 
quaternion  can  be  written  as 

q  =  qo  +  q\  i  +  qij  +  q^k  (2.7) 

where  qo,  qi,  q2,qi  €  1R  and  i ,  j  and  k  satisfy  the  defining  relations 

t2=  j2  =  fc2  =  - 1,  (2.8) 

ijk=- 1  (2.9) 

The  conjugate  of  a  quaternion  q  is  defined  as 

q*  =  q0  -  q\i  -  q2j  -  q^k.  (2.10) 


The  set  H  is  closed  under  addition,  multiplication,  inversion  and  conjugation.  One  can  compute 
with  the  quaternions  using  these  relations  just  as  one  does  with  complex  numbers  using  i2  — 
1  The  quaternions  have  a  norm  defined  by 

IMI  =  v'o  =  Y  <7o  +  9?  +  <72  +  <73-  (2-1 1) 

Wc  can  also  define  the  real  part  of  a  quaternion  to  be  $t(q)  =  qo  and  note  that  for  any  quaternion 


q’  +  q  =  2'J  )(q). 


(2.12) 


and 


q*  -  q  q  =  R{q)-  (2.13) 

There  is  a  faithful  representation  of  H  in  terms  of  2  x  2  complex  matrices  defined  by 


qo  +  q\i  +  q2j  +  q$k  — >  ^ 

The  2x2  matrix  here  has  the  Alamouti  space-time  block  code  structure 


qo  +  iq i  <72  +  % 
-q2  t  w?3  qo  —  W\ 


y 


(2.14) 


f  h2\ 

2  h\) 


(2.15) 


with  h\  =  qo  +z‘<7i ,  h2  =  q2  +  iqi  €  C.  Note  that  in  (2.14)  we  identified  the  real  quaternions  q0 
with  real  multiples  of  the  identity  <7ol2  burthermorc,  it  can  be  seen  the  quaternion  conjugation 
corresponds  to  Hermitian  conjugation  for  the  matrix  representation.  If  H  is  a  matrix  of  the 
form  (2. 15).  the  quaternion  properties  (2  1 1 )— (2. 1 3)  correspond  to: 
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[Ql]  l|//f  =  5Tr (//*//). 

[Q2]  //+//  =  ////t  =  ||//||2I2. 

[Q3]  //  +  IP  =  2»(/ti)I2. 

[Q4]  II  =  IP  <=>•  //  =  J?(/i,)I2. 

2.3  Properties  of  Block  Quaternion  Matrices 


Here  we  consider  some  properties  of  vectors  and  matrices  with  entries  that  are  quaternions. 
For  example  a  2  x  2  quaternion  matrix  takes  the  form 


(2.16) 


where  si,s2,S3iS4  €  H  If  such  a  quaternion  matrix  is  Hermitian  symmetric,  i.e  S  =  SK 
we  refer  to  S  as  a  quaternion  symmetric  matrix  Such  a  matrix  has  real  quaternions  on  the 
diagonal. 

Theorem  2  A  quaternion  symmetric  matrix  S  has  the  following  properties: 

[BI ]  If  q i,  72*  *  *  •  i  <7.v  £  H  are  quaternions  then 


(q'i  92  •••  q’N)S(q\  92  9.v)r  =  a, 


(2.17) 


for  some  real  quaternion  a. 

IB2)  S  has  real  eigenvalues  and  complete  set  of  orthonormal  quaternion  eigenvectors. 

IB3]  If  S  is  nonsingular  then  S~  1  is  also  a  quaternion  symmetric  matrix. 

We  note  that  in  the  2  x  2  matrix  representation  defined  in  (2.14)  quaternion  matrices  are  com¬ 
plex  matrices  built  from  2x2  quaternion  blocks  of  the  form  (2. 1 5)  and  that  the  real  quaternions 
are  real  multiples  of  I2.  Since  the  represention  (2  14)  is  faithful  we  can  use  it  to  prove  algebraic 
relations  (as  we  do  below)  which  are  true  for  the  abstract  quaternions. 

Proof:  Property  B I  follows  from  the  fact  that  the  left  hand  side  of  (2.17)  evaluates  to  a 
quaternion  a  which  satisfies  a*  =  a  and  so  must  be  real 

Property  B2  follows  from  the  fact  that  S  is  a  complex  Hermitian  symmetric  matrix.  Define 
column  vectors  h  \  and  h 2  by 


(h\  h2)  =  (H\  f?2  •  //n)1  . 


(2.18) 


The  vectors  /ij  and  h2  are  orthogonal  and  satisfy  ||/i]  ||  =  ||Ji2||  We  then  have  by  Property 
B 1 ,  that  the  Rayleigh  quotient  is 


h\sh\  h^Shi 

hi  hi  h\li2 


(2.19) 
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Thus,  if  hi  is  an  eigenvector  of  S,  then  /i2  is  also  an  eigenvector  corresponding  to  the  same 
eigenvalue.  Overall  we  can  write  S  as 


N 

S  =  Y^HjKjK],  (2.20) 

J-l 

where  the  k3  are  the  real  eigenvalues  of  S  and  where  the  Kj  are  the  corresponding  normalized 
quaternion  eigenvectors.  The  quaternion  eigenvectors  satisfy  l\\K(  =  I2.  Finally,  if  S  is 

non-singular,  then  hy  (2.20),  1  is  block  quaternion  symmetric.  ■ 

We  note  that  the  relation  (2  19)  explains  the  computational  “miracle"  reported  in  [36]  where 
the  filter  that  estimates  c  is  independent  of  the  interfering  signal  s. 


2.4  Bayesian  Detection 

The  model  (2.5)  deserihed  above  implies  that  the  likelihood  function  of  codewords  c  and  s 
given  the  received  signal  r  is  given  hy 

p(r\c,s)  =  ^exl>(-^Hr-CW-^||2).  (2.21) 

Our  prior  knowledge  of  symbols  xt  is  that  they  are  selected  from  the  same  constellation  C 
independently,  and  that  each  symbol  in  C  is  transmitted  with  equal  probability.  We  will  also 
assume  that  the  constellation  is  centered  and  normalized  so  that 


(2.22) 


E{c}  =  0  and  E{c*c}  =  I2, 

E{s}  =  0  and  E{s^  s}  —  I2. 

From  Bayes'  rule  the  posterior  probability  of  c  and  s  given  the  received  data  r  is  given  by 

p(c)p{s)p{r\c  s)  _ 

:.  s) 

(2.23) 


p(c.  s|r)  = 


£c,.ec*  p{c)p(s)p(r\c.  s) 


Er  «€C2  p(rlc-s)' 

In  the  Bayesian  approach  to  decision  theory  we  introduce  a  loss  function  T(c,  s|c,  s)  which 
quantifies  the  loss  incurred  in  deciding  the  transmitted  symbols  are  c  and  s  when  c  and  s  were 
transmitted  [38]  Bayesian  decision  theory  then  seeks  to  minimize  the  expected  loss 


p(c,  s)  —  £(c,  s|c,  S)p(c,  s|r). 

c,sec2 


If  we  take  L  to  be  the  0-1  loss  function 

L(c,slc,s)  =  |j’  ' 


if  c  =  c  and  s  —  s, 
otherwise. 


(2.24) 


(2.25) 


then  p(c,  s)  is  minimized  by  the  maximum  a  posteriori  (MAP)  estimate 


(c.  s)  =  argma xp(c,  s|r), 

c,*CC2 


(2.26) 
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and  the  expected  loss  is 


(2.27) 


P(c,  s)  -  JZ  p(c,  s|r), 

which  is  the  posterior  probability  of  error  If  as  we  have  assumed  p(c)  —  1/|C|2  and  p(s)  — 
1/|C|2,  then  (2.26)  is  equivalent  to  the  maximum  likelihood  estimate 

(c,  s)  =  arg  maxp(r|c,  s).  (2.28) 

c,a£C2 


Let  x  —  (c  s),  A  —  (H  Q)1 ,  then  from  the  likelihood  function  p(r\c,s)  given  in  (2.21),  we 
have 

7>(r|c,  s)  a  exp  f  ~(x  -  x)/?  *(x  — x)M,  (2.29) 


where 


±  =  rA'(AA'y\ 


and 


R  1 


=  -rAT/t 


WIIHH2i2 
a2  V  GW' 


we*  \ 

lien2 1 J’ 


(2.30) 


(2.31) 


where  ||- 1|2  =  2  ||  ||p.  For  a  given  constellation  C,  the  performance  of  (2.28)  in  decoding  c 
and  s  is  determined  by  the  determinant  of  R~] .  We  have  from  [39], 


det(/?  ') 


l|H||4||g||4 

rr& 


(1-|A|2)2. 


(2.32) 


where 

HG' 

m  neir 


(2.33) 


The  parameter  A  is  an  inner  product  of  two  unit  quaternion  vectors,  and  measures  the  angle 
between  the  desired  signal  channel  vector  7i  and  the  interference  signal  channel  vector  Q. 
Wc  show  in  Section  5.3  [40]  that  this  parameter  is  fundamental  in  the  analysis  of  detection 
performance  for  multiple  Alamouti  schemes  involving  polarization. 

Clearly,  the  MAP  detector  in  (2.28)  involves  a  complex  search  for  a  large  constellation,  as 
(c,  s)  E  C4.  For  example,  a  24  sy  mbol  constellation  involves  216  evaluations  of  the  likelihood 
function  This  motivates  the  development  of  a  number  of  alternative  sub-optimal  solutions 
which  we  will  discuss  from  a  Bayesian  perspective. 


2.5  Bayesian  Interference  Cancellation 

The  motivation  behind  interference  cancellation  is  that  the  computation  of  full  MAP/ML  detec¬ 
tion  can  be  reduced  if  one  of  the  signals,  s  say,  can  be  canceled  out.  We  would  then  only  need 
to  search  over  c  E  C2 .  This  implies  that  we  should  attempt  to  marginalize  the  joint  posterior 
distribution  for  c  and  s  with  respect  to  s  This  leads  to 

p(c|r)  ex  ^  p(r\c  s)  (2.34) 

J  l  .Saco 
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and  the  MAP  decision  rule  is 


(2.35) 


c  =  arg  max  p(c|r). 

c€C2 

It  is  evident  that  this  does  not  help  as  the  sum  over  s  cannot  be  evaluated  analytically,  and 
so  the  evaluation  of  (2.34)  requires  just  as  many  likelihood  function  evaluations  as  (2.26)  or 
(2  28)  However,  we  notice  that  if  the  sum  in  (2.34)  were  to  be  replaced  by  a  Gaussian  integral 
the  marginalization  eould  be  computed  analytically.  In  Bayesian  terms,  the  trick  is  to  forget 
something  that  we  know  in  order  to  reduce  computation.  Instead  of  using  the  fact  that  s  lies  in 
some  constellation,  we  simply  recall  from  (2.22)  that 

E{s}  —  0  and  E{si  s}  —  I2.  (2.36) 


The  prior  for  s  is  then  taken  to  be  the  maximum  entropy  distribution  satisfying  the  constraints 
(2.22).  This  is  the  Gaussian  distribution  with  zero  mean  and  unit  variance,  i.e  p(s)  oc 

exp  )  The  new  prior  for  s  is  perfectly  consistent  with  our  partial  prior  knowledge 

of  s,  it  just  doesn't  represent  all  that  we  know.  Overall  this  means  we  allow  some  increase  in 
the  probability  of  error  in  delecting  codewords  c  (and  s),  for  reduced  computational  load.  We 
substitute  (2.21)  and  the  prior  p(s)  into  (2  26)  and  marginalize  the  posterior  distribution  with 
respect  to  s  to  obtain 


where 


P(c|r)  oc 


L 


/c 2 
(x  exp 


exp 


2  a  2 


|| r  —  cH  -  sQ || 2  )  exp 


H(c- 


c)rur  ln](c  -  c 


R~' 


h  _  &G_ 

"2  vH\m2  +  <r2)' 


(2.37) 


(2.38) 


and 


c^rR  1Ht(H/TlWt)‘1' 


(2.39) 


and  we  have  dropped  terms  independent  of  c. 

Since  R  1  is  a  quaternion  symmetric  matrix,  we  have  from  Theorem  2, 

TiR  lni  = /32h,  (2.40) 

where  (3 2  is  real  and  positive.  This  is  the  projeelion  of  interference  plus  noise  onto  the  signal 
subspaee  Substituting  (2.40)  into  (2.37)  we  obtain  two  independent  maximization  problems 
for  decoding  symbols  x\  and  x?  as 


where 


/  \  \  (  $c.  I  -  1 2  0c  I  -  (2^ 

p(c|r)  a  exp  I  -  —  \x\  -  xi\  -  —\xi  -  x2\  U 


(2.41) 


(2.42) 
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and 


_ !Ai!_  ) 

i  +  ^2/iienv 


(2.43) 


This  captures  the  solution  obtained  tn  [36],  [32],  and  will  be  called  the  1C  solution  for  conve¬ 
nience.  Note  that  if  Q  —  0,  then  the  performance  described  by  (2.43)  will  be  reduced  to  the 
performance  of  a  single  STBC  user  with  2  receive  antennas,  that  is 


&  ii^ii2  iiw.f  ,  HW2I12 

PlSTftC  *  ~  3  +  7  • 


(2.44) 


If  either  Hi  =  0  or  W.2  =  0  then  we  obtain  the  performance  of  a  single  STBC  user  with  one 
receive  antenna. 

Assuming  that  symbols  c  from  the  first  terminal  have  been  decoded  correctly,  the  receiver 
then  cancels  the  contribution  of  the  first  terminal  in  the  received  signal  vector  r.  The  receiver 
decodes  symbol  s  by  applying  MAP  decoding  (2.28)  to  the  received  signal  vector  after  cancel¬ 
ing  signals  from  the  first  terminal.  This  gives 


p(s\r')  oc  exp 


(2.45) 


where 


r  =  r  -  cH ,  (.T3,  X4) 


1 

W 


r'g] 


and 


(2.46) 


Assuming  that  the  symbols  n  and  xi  have  been  decoded  correctly,  we  can  see  from  (2.46) 
that  the  performance  for  the  decoded  symbols  s  is  equivalent  to  that  of  1  user  with  2  transmit 
and  2  receive  antennas  as  previously  noted  in  [32].  The  performance  for  decoding  symbols  c 
and  s  is  given  by 


,2  _  iiwn4  lien4 
a'ct  -  73 


v  l+^/iieiiv 


(2.47) 


At  this  point  we  note  that  it  is  clearly  advantageous  to  have  chosen  to  cancel  the  Alamouti 
channel  with  the  lower  SNR  Thus,  if  ||(J||  >  ||7Y||,  we  would  exchange  the  roles  of  H  and 
Q  in  the  above  algorithm  and  analysis.  We  further  note  that  x\  and  X2  are  assumed  to  have 
been  decoded  correctly,  therefore  (2.47)  provides  an  upper  bound  for  the  performance  of  the 
IC  solutions 

Alternatively,  symbols  s  could  be  decoded  in  a  similar  fashion  to  the  decoding  of  symbols 
c.  In  this  case  the  posterior  probability  for  s  given  r  becomes 


p(s|r)  (x  exp 


yN-*3|2 


(2.48) 


where 


and 


« =  4 irR  "s1. 

Pi 

=  mb!  (\  iai2  n 

V  \+CT*I\\H\\2)' 

r ,  =  I4  H'H  _ 

^  '  <T*(\\n\\2+o*y 


(2.49) 

(2.50) 

(2.51) 
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The  performance  is  then  given  by 


02  =  MW 

<7S 


(i _ _ !*! !_  _V 

v  l  +  ^/liaii2/  V  i  +  ff2/liwiiv 


(2.52) 


2.6  The  Zero-Forcing  Decoder 


Naguib  and  Scshadri  [37]  show  that  if 

W  = 

then 


-WH\ 


9t\ 

"-(«)- (5) -(??)©+ (5)- 


where 


W  =  Ht  -  QiGi'n 2  and  5'  =  Q2  -  HiHT'Gi. 


(2.53) 

(2.54) 

(2.55) 


The  matrix  W  transforms  the  problem  of  joint  deteetion  of  two  co-channel  users  into  separate 
detection  of  two  individual  space-time  users.  Furthermore,  the  algebraic  structure  of  the  block 
space-time  code  (closure  under  addition,  multiplication  and  taking  inverses)  implies  that  the 
matrices  H!  and  Q'  have  the  same  structure  as  the  matrices  7i\,7i2,  Q\  and  £2-  Now  the  new 
noise  vectors  n\  and  n2  are  correlated,  and  if  we  take  these  correlations  into  account  and 
design  the  optimal  detector  we  would  just  be  back  to  (2  26).  The  key  to  reducing  computation 
is  to  again  forget  something  that  we  know  In  this  ease  we  forget  that  the  noise  veetors  n[  and 
n2  are  correlated  and  replace  the  joint  distribution  of  n\  and  n2  by  the  higher  entropy  product 
of  marginals 

This  reduces  the  problem  of  joint  deteetion  of  codewords  c  and  s  to  four  independent  max¬ 
imization  problems.  That  is  we  have  the  likelihood  of  the  codew  ords  given  the  received  signal 


where 


with 


P(ri  |c)  a  exp  y|xi  -  ill2  -  y  |x2  -  x2|2  j  , 
/)(r'2|s)  a  exp  |x3  -  x3|2  -  y  |x4  -  x4|2^  , 

r'.W't  rtf'  ^ 


(c,s)  = 


l|W||2(l  -|AP)’||C?||2(1  -|AP), 


and 


The  performance  of  the  7ero-forcing  decoder  for  symbols  c  and  s  is  given  by 


(2.56) 

(2.57) 


(2.58) 


(2.59) 


A- 1 «(!-!# 


(2.60) 


22 


Figure  2.  I:  Behaviour  of  fC  for  two  users  of  Alarnouti  code ,  at  different  S/NR 


2.7  Performance  Analysis 

First  note  that  the  expressions  in  (2.32),  (2,47),  (2.52)  and  (2.60)  give  the  effective  SNR  when 
estimating  c  and  s.  This  effective  SNR  includes  both  noise  and  interference  Now  consider 
the  performance  expression  described  in  (2  43)  or  (2.50),  if  the  power  of  the  interfering  path  is 
small  i.e.  if  ||£||2  — *  0  then  ft2  — ►  \\H\\2  /a2.  That  is,  the  performance  improves  as  the  inter¬ 
ference  power  decreases  and  converges  to  that  of  a  single  space-time  user  and  a  receiver  with 
two  antennas  [36].  Figure  2.1  shows  the  simulation  results  which  demonstrate  this  behavior. 

To  compare  the  performance  of  the  three  different  solutions,  we  first  determine  the  effective 
SNR  of  a  symbol  for  each  method  by  taking  a  fourth  order  root  in  (2.32),(2.6Q)  and  (2  52),  to 
obtain 


SNRMai> 

•SNR/,: 


SNRir 


If  the  channels  H  and  Q  are  orthogonal  i.e.  |A|2  =  0  then  there  is  no  difference  in  the  perfor¬ 
mance  of  the  three  solutions.  From  (2.61  )-(2  63),  it  is  clear  that  when  |A|2  ^  0,  the  maximum 
likelihood  solution  outperforms  the  zero-forcing  and  1C  solutions.  Furthermore,  the  perfor¬ 
mance  gap  (the  ratio  or  difference  in  dB  for  a  given  error  rate)  between  zero-forcing  ,  1C 


ll^ll  11011, ,  |  \  ,2\ 


<7* 

iiwii  lian 


L(1  -  |A|"), 


l(l  ^  )  (j  lAl2  ^ 
VV  1  +  ^7  nail2  A  i  +  a2/ \\h\\2 / 


(2.61) 

(2.62) 

(2.63) 
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Figure  2.2  Performance  comparison  of  the  three  different  detection  algorithms  for  two  users 
of  Alamouti  code 


solution  and  maximum  likelihood  solution  can  be  explicitly  described  as  follows: 


SNR/). 

SNRmai1 


SNRic 

SNRmai* 


V  1  -  I'M2- 


\ 


(\  -  i*i2  , \  (\  _  w*  'i 

V  HgTWA1 


(1  -  |A|2) 


(2.64) 

(2.65) 


Figure  2.2  shows  the  simulation  results  of  interference  cancellation  for  two  co-channel  users, 
each  using  the  Alamouti  scheme  with  8PSK  modulation  comparing  the  three  different  solu¬ 
tions.  The  channels  Ft  and  Q  are  taken  to  have  the  channel  separation  | A| 2  =  0.8  with  a 
tolerance  of  0.05,  or  in  other  words,  the  angle  between  the  two  channels  is  ^  37°.  In  order  to 
show  that  the  performance  expressions  obtained  above  arc  valid.  Figure  2.3  shows  the  result  of 
shifting  the  zero-forcing  and  IC  solutions  shown  in  Figure  2.2  by  the  performance  gap  given  in 
(2.64)  and  (2.65)  respectively  As  shown  all  curves  arc  now  aligned. 
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Figure  2.3.  Performance  of  the  three  different  detection  algorithms  agree  after  the  SNR  is 
shifted  by  the  differences  given  in  (2.64)  and  (2.65) 


2.8  Generalization  to  N  Alamouti  Users 

We  now  extend  the  analysis  of  Section  2  to  detection  of  N  co-channel  users  using  N  receive 
antennas,  where  eaeh  user  employs  the  Alamouti  STBC.  The  generalization  of  (2.5)  is  straight¬ 
forward,  and  the  reeeived  signal  ean  be  written  as 

N 

r  =  x„G„  +  n  (2.66) 

n  —  1 

where  Gn  =  (Gn  i  ,  Gn. 2,  •  •  •  ,  £n,/v)  are  the  channel  matrices  from  user  n  to  receiver  antennas 
1, 2,  •  •  •  ,  N  and  x»  =  (x„.i  xn,2)  is  the  codeword  transmitted  by  user  n. 


2.8.1  MAP/ML  Detection 

Let  x  =  (x i ,  X2 ,  *  *  *  ,xj v)  and  A  =  (Gi,G2,  •••  ,G^)!  then  the  posterior  distribution 
p(x i,  --  ,  Xiv|r)  is  given  by 

p(xi,  •  •  •  ,  xjv|r)  oc  exp  ^-(x  —  x)R  1  (x  -  x)^  .  (2.67) 

where  R~x  —  AA*  and  is  obtained  by 


/||Gi  ||2 12 

GlG^ 

GiG|v  \ 

R 1  =  -L 

g2g{ 

l|G2||2l2  ■ 

GiG^ 

2a2 

\  GVG| 

G  \  G2 

•  ||G^||2I2/ 
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The  performance  of  the  MAP/ML  detector  is  determined  by 


0ml  =  flat  R  ' 


I|G1||4||G2||4---||Ga11i,2 

a*N 


(2.69) 


where 


with 


1/2 


/  12  Ai,2  ■  ■  . 

a*,2  i2  ... 

Ai.at\ 

A  2,N 

det 

\Aj  n  a2 n  . . . 

I2  / 

\ 

GmGj, 

11  l|Gm||  ||G„|| 

(2.70) 


(2.71) 


Writing  {gn  \  ,  gn2)  =  Gn/  ||G„||  we  see  that  the  quantity  V  is  the  volume  of  the  parallclotope 
generated  by  the  complex  vectors  {gri  x  gn  2\n  =  1  ■  •  ,  Ar}.  It  is  a  function  of  the  pairwise 

quaternion  angles  between  the  users.  Since  the  vectors  gn  ,  are  normalized,  0  <  V  <  1.  We 
refer  to  V  as  the  normalized  channel  volume  In  terms  of  V\  we  can  define  the  signal  to  noise 
ratio  per  user  symbol  as  the  (2Ar)*h-root  of  (2  69),  that  is, 


SNR  ml  =  0(SNRi,SNR2,"‘  ,SNRjV)  V]/N ,  (2.72) 

where  SNRn  =  ||G„  ||2  /cr2,  is  the  SNR  of  the  nth  user  and  Q  denotes  the  geometric  mean.  For 
the  two  user  case  described  in  Section  2.5,  the  normalized  channel  volume  V  =  1  -  |A|2. 


2.8.2  Bayesian  Interference  Cancellation 

Bayesian  interference  cancellation  in  the  N  user  case  is  analogous  to  the  two  user  case,  al¬ 
though  it  becomes  tedious  to  marginalize  the  posterior  probability  p(x \,X2,  ■  ■  •  ,  x/v|r)  with 
respect  to  x  because  of  the  multiple  integrations  However,  since  we  assume  the  prior  distribu¬ 
tion  for  X2 ,  £3,  *  *  *  » &i\  lo  be  Gaussian  with  zero  mean  and  unit  variance,  we  have 


Af(0,  fl), 


where 


N 

ll  =  4  N  ^  Gp  Gn. 

n  2 


Thus,  we  can  write  down  the  posterior  distribution  p(x\  |r)  directly  as 


p(x\ |r)  a  exp 


aT.i 


(2.73) 


(2.74) 


(2.75) 


where 

5i  =  Trfl->Gi.  (2.76) 

and  0  will  be  computed  below.  Note  that  the  reduced  MAP  decoding  now  consists  of  two 
independent  maximization  problems  for  decoding  symbols  S|,i  and  s i.2. 
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In  order  to  judge  the  performance  of  this  reduced  MAP  decoding  scheme  we  expand  the 
positive  semi-definite  block  quaternion  symmetric  covariance  matrix  /?'  given  by 


N 

n'  =  Y,  G»G"  (2-77> 

n  — 2 

in  terms  of  its  block  quaternion  eigenvectors,  (See  Appendix),  so  that 


X 

R  =  cr2l2N  +  £kJ;k„,  (2.78) 

»i=  2 

where  K^Kj  =  0,  if  n  ^  j ,  and  we  have  normalized  the  eigenvectors  KT1  so  that  ||/\n  ||  2  =  k„, 
their  corresponding  eigenvalues  In  general  some  of  the  k}  could  be  zero  and  so  they  would  not 
contribute.  Thus,  we  have  replaced  the  Alamouti  signaling  interferes  Gn  hy  an  equivalent  set 
of  mutually  orthogonal  Alamouti  signaling  interferers  K„.  Interference  cancellation  can  now 
be  analyzed  in  terms  of  these  equivalent  interferers. 

Applying  the  Matrix  Inversion  Lemma  [41],  the  inverse  of  R  is  seen  to  be 


and  so  we  have 


where 


R~l 

02  = 


1  l  N 


k„kJ, 


G, 


>-E 


|A„|2 


1  +<r2/l|K„||‘ 


A,.  = 


Gi  Kj, 

iiGiii  iiK„ir 


(2.79) 

(2.80) 

(2.81) 


|  An  | 2  measures  the  angle  hetween  the  desired  signal  channel  Gi  and  the  nlh  equivalent  mtcr- 
ferer  channel  Kn. 

As  with  the  1C  algorithm  descrihed  in  Section  2.5,  interference  cancellation  is  first  earned 
out  to  separate  the  user  G 3  with  the  highest  signal  to  noise  ratio.  The  transmuted  symhols  for 
this  user  are  estimated  (detected)  and  then  substituted  in  (2.66)  to  ohtain 


x 

r—r-  SjGj  =  ^  sn Gn  +  n  (2.82) 

n  1 
n*J 


The  above  procedure  is  then  repealed  with  the  Alamouti  channel  having  the  next  highest  SNR 
and  so  on.  This  algorithm  is  the  counterpart  of  Foschim’s  scheme  [42]  for  space-time  multi 
plexing  one  dimensional  signals.  Here  we  are  multiplexing  2x2  matrix  signals. 


2.9  Summary 

Our  Bayesian  analysis  of  interference  cancellation  provides  a  unified  framework  with  which 
to  understand  the  relationship  between  various  techniques  for  signal  detection.  We  ohtained 
simple  new  results  for  performance  analysis  in  terms  of  SNR.  The  analysis  also  leads  to  a 
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new  parameter  which  measures  the  angle  between  a  pair  of  block  quaternion  channel  vectors. 
This  parameter  provides  new  theoretical  insight  to  understanding  and  being  able  to  predict 
the  performance  of  different  decoding  algorithms  as  a  function  of  signal  to  noise  ratio.  In 
Chapter  we  apply  this  parameter  to  analyse  the  performance  of  various  coding  schemes  for 
dual  polanmetric  antennas. 
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Chapter  3:  Conditional  Optimization  and  Fast  Decoding 
Algorithms  for  Space-Time  Block  Codes 


The  drawback  of  the  full-rate  full-diversity  STBCs  listed  in  Section  1.5  is  the  decoding  com¬ 
plexity.  The  sphere  decoder  was  developed  [43,44,45]  to  reduce  the  computational  complexity 
by  limiting  the  search  to  only  within  a  sphere  of  given  radius.  However,  the  initial  radius  of 
the  sphere  depends  on  the  signal-to-noisc  ratio  and  the  channel  condition  number.  When  the 
channel  matrix  is  close  to  singular,  the  preprocessing  stage  of  the  sphere  decoding  algorithm 
yields  a  plane  of  possibilities  rather  than  a  single  initial  estimate.  When  this  occurs,  the  lattice 
point  search  degenerates  to  an  exhaustive  search,  and  consequently  the  overall  complexity  of 
the  sphere  decoder  is  no  belter  than  the  exhaustive  search.  In  wireless  communication,  when 
the  channel  between  base  station  and  terminal  is  line-of-sight,  the  induced  channel  is  rank  1 
This  has  motivated  us  to  develop  fast  decoding  algorithms  for  a  number  of  important  codes 
which  have  fixed  low  complexity  across  channel  conditions,  i  e.  across  propagation  conditions 
from  pure  line-of-sight  to  rich  scattering  environments. 

Our  fast  decoding  algorithms  are  based  on  a  technique  called  conditional  optimization  This 
technique  is  used  widely  in  statistical  estimation  and  signal  processing  to  reduce  the  search 
space  of  an  optimization  problem  by  taking  advantage  of  the  possibility  of  analytically  op¬ 
timising  over  some  subset  of  the  parameters  conditioned  on  the  remaining  parameters.  The 
approach  is  applied  to  a  large  class  of  full  rate,  full  diversity  space-time  block  codes.  This  ap¬ 
proach  can  also  be  applied  in  an  approximate  form  to  the  Golden  Code  and  perfect  space-time 
block  codes  to  obtain  essentially  ML  performance  with  reduced  complexity  [46,  47]  compared 
to  exact  ML  decoding. 

3.1  Two  Users  of  Alamouti  Signaling 

In  Chapter  2  we  considered  an  interference  cancellation  approach  to  the  detection  of  two  users 
of  Alamouti  codes,  which  provides  linear  complexity  processing  with  a  corresponding  trade¬ 
off  in  performance.  We  show  here  that  such  systems  have  the  structure  that  allows  fast  ML 
decoding  with  complexity  0(N 2),  wherre  N  is  the  size  of  the  underlying  constellation. 

From  (2.3),  the  received  signal  for  a  two  Alamouti  user  system  is  given  by 


r  i  -  cH\  +  sQ i  +  ri| 
T2  —  cJi.2  +  3^2  +  TI2 


(3.1) 


where  c  =  =  (2:3, X4), 


(3.2) 


Interference  cancellation  as  described  in  Chapter  2  is  a  suboptimal  approach  with  linear  com 
plexity,  [48,  37].  Here  we  show  that  two-user  ML  detection  of  Alamouti  signals  is  possible 
with  complexity  C(N2). 
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3.2  Code  Multiplexed  Alamouti  Blocks 


An  approach  to  constructing  high  data  rate  codes  is  to  multiplex  existing  codes.  For  example, 
a  high  rate  space-time  code  might  he  constructed  as  a  "linear”  combination  of  Alamouti  blocks 
[27,  30,  31, 49]  of  the  general  form 


X  = 


(3.3) 


where  the  matrix 


is  usually  taken  to  be  a  unitary'  matrix  [50].  This  restriction  on  A  leads  to  code  properties  such 
as  eubic  shaping  and  information  losslessness  (see  [7]  for  details  on  eubie  shaping)  However, 
we  note  our  fast  decoding  algorithms  do  not  depend  on  this  property  and  will  work  for  arbitrary 
non-singular  A. 


Remark  1  If  A  is  a  unitary  matrix  of  the  form 


A  = 


(3.4) 


then  (3.3)  reduces  to  Alamouti  signaling  with  a  non-standard  constellation  and  decoding  is 
essentially  linear  in  complexity.  One  example  is  Alamouti  signaling  with  16-QAM  viewed  as 
a  multiplexing  of  Alamouti  blocks  employing  4-QAM.  However  a  disadvantage  of  this  con¬ 
struction  is  that  the  value  of  adding  a  second  receive  antenna  does  not  extend  beyond  noise 
averaging,  in  contrast  to  codes  such  as  the  Golden  Code  where  it  provides  additional  diversity 
gain. 


A  number  of  codes  proposed  in  the  literature  are  of  this  multiplexed  form.  We  discuss  these 
examples  in  turn. 

Hottinen-Tirkkonen-Kashaev  (49,  27]  These  authors  considered  the  family  of  2  x  2  full- 
rate,  full-diversity  STBCs  as  described  in  (1.52).  This  code  ean  be  expressed  in  the  form  of 
(3.3)  with 


1  /  lTi  -1  +2i\ 
y/7  V-1  ~2t  ~1  +  i  J  ' 


(3.5) 


At  the  receiver,  the  received  signal  is  represented  as  in  (3.1)  where  Hi,  H2  are  matrices  as 
defined  in  (3.2),  and 


r  —  (S11  ~oh\  r  (q\2  —922\ 

Gl  ~  U>  Six)'  G2-\922  9h)  (16) 

where 

(#ll»#2l)  =  (/*11>/*2l)A,  (#12,422)  =  (h .12,  /l22)A.  (3.7) 

Sezgincr-Sari  [30]  Sezginer  and  Sari  proposed  a  2  x  2  full-rate,  full-diversity  STBC  given  in 
(1 .55)  which  can  be  written  in  the  form  of  (3.3), 


X  = 


/a.ri 


dx3~\ 
-<1*1 )’ 


(3.8) 
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and  also  leads  to  the  received  signal  expressions  as  in  (3. 1 )  with 


H  i  - 
«a  = 


(3.9) 


Rabiei-Al-I)hahir  [31]  Rabiei  and  Al-Dhahir  propose  the  code  given  in  (1.58)  which  can  be 
rewritten  as 

_  /  Tiai  x4ai\  / 0  1\  (  x3a2 

\-X402  x\02)  +  yl  0 J  y-x^i 


^2Of2\ 

*301 )  ’ 


(3.10) 


and  the  received  signal  can  also  be  expressed  as  in  (3.1)  with 


c  =  (xi.xj),  s  =  (x3,  X2  ) , 
r  1  =  (ru,rl2),  r2  =  (r2i,r22)< 


n2 


(  hii^i  h^i02\  r  _  (  h2i<*2 

\-fi2\ih  'yi  ~  \-hu0i  /i-2ia2y’ 

(  hl2&l  ^2202\  q  f  h.22<*2  hy20l\ 
\-h 2202  fihcki)  ’ y2  ~  \-hi20l  /l22<*2/ 


(3.11) 


3.3  Fast  Optimal  Decoding 

As  shown  in  the  last  section,  the  decoding  problem  for  all  cases  takes  the  form 

r  1  =  cH\  +  sQ\  -f  rii 

(3.12) 

r2  =  cH2  +  sQ  2  + 

where  c  and  s  are  vectors  of  a  pair  of  symbols  x,  and  H\ ,  H2,  Q\  and  Q2  are  Alamouti  blocks. 
Assume  that  symbols  x,,  z  =  1  ...  4  are  selected  from  some  QAM  constellation  C  of  size 

N.  We  now  show  that  exact  ML  can  be  implemented  with  complexity  0(N2),  through  a  very 
simple  algorithm  hirst  write  (3. 12)  as 

r  -  cH  +  sQ  4-  n  (3.13) 

where  r  =  (ri,r2),  H  =  (H1/H2),  Q  =  {G\<G2)  and  ti  =  (m,n2).  The  likelihood 
function  for  (3. 1 3)  is 

;>(r|s,c)  ocexp  ^  9^2  Hr  cH  ~  sG\\2^j  (3.14) 

and  the  ML  estimate  is  given  by 


(c,  s)  =  arg  max  p(r\s.  c).  (3.15) 

c,s€C2 

At  first  glance  this  optimization  seems  to  involve  a  search  over  C4.  However,  conditional 
optimization  can  be  applied  to  optimize  (3.14)  exactly  as  follows.  First  maximize  (3.14)  with 
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respect  to  c  given  s,  by  expanding  (3.14)  and  completing  the  square,  we  obtain 


p(r |s,  c)  oc  exp  (  -  ^  (r  -  sQ)  (  I4  - 

jjgf 

2<T2 


H*H  \ 

IIWIIV 


(r  -  st?) 


') 


x  exp  I  — 


|c  -  c(s) 


where 


c(s)  = 


rH1  -  sQV} 

IIWII2 


and  m||2  =  5  mi|2-.  Thus  the  ML  estimate  of  c  given  s  is 

c(s)  =  Q(c(s)) . 

Substituting  (3.18)  into  (3. 14),  we  obtain  the  optimization  problem  for  s  alone 


(3.16) 


(3.17) 


(3.18) 


s  =  argmax  \\r  -  c(s)H  -  sQ\\  . 
sec2 


(3.19) 


Thus,  (3.19)  and  (3  18)  provide  an  algorithm  for  obtaining  the  ML  solution  of  c  and  s,  which 
involves  at  most  |C|2  =  N2  evaluations  of  the  likelihood  function.  That  is,  we  have  an  0(N2) 
algorithm  for  estimating  s  and  c. 

Remark  2  Above  we  have  assumed  that  a  QAM  constellation  is  used  This  is  not  necessary 
for  conditional  optimization  to  provide  a  computational  benefit ,  but  does ,  along  with  certain 
other  lattice  based  constellations,  provide  maximum  computational  benefit.  For  an  arbitrary 
constellation  the  quantization  step  (3.18)  is  replaced  by  a  search  which  is  at  most  0(iC|), 
while  if  the  constellation  is  a  Cartesian  product  of  two  real  constellations,  i.e  C  —  TZ  x  7Z 
then  this  search  is  at  most  0(  \/TC|)-  Finally,  if  the  constellation  is  a  subset  of  one  of  a  number 
of  lattices,  such  as  the  QAM  or  HEX  constellations,  then  the  quantization  step  is  0(1). 


3.4  Conditional  Optimization 

In  the  previous  section  we  used  a  conditional  optimization  approach  to  develop  fast  ML  de¬ 
coding  algorithms  for  multiplexed  Alamouti  codes.  Here  we  consider  conditional  optimization 
from  a  more  general  perspective  and  discuss  its  application  to  fast  decoding  of  space-time 
codes. 

Conditional  optimization  is  a  technique  widely  used  in  statistical  estimation  and  signal  pro¬ 
cessing.  The  goal  of  this  technique  is  to  reduce  the  search  space  of  the  optimization  hy  taking 
advantages  of  the  possibility  of  analytically  optimizing  over  some  suhset  of  the  parameters  con¬ 
ditioned  on  the  remaining  parameters,  or  at  least  more  efficiently  than  an  exhaustive  search  An 
archetypal  example  of  this  approach  is  Rife  and  Boorstyn’s  [51]  This  reduces  the  parameter 
estimation  problem  for  a  single  tone  in  noise,  which  involves  the  estimation  of  three  parame¬ 
ters,  to  an  optimization  prohlem  involving  only  the  frequency  of  the  tone. 

In  general,  suppose  that  we  wish  to  maximize  a  likelihood  function  of  the  form 

p(r|0)  (3.20) 
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where  r  is  some  data  and  0  is  a  set  of  parameters  we  need  to  optimize  over  If  the  parameter  set 
can  be  split  0  =  (0i,  02)  such  that  optimization  overtfi  given  02  can  be  carried  out  analytically 
or  at  least  very  efficiently,  then  the  optimization  problem  can  be  efficiently  carried  out  as: 

02  =  argmaxp(r|02,0i(02)),  (3.21) 

*2 

where 

0i (02 )  =  argmaxp(r|02,0i),  (3.22) 

0i 

and  —  0i  (02). 

We  now  consider  sufficient  conditions  under  which  conditional  optimization  leads  to  a  re¬ 
duction  in  decoding  complexity  of  space-time  codes.  Assuming  that  perfect  channel  state 
information  is  available  at  the  receiver,  the  received  signal  is  given  by 

Y  =  HX  +  W  (3.23) 

where  X  is  the  transmitted  STBC  codeword  with  code  length  /  and  entries  that  are  information 
symbols  drawn  from  an  A'-QAM  constellation,  H  is  the  matrix  of  channel  gains  from  the 
transmit  to  the  receive  antennas  and  n  is  i  i.d  Gaussian  noise  with  zero  mean  and  covariance 
2cr2I  „,T. 

The  received  signal  can  be  written  as 


r  —  xTi  +  n  (3.24) 

where  x  is  the  transmuted  information  symbol  vector  and  H  is  the  induced  (equivalent)  channel 
matrix 

The  likelihood  function  of  symbols  x  given  the  received  signal  r  is  given  by 

;j(r|x)  oc  exp  ~  ||r  -  xH\\2^j  .  (3.25) 

Taking  the  prior  distribution  of  symbols  x  to  be  uniform  on  the  constellation  Cs,  we  obtain 
the  ML  estimate 

x  =  arg  maxp(r|x).  (3.26) 

ML  decoding  is  certainly  achieved  with  |C|A  computations  of  likelihood  (3.25),  but  if  the 
symbols  are  taken  from  an  Af-QAM  constellation  then  dramatic  reductions  in  complexity  are 
possible. 

Theorem  3  If  the  induced  channel  matrix  7i  has  m  rows  which  are  mutually  orthogonal,  for 
all  channels  H,  then  exact  ML  decoding  can  he  implemented  with  complexity  0(Nh  ~Tn). 

Remark  3  It  follows  that  in  all  cases  the  complexity  of  exact  ML  decoding  is  at  most  0(N ^  1 ). 

Proof  Let  H  i  be  the  sub  matrix  consisting  of  m  mutually  orthogonal  rows  of  the  induced 
channel  matrix  7i  and  let  H2  be  the  matrix  consisting  of  the  remaining  rows.  We  can  then 
rewrite  (3.24)  as 

r  =  XiLii  +  xfHi  +  n  (3.27) 
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where  X\  is  the  1  x  m  vector  of  complex  symbols  corresponding  to  the  rows  in  H\  and  X2 
consists  of  the  h  in  remaining  symbols.  The  likelihood  function  associated  with  (3.27)  is 

p(r|xi,xa)  nc  exp  ^“^2  Hr  “  x*Ui  “  •  (3.28) 

We  have  isolated  xi,  as  our  approach  is  to  maximize  the  likelihood  function  (3.28)  with  respect 
to  x\  given  x2,  and  then  maximize  the  resulting  partially  optimized  likelihood  function  with 
respect  to  x2.  Now  since  the  rows  of  Hi  are  mutually  orthogonal 

H\H\  =  D,  (3.29) 

where  D  =  diag(df,  •  •  •  .  </^„ ).  and  dj  is  the  norm  of  the  j'h  row  of  Hi.  The  likelihood 
function  (3.28)  can  then  be  written  as 


p(r|xi,  ij)  (X  exp  (- -  *2X2)  (l-W{l>  ’Wi)  (r  -  X2H2)' 


PXp  (  “2^2  -M*2)|: 


j-i 


where 

Thus,  given  x2, 


x\(xz)  =  (r  x2H2)H{D  l. 


(3.30) 


(3.31) 


(3.32) 


X\(x2)  -  (Q(xi(x2)),-  •  ,  Q(xm(x2))), 

=  Q(xi(x2)), 

maximizes  (3.30)  Substituting  (3.32)  into  (3.30)  wc  obtain  the  optimization  problem  for  x2: 


x2  =  argmin  \\r  -  X\ (x2)Hi  —  x2H2||2 


(3.33) 


which  can  be  substituted  back  into  (3.32)  to  obtain  the  ML  estimate  of  xj 


x  i  =  Xi(x2).  (3.34) 

Thus,  equations  (3.33)  and  (3.34)  provide  an  algorithm  for  obtaining  the  ML  estimate  ofx*,  i  = 
1, . . . ,  K%  which  involves  at  most  Nh  ~m  evaluations  of  the  right  hand  side  of  (3.33).  ■ 

Wc  see  that  for  codes  satisfying  the  conditions  of  Theorem  3,  conditional  optimization  re¬ 
duces  the  search  for  the  ML  estimate  to  optimization  of  the  objective  function  in  (3.33)  over  the 
symbols  in  x2  This  reduced  search  can  be  carried  out  by  any  means,  in  parallel,  sequentially, 
by  tree  search  and/or  using  sphere  decoding.  Wc  will  further  elucidate  its  relation  to  sphere 
decoding  in  Section  3.5.  For  non  square  QAM  see  Remark  2 

Remark  4  If  X\  and  x2  are  vectors  of  real  symbols,  then  the  above  analysis  holds  with  the 
following  modifications.  Firstly,  in  Theorem  3  “mutual  orthogonality  "  of  the  rows  of  Ti  (see 
Equation  (3  29))  is  replaced  by  the  condition 

HiH\+HYH\  =  D.  (3.35) 

Secondly ,  Equation  (3.31)  is  replaced  by 

x,(*2)  =  (rH\  +  rH\  -  x2(H2H\  +  H2H\)\  D~l.  (3.36) 

(Here  denotes  complex  conjugate.) 
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3.5  Relationship  to  Sphere  Decoding 

The  idea  of  sphere  decoding  is  to  search  only  within  a  sphere  or  ellipsoid  of  a  certain  size 
defined  by  the  induced  covariance  matrix  and  the  SNR.  Suppose  for  a  given  decoding  we  have 
determined  that  wc  need  only  search  symbols  that  lie  within  the  ellipsoidal  region  S  C  C\ 
defined  by  the  equation 

5  :  <  py  (3.37) 

where  R  is  a  diagonal  positive  definite  Hermitian  (symmetric,  if  x  is  real)  matrix  and  p  >  0. 
Having  decided  on  a  sphere  the  reduced  decoding  problem  becomes 

x  =  arg  maxp(r|x).  (3.38) 

x€Ck ns 

In  Theorem  3  we  split  the  code  vector  x  into  two  parts,  x\  associated  with  the  ni  mutually 
orthogonal  rows  of  H ,  and  X2  associated  with  the  other  K  -  m  rows  of  H.  Let  x  be  a  general 
vector  in  CK  and  define  U2  :  Ch  — >  CK~m  by 

xU 2  -  X2,  (3.39) 

so  that  II2  =  is  the  orthogonal  projector  on  the  “X2”  subspacc.  Similarly,  we  can  define 

an  orthogonal  projection  II 1  =  V\ U\  onto  the  “xi"  subspace.  If  x  6  5,  then  x2  is  contained 
in  the  region 

Sn2  •  ^2  ^2 ^  *^2)  x\  <  p,  (3.40) 

where  (/^R  XU2  €  C{K~m)x{h  ~m). 

The  optimization  problem  for  X2  thus  reduces  to 

X2  =  arg  mm  ||r  xi(x2)7fi  x27f2||2.  (3.41) 

x2€Ch  ~  ,nnS n2 
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The  situation  is  illustrated  in  Fig.  3.1,  where  the  axes  laheled  x\  and  x2  represent  the  two 
suhspaces  dehned  by  the  orthogonal  projections  II i  and  02.  The  grey  points  represent  the 
projection  of  the  full  constellation  Ch  onto  the  II2  subspace.  The  projected  constellation  is 
searched  over  the  projection  of*  the  sphere  S  onto  the  subspace  n2,  that  is  5n2. 

Clearly,  the  method  of  conditional  optimization  can  he  naturally  integrated  with  sphere  de¬ 
coding  in  a  simple  way  leading  to  a  sphere  decoding  algorithm  with  worst  case  complexity 
corresponding  exactly  to  that  of  our  direct  algorithm  in  (3.33)-(3.34),  that  is  0(\C\K~m). 


3.6  Fast  Decoding  of  the  Golden  Code 

It  has  been  shown  recently  that  the  ML  decoding  complexity  of  the  Golden  Code  is  cubic  in 
the  size  of  the  underlying  QAM  constellation  for  an  arbitrary  iV-QAM  constellation  [46,  52]. 
In  Chapter  4  we  show  that  if  a  square  QAM  constellation  is  used  then  the  complexity  of  exact 
ML  decoding  is  further  reduced  to  Q(N2 y/N).  Here  we  apply  the  conditional  optimization 
approach  to  the  Golden  Code  to  obtain  a  simple  approximate  quadratic  complexity  decoding 
algorithm  with  essentially  ML  performance.  For  a  large  constellation,  this  is  a  significant  re¬ 
duction  in  complexity  when  compared  to  0(N2y/N).  The  algorithm  can  be  employed  by 
mobile  terminals  with  either  one  or  two  receive  antennas  and  it  is  resilient  to  near  singular¬ 
ity  of  the  channel  matrix.  Dual  use  is  an  advantage,  since  there  will  likely  be  some  IEEE 
802.16  mohile  terminals  with  one  receive  antenna  and  some  with  two  antennas  The  key  to 
the  quadratic  algorithm  is  a  maximization  of  the  likelihood  function  with  respect  to  one  of  the 
pair  of  signal  points  conditioned  on  the  other  This  choice  is  made  by  comparing  the  determi¬ 
nants  of  two  covariance  matrices  (or  equivalently,  the  norm  of  two  channel  vectors),  and  the 
underlying  geometry  of  the  Golden  Code  guarantees  that  one  of  these  choices  is  good  with 
high  probability. 


3.6.1  Model  and  Decoding  Problem 

The  Golden  Code  is  a  2  x  2  block  space-time  code  that  employs  2  transmit  and  2  receive  an¬ 
tennas  and  encodes  four  complex  QAM  .symhols  over  two  time  slots  yet  achieves  full  diversity 
[22,  53,  54]  and  is  incorporated  in  the  IEEE  802. 16  standard.  As  described  in  Section  1.5.1  the 
codewords  in  the  Golden  Code  lake  the  form 


/a  0\  /  x\  +  TX2  :r3  +  t: r4\ 

yO  aj  yi(x3  +  fix 4)  X\  +  fix 2 J 


\  \(x « 

*3\  4. 

(r  O' 

\  fX2 

aux! 

j  [{^3 

*1)  + 

1° 

j  \ix 4 

X2)\ 

(3.42) 


where  {x}l  |  E  C  C  Z [i]  are  transmitted  symhols,  and  C  is  a  signal  constellation  taken  to  be 
2m-QAM.  The  parameters  t  =  (1  4-  \/5)/2  and  fi  —  {\  —  \/5)/2  are  the  Golden  Ratio.  The 
diagonal  matrix  diag[a,a]  where  a  =  (1  +  ifi)/\/5,  ft  =  (1  4-  ir)/y/ 5  serves  to  equalize 
transmitted  signal  power  across  the  two  transmit  antennas. 

Let  (n  1 .  ri2)  and  (7*21 .  r2 2)  he  the  two  received  signal  vectors  at  the  first  and  second  receive 
antennas  and  the  components  are  the  signals  received  over  two  consecutive  time  slots.  The 
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received  signal  can  be  represented  as 


(Xi  X3\  (X  2  X4\ 

ix 3  XI  j  +  (rhiuvJtoi)  ^-l4  X2J  +  (nn,ni2), 

(Xi  X3\  (  X2  X4 \ 

ix 3  xi )  +  (^12,^22)  (yx4  X2j  +  (n12.«2a), 


(3.43) 


where  ntJ  are  complex  Gaussian  random  variables  with  zero  mean  and  variance  2rr2.  Equa¬ 
tion  (3.43)  can  be  rewritten  as. 


(rn,ria)  =  (x,,*3)(I^,1 
(r21,r22)  =  (*i,x3) 


oJl2l\ 
a/in  J 

d7i22\ 
a/112  J 


.  ,  ,  (  rhn  nfi2\\ 

+  (X2'X4)  U«/*2,  rh„J 

,  /  .  v  /  rhn  lih22\ 

-h  (a  2, 0*4)  r/ll2  J 


+  (mi,  m2), 

+  (7121,1122). 


(3  44) 


Given  that  the  channel  gains  htJ  are  known  at  the  receivers  and  each  symbol  is  transmitted 
with  equal  probability,  optimum  decoding  is  provided  by  the  maximum  a  posterior  MAP/ML 
estimate  as  follows.  We  can  rewrite  (3.44)  as 


r  =  cH  +  sQ  +  n 


(3.45) 


where 


r  =  (rii1r12lrai,r22)»  n  =  (nu,ni2,u2i,n22),  c=  (xi,x3),  s  =  (x2,x4). 

ah2\  r*/u2  a/i22\  ^  /or /in  (\fih2 1  ar/n2  a/i/i22\ 

n  ~  l  ia/i2i  a/in  ia hw  ahn y  ’  ^  —  l  iQ/1/121  ar/in  iot/i/122  ar/ii2  J  * 


(3.46) 


The  likelihood  function  of  codewords  c  and  s  given  the  received  signal  r  is 


p(r  |c,  s)  x  exp 


("2^l|r- 


cH-sQ  ||2  . 


(3  47) 


Taking  the  prior  distribution  of  the  symbols  c  and  5  to  be  uniform  on  the  constellation  C,  we 
obtain  the  maximum  likelihood  estimate: 


(c,  5)  =  argmaxp(r  c,  s). 
C.aec2 


(3.48) 


3.6.2  Quadratic  Decoding  of  the  Golden  Code 

In  this  section  we  show  that  decoding  the  Golden  Code  with  essentially  ML  performance  can 
be  achieved  with  a  simple  algorithm  which  is  quadratic  in  size  of  the  underlying  QAM  signal 
constellation.  Our  approach  is  an  extension  of  ideas  used  to  derive  the  fast  optimal  algorithm 
for  multiplexing  orthogonal  designs  described  in  the  Section  3.3. 
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Optimal  Decoding:  0(N3) 

We  begin  by  rewriting  (3.44)  as 


r  =  x\h  +  xH  +  n  (3.49) 

where  r  isdehned  in  (3.46),  x  =  (.t2,x3,.T4),  h  =  (a/in,  a/i2i,  a/»i2, a/122)  and 

(ar/in  aph2\  atrh\2  Qf.ih2 2\ 

ia/?2i  a/iii  ia/i22  a/ii2  J  .  (3.50) 

iaph2\  arhu  iapd i22  ar/ii2/ 

The  likelihood  function  associated  with  (3.49)  is 

p(r|xi,a:2,X3,a,4)  oc  exp  r  -  a-i/i  -  xh||  j.  (3.51) 

Wc  have  isolated  a*i,  as  our  approach  is  to  maximize  the  likelihood  with  respect  to  xi  given  x , 
and  then  maximize  the  resulting  partially  optimized  likelihood  with  respect  to  x.  To  this  end, 
we  rewrite  (3.51 )  in  the  form 


;j(r |ati,*)  oc  exp  ^-^(r  -  xH)  ^I4  -  jj^p)  (r 

X  CXP  II11  "  •’Jlf;r)l|2j  ’ 

where  \\h\\2  is  a  Euclidean  norm  and 


Thus,  given  x, 

a-i(x)  =  Q(xi(x)), 

maximizes  (3.52). 

Substituting  (3.54)  into  (3.51)  we  obtain  the  optimization  problem  for  x: 


x  =  argmin  r  -  x7i  -  xi(x)h 

x£C3  ' 


(3.52) 


(3.53) 

(3.54) 


(3.55) 


The  ML  solution  is  then  (xi,x2,X3? X4)  =  (xi(x),x).  Thus  (3.54)  and  (3.55)  provide  an 
algorithm  for  obtaining  the  ML  estimate  of  xt,i  =  1, ...  ,4,  which  involves  at  most  N3 
evaluations  of  the  right  hand  side  of  (3.55)  since  there  are  N3  possibilities  to  choose  the  symbol 
set  (x2, X3,.T4)  and  (2(1)  complexity  for  calculating  the  symbol  xi  from  (3.53)  and  (3.54). 
Thus  we  have  an  algorithm  with  ML  performance  with  a  complexity  G(N 3),  where  N  is  the 
size  of  the  underlying  QAM  constellation. 


Optimal  Decoding:  0(N2yfN) 
See  Section  4.4.1 
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Essential ly  Optittial  Quadratic  Decoding:  Q(N2) 


Consider  the  likelihood  function  given  in  (3.47).  We  am  optimize  (3.47)  in  two  steps,  i.e  first 
maximize  with  respect  to  s  given  c  or  vice  versa.  Suppose  we  maximize  with  respect  to  s 
given  c.  We  first  write  (3.47)  as 


p(r|c,  s )  oc  exp 
x  exp 


( 

(- 


2~2  (r  cW)(I,  )(r 

2^2  (*  -  Hc))gg'(s  -  stc))’ j 


(3.56) 


where 

3(c)  =  (r  -  cH)G'(GG')~1  (3.57) 

We  now  make  what  is  essentially  a  zero  forcing  approximation,  since  QQ[  is  not  generally  a 
multiple  of  the  identity.  We  take 


a(c)  =  Q(S(c)) 

=  (Q(Si(c)),Q(S2(c))), 


Substituting  (3.58)  into  (3.47)  we  thus  estimate  c  and  s  as  follows: 

c  =  argmin  \\r  -  cH  -  -s(c)^||2  , 

c€C2 


(3.58) 

(3.59) 


s  =  Q  (s(c)) ,  (3.60) 

where  s(c)  are  given  in  (3.57). 

Alternatively,  if  we  reverse  the  roles  of  c  and  s,  we  obtain  the  estimate 

s  =  argmin  ||r  -  c{s)H  sQ\\2  ,  (3.61) 

«€C2 

c  =  Q(c(s)) ,  (3.62) 

where 

c(s)  =  (r  -  sG)H\HH.'y\  (3.63) 

and 

c(s)  =  Q(c(s)).  (3.64) 

Thus,  we  have  two  possible  decoding  solutions  (3.59)  and  (3.61).  Of  course,  if  HEO  and 
QQ *  were  multiples  of  the  2  x  2  identity  matrix,  both  optimizations  (3.59)  and  (3.61 )  would  be 
exact  ML,  and  we  would  not  need  to  make  a  choice.  However,  as  we  are  making  a  zero  forcing 
approximation,  we  need  to  choose  the  best  alternative  for  each  channel  In  order  to  make  the 
best  choice,  let's  first  consider  the  following  covariance  matrices 


™t=(&*  a) ab)  065) 

where 

a  =  |«|2||/l,||2+|d|2||M|2  (3.66) 

o  =  l«l2  IIM2  +  M2  IIM2 

b  =  -ma*(huh2i  +  *12/122)  +  aa*(/i2i*Ji  +  *22*12) 
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where  ||h,||2  =  |/m| z  -f  |/i *2 12  The  accuracy  of  the  quantization  in  (3.60)  and  (3.62),  de¬ 
pends  on  both  the  determinant  (which  determines  the  SNR)  and  the  condition  number  (which 
determines  the  accuracy  of  the  zero  forcing  approximation)  of  or  QQ^ .  The  following 
compares  the  determinant  and  the  condition  numher  of  TiTi *  and  QQ^ . 


HH* 

QQ] 

Determinant 

a- 2  ~  W 

a2  -  |6|2 

Condition  Number 

y  -  *±IM 

7  “  a-|6| 

~ 

7  o-|*>l 

Fortunately,  it  can  he  seen  that  for  'H'hfl  and  the  matrix  with  the  largest  determinant 
also  has  the  smallest  condition  numher.  Thus,  wc  have  a  clear  choice:  i.e.  if 

(let (gg])  >  dot (3.67) 

then  we  estimate  c  and  s  using  (3.59),  otherwise,  use  (3.61). 

The  condition  (3.67)  is  equivalent  to  a  >  a,  which  in  turn  reduces  to 

IIM2  >  IIM2.  (3.68) 

The  simulation  results  presented  in  Fig.  3.3  show  use  of  our  fast  decoding  algorithm  involves 
little  loss  in  performance  compared  to  the  optimal  ML  decoder.  Wc  gain  some  understanding 
of  why  this  is  by  examining  the  joint  behavior  of  the  condition  numbers  for  and  QQ\  7 
and  7.  We  find  empirically  that  for  i  i.d  Gaussian  channel  coefficients,  although  the  condition 
numbers  7  and  7  can  individually  be  large,  the  minimum  of  the  two  (111111(7, 7))  has  a  very  high 
probability  of  being  small.  Fig.  3.2  shows  the  distribution  of  min(7, 7),  for  10<x  realizations. 
The  largest  value  of  11101(7,  7)  obtained  in  the  10*'  realizations  was  approximately  1 7.5.  This 
corresponds  to  a  ratio  of  the  lengths  of  the  major  and  minor  axes  of  the  noise  ellipse  of  v  17.5  % 
4.18. 

Note  that  the  minimum  condition  number  is  not  hounded  above.  If  equality  holds  in  (3.68) 
then  there  is  a  common  condition  number,  and  when  the  magnitudes  of  the  channel  gains  are 
approximately  equal,  and  their  phases  are  aligned,  then  this  condition  number  can  be  made 
arbitrarily  large  This  however  is  a  very  improbable  event 

To  sum  up,  the  algorithm  involves  one  of  two  possibilities: 

if  IIM2  >  IIMI2 

c  =  argmin  || r  —  cTi  s(c)<J||2, 

r€C2 

*  =  Q  (*(«)). 

otherwise, 

s  —  argmin  ||r  -  c(s)7i  —  sQ ||2  , 
sec2 

c  =  Q  (c(s)), 


where 


s(c)  =  (r-cH)gf(QQ1r\ 

c(s)  =  (r -sQ)V} 
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Figure  3.2:  Empirical  distribution  of  mi  11(7,  7) 


This  algorithm  involves  only  N2  evaluations  of  the  argument  of  (3.59)  or  (3.61 )  compared  to 
the  0(N3)  ML  algorithm.  As  the  simulation  result  shows  in  Figure  3.3,  the  loss  in  performance 
compared  to  ML  is  negligible. 

The  algorithm  works  equally  well  when  there  is  only  one  receive  antenna.  In  this  ease  the 
decoding  problem  becomes 

r  =  ch  +  sg  -F  n,  (3.69) 

where  r  =  (7*11,7*12),  n  =  (m  1,7112)  and  h,g  are  defined  in  (3.44).  The  algorithm  applies 
with  /i,  (j  in  place  of  7 ~C,  Q. 


3.6.3  Simulation  Results 

We  compare  the  performance  of  the  fast  decoding  algorithm  {Quadratic  Decoder)  described 
in  the  last  section  with  the  ML  decoder.  The  channel  is  assumed  to  be  known  at  the  receiver. 
The  elements  of  the  channel  matrix  are  modelled  as  samples  of  independent  complex  Gaussian 
random  variables  with  zero  mean  and  variance  0.5  per  real  dimension.  The  noise  is  complex 
white  Gaussian  with  zero  mean  and  variance  2<r2  The  signal  to  noise  ratio  at  a  receive  antenna 
is  defined  as 

SNR(dB)  =  101og10  (Jffj  ,  (3.70) 

where  Ps  is  the  (average)  signal  power  per  symbol  at  a  receive  antenna  which  is  defined  as 

P,  =  +  Hell2),  (3.71) 

and  Es  is  the  average  energy  per  symbol. 

Figure  3.3,  shows  the  symbol  error  rate  for  4-QAM  and  16-QAM  as  a  function  of  SNR.  The 
simulation  shows  that  in  both  cases  the  performance  of  the  Quadratic  Decoder  is  essentially 
optimal  (ML  decoder). 
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Figure  3.3:  Performance  comparison  between  the  MAP/ML  0(N3)  decoder  and  the  quadratic 
decoder  for  4 -QAM  and  16-QAM 

3.7  Fast  Decoding  of  3  x  3  Perfect  STBC 


Oggier  ct  al  [7]  introduced  the  perfect  space-time  block  codes  w  hich  satisfy  all  of  the  following 
criteria-  full-rate,  full-diversity,  non-vanishing  determinant,  good  shaping  and  uniform  average 
transmitted  energy  per  antenna.  These  codes  are  constructed  for  2  x  2, 3  x  3, 4  x  4  and  6x6 
Ml  MO  systems.  An  example  of  a  2  x  2  perfect  STBC  is  the  Golden  Code  [22]  described  in 
Section  3.6.  The  conventional  ML  decoder  for  perfect  STBCs  with  an  Ar-QAM  or  Ar-HEX 
constellation  is  based  on  an  implementation  of  sphere  decoding.  Although  there  is  no  report 
on  the  decoding  complexity  for  other  perfect  codes,  it  is  expected  that  when  the  channel  matrix 
is  close  to  singular,  the  preprocessing  stage  of  the  sphere  decoding  algorithm  will  yield  a  plane 
of  possibilities  rather  than  a  single  initial  estimate.  When  this  occurs,  the  lattice  point  search 
degenerates  to  an  exhaustive  search. 

In  this  Section  we  show  that  the  approach  described  in  Section  3.6  ean  be  applied  to  ohtain  a 
fast  decoding  algorithm  for  3  x  3  perfeet  STBCs  which  gives  essentially  ML  performance  with 
substantially  reduced  complexity.  This  approach  ean  be  also  applied  to  other  perfeet  STBCs  to 
obtain  low  complexity  decoding  with  essentially  ML  performance. 


3.7.1  Model  and  Decoding  Problem 

Consider  the  codeword  matrices  of  the  perfeet  STBCs  for  3  x  3  MI  MO  systems  w  hich  we  write 
in  a  form  that  will  assist  in  the  development  of  our  algorithm.  The  perfeet  3x3  STBC  transmits 
nine  complex  (usually  N- HEX  constellation)  information  symbols  {.Ti}?  ,  over  three  time 
slots  from  three  transmit  antennas.  The  transmit  codewords  of  the  3  x  3  perfect  STBC  can  be 
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expressed  as 


2 

X  = 

I  0 


/  *3i+l 
JJ3l  |  3 

yj f  2 


#34  +  2 


#3i  +  I 
J#3i  \  3 


#3i  +3^ 
#3*  4  2 

#3i  4  1  J 


where  the  diagonal  matrices  Bx  are 

B\  =(i  +  j)/3  +  e 


B2  =  (-l-2j)73  +  j02 

b3  =  (-  i  2j)/3  +  (i  +j)o  +  (i  +j)e2 


with  0  =  diag(0i,02,03),0i  =  2cos(2*7t/7),  j  =  c2,r,/3  and  B,  satisfy 

Tr(BJB,tn)  =  76]m. 


(3.72) 


(3.73) 


(3.74) 


Assume  that  the  channel  state  information  is  available  at  the  receiver.  Let  hXJ  be  the  channel 
gain  from  transmit  antenna  i  to  a  receive  antenna  j  then  the  received  signal  is  given  by 


Y  =  HX  +  W 


where 


(/in  /i2i  h$\\ 

h\  2  h 22  h  32 

h\3  /i23  7*33  j 


Equation  (3.75)  can  be  rewritten  as 

r  =  (x\yX2>x$)'H\  +  (x4,xs,xe>)?i2  +  (#7,  #8,  #9)7*3  +  n 


(3.75) 


(3.76) 


(3.77) 


where  r  =  (n,  r2  r3)  contains  the  three  received  signal  vectors,  r\  —  (rx\ ,  r*2i  rt3)  with  the 
component  rtJ  representing  the  received  signal  at  antenna  i  in  lime  slot  j.  The  noise  n  is  i .i  d 
Gausian  noise  with  zero  mean  and  covariance  2 <t2/3  and 


n  i 


(IIuGuC^ni  =  (//2,G2,r2).7*3 


(/73,  C73lC  3), 


(3.78) 


where  77 M  G,  and  are  induced  channel  matrices  from  the  three  transmit  antennas  to  the  first, 

second  and  the  third  receive  antenna  respectively.  Explicitly 


IIi  = 


/  biihu  bi2h2\ 
jbisli3 1  i  /?  1 1 
\^jbi2^2l  J  bi3 /l3  I 


bt3fl3l\ 
^2^21 
6,  l  /l  1 1  y 


(3.79) 


and  similarly  for  Gt  and  Ct.  The  induced  channel  matrices  7*i,7*2,7*3  have  the  following 
property  which  is  the  basis  of  our  fast  decoding  algorithm 

-H.-Hl  +  HtH't  +  HsHl  =  7  ||H||2, 1.j.  (3  80) 

Thai  is  Z*  ,  7 i,nj  is  a  multiple  of  the  identity  A  similar  property  of  the  induced  channel 
matrices  holds  for  all  of  the  perfect  STBCs,  including  the  Golden  Code.  In  fact,  our  fast 
decoding  method  will  apply  to  any  STBC  with  structure  giving  rise  to  a  relation  of  the  form 
(3.80). 
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3.7.2  Low  Complexity  Decoding  Algorithm 

Let  c  =  (xi,  X2,x3),  s  =  (X4,  X5,  x&),  y  =  (X7,  xg,  X9),  we  can  rewrite  (3.77)  as 


r  =  cH  1  +  s'ht 2  4  +  n  (3.81) 

The  likelihood  function  associated  with  (3.81)  is 


p{r\c,  s,  y, )  oc  oxp  I |r  -  cHi  -  sTii  -  t/7-^3||2^  .  (3.82) 

Based  on  the  conditional  optimization  described  in  Section  3.4,  we  first  maximize  (3.82)  with 
respect  to  s  and  y  given  c. 


p(r|c,  s,  y)  oc  exp 
x  exp 


("2 j^r'(l9-*l(W,«lr,W,)r,t) 

^-5^2  (c  c(s,y))'H\'H\  (c  -  c(s,y))f  j 


(3.83) 


where  r'  =  r  —  sTiz  —  y7i3,  and 

c(s,  y)  =  (r-sHi  -  .  (3.84) 

We  now  make  what  is  essentially  a  zero  forcing  approximation,  since  T-l\T-l\  is  not  generally 
a  multiple  of  the  identity.  We  take 

c(«,  y)  =■  Q(c(s,  y))  =  (Q(xt(s  y)),Q(f2(s,y)),Q(x3(s.y))).  (3.85) 

Substituting  (3.85)  into  (3.82)  we  thus  estimate  c,  s  and  y  as  follows: 

(s,  y)  =  argmin  ||r  -  c(s,  y)Wi  sW2  -  yW3||2 

’■y^c  (3.86) 

c  =  Q(c(s,y)) 


where  c(s  y)  is  given  in  (3.84). 

If  we  first  maximize  (3.82)  with  respect  to  c  and  y  given  s  we  obtain  the  estimate 

(c,  y)  =  arg  mill  ||r  -  cTi  1  -  s(c,  y)  -  yfi-t ||2  , 

r«ec  (3.87) 

s  =  Q  (S(c,y)) 


where 


3(c,y)  =  (r  rTi,  -  y'H3)'H\{'H.2'H\y\ 
Alternatively,  if  we  maximize  (3.82)  with  respect  to  c  and  s  given  y  we  obtain 

(c,  s)  =  arg rnin  ||r  -  cH\  -  s'H2  -  y(c.  s)7i^\\2  , 
c.a€C 

y  =  Q(y(c,s)) 


(3.88) 


(3.89) 


where 

y(c  a)  =  (r  -  cTf,  -  (3.90) 
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Figure  3.4:  Empirical  distribution  of  minimi ,  72, 73) 


Equations  (3.86),  (3.87)  and  (3.89)  each  provide  an  algorithm  for  obtaining  the  estimate  of 

xt1i  —  1 . 9,  each  of  which  involves  at  most  /V6  evaluations  of  the  right  hand  side  of 

one  of  (3.86),  (3  87)  and  (3.89).  Now,  we  have  three  possible  decoding  solutions.  Of  course 
and  Ht.'H.I  were  multiples  of  the  identity  matrix,  all  of  the  optimizations 
(3.86),  (3.87)  and  (3.89)  would  have  exact  ML  solutions  and  we  would  not  need  to  make  a 
choice.  However,  as  we  are  making  a  zero  forcing  approximation,  we  need  to  choose  the  best 
alternative  for  each  channel.  One  approach  is  to  compute  all  three  alternatives  and  take  the 
alternative  which  maximizes  the  likelihood  The  key  to  the  current  algorithm  is  that  due  to  the 
structure  of  the  code  one  of  the  three  estimates  is  good,  i.e.  essentially  ML,  with  very  high 
probability  The  accuracy  of  the  quantization  depends  on  both  the  determinant  (which  deter¬ 
mines  signal  to  noise  ratio)  and  condition  number  (which  determines  the  accuracy  of  the  zero 
forcing  approximation)  of  Ft \Fi\  FtvFt\  or  >£3^3.  Fig. 3.4  shows  the  empirical  distribution 
of  min(7i ,  72, 73)  for  i  i.d  Gaussian  channel  coefficients,  where  71,72,73  represent  the  con¬ 
dition  numbers.  This  shows  that  although  the  condition  numbers  can  individually  be  large,  the 
minimum  of  the  three  has  a  very  high  probability  of  being  small 
We  can  reduce  the  computation  by  a  factor  of  three  by  deciding  on  one  of  the  three  estimates 
based  on  the  channel.  A  possible  criterion  is  to  choose  to  quantize  the  variables  corresponding 
to  the  7ij  with  the  largest  value  of  det (FijFt^).  Another  choice  is  to  quantize  the  variable 
for  which  the  corresponding  matrix  Ft}Ft]  has  the  smallest  condition  number.  For  the  Golden 
Code  these  two  criteria  are  equivalent,  but  here  they  are  not.  We  have  found  experimentally  that 
the  former  criterion  is  just  slightly  better  and  obviates  the  need  to  compute  all  three  estimates 
The  algorithm  can  be  summarized  as  follows: 


45 


Figure  3.5  Performance  of  the  0(N6)  decoder  and  the  ML  0(N8)  decoder 


Let  det('H)  =  max[det(7i,7i|),det(7i2'Ht2),det('H3'W^)]. 
Ifdet(7i,7it1)  =  det(7t) 

(s,y)  =  argmin  ||r  -  c(s,y)?ii  -  sH?  -  y'Ha ||2 
i  wee 

c  =  Q (c(s.  y)) 
clscifdet(W27i^)  =  det(7i) 

(c,y)  =  arginin  ||r  -  c7i i  s(c.y)W2  y7^3 1|2 

C.W€C 

s  =  Q  (5(c,y)) 


otherwise 

(c,  s)  =  argmin ||r  -  cT-C}  s'H2  -  y(c,  s )7-^ ^ j| 2 

c,s€C 

y  =  Q(y(c,  s)) 

where  c,  s,  y  are  given  in  (3.84).  (3.88),  (3.90)  respectively. 

The  perfect  STBCs  are  constructed  in  terms  of  information  symbols,  either  a  QAM  or  HKX 
constellation.  This  means  that  the  computational  complexity  of  the  quantization  step  is  0(1). 
Therefore  our  algorithm  involves  at  most  N6  evaluations  of  the  likelihood  function.  We  com¬ 
pare  the  performance  of  the  fast  decoding  algorithm  described  above  with  the  ML  decoder  for 
a  3  x  3  perfect  STBC.  Figure  3.5  shows  the  symbol  error  rate  as  a  function  of  SNR  using 
a  4-HEX  constellation  for  a  Rayleigh  fading  channel  model.  The  result  shows  that  our  fast 
decoder  is  essentially  a  ML  decoder  with  complexity  0(N6). 
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Chapter  4:  Fast  Decodable  Space-Time  Block  Codes: 
General  Theory 


In  this  chapter  we  consider  how  to  write  linear  dispersion  codes  and  in  particular,  the  perfect 
STBCs,  in  a  way  that  makes  the  decoding  structure  clear.  It  will  become  clear  that  the  per¬ 
fect  STBCs  can  be  considered  as  multiplexed  quasi-orthogonal  designs.  Application  of  the 
conditional  optimization  approach  developed  in  Chapter  3  lead  to  fast  exact  ML  decoders. 

Consider  a  MIMO  communication  system  with  nt  transmit  antennas  and  nr  receive  anten¬ 
nas.  The  transmitted  signal  matrix  X  €  Cn'  xT  and  the  received  signal  matrix  Y  €  Cn,xT 
are  related  by 

Y  =  HX  +  W  (4.1) 

where  H  €  CrirXn*  denotes  the  ehannel  matnx  with  entry  fitj  representing  the  fading  coef¬ 
ficient  associated  with  the  transmit  antenna  i  and  receive  antenna  j  The  channel  coefficients 
are  i  i.d.  circularly  symmetric  complex  Gaussian  random  variables  with  zero  mean  and  unit 
variance.  W  €  Cn,xT  denotes  the  complex  additive  white  Gaussian  noise  with  i.i.d.  en¬ 
tries  W tJ  ~  CA/"(0,  2<t2).  We  assume  quasi-static  flat  Rayleigh  fading  where  the  channel 
coefficients  are  fixed  for  the  codeword  duration  i.e.  for  T  symbol  periods. 


4.1  Linear  Dispersion  Codes 

Let  be  a  set  of  real  scalar  symbols  which  are  selected  in  pairs  from  some  constella¬ 

tion  A  linear  dispersion  (LD)  code  [34]  is  a  space-time  code  in  which  codeword  matrices  are 
obtained  as  linear  combinations  of  certain  basis  matrices  according  to 

2  K 

X  =  ]T:r*.4*,  (4.2) 

k=\ 

where  €  Cn>  xT  is  a  set  of  complex  matrices  If  the  Xk  roam  over  all  of  R  this  is  the 

real  linear  span  of  the  matrices  Ak  and  as  such  is  a  real  vector  subspace  of  Crif  xT.  That  is,  a 
subspace  of  Cn'  x  i  which  is  a  vector  space  over  R.  We  denote  this  space  by  V. 

Note  that  LD  codes  encompass  all  possible  linear  space-time  codes.  Substituting  (4.2)  into 
(4. 1 ),  the  received  signal  vector  y  can  be  w  ritten  as 


where  x  € 
matnx 


y  ~  xTi  4  n  (4.3) 

is  the  transmitted  information  symbol  vector,  and  7i  is  the  induced  channel 

(  h\A\  /i2^i  -  •  hnt  A\  \ 


n  - 


(4.4) 


\h\A2K  h2A'2K  •••  hn,  A'lK  j 


where /i,  =  (/iw,  h2i, . . . ,  /?„,*). 
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4.2  Algebraic  Properties 


In  this  section  wc  look  at  some  algebraic  properties  of  an  LI)  codes  which  leads  to  a  fast  ML 
decodable  property.  We  begin  by  recalling  the  definition  of  an  orthogonal  design. 

Definition  2  If  an  LD  code  (4.2)  is  an  orthogonal  design  then  for  all  k,  l  —  1, . . . ,  2  K, 

AkA\  +  A,A[  =  ak6uI,  ak  >  0.  (4.5) 

There  are  many  codes  which  arc  not  themselves  an  orthogonal  design  but  do  have  subcodes 
which  are  orthogonal  designs.  Following  [9]  we  refer  to  these  codes  as  having  embedded 
orthogonality.  More  precisely: 

Definition  3  An  LD  code  (4  2)  is  said  to  have  embedded  orthogonality  of  order  m  if  some 
subset  of  size  in  of  the  mat  rices  Ah  satisfies  (4.5).  For  such  a  code  we  can  relabel  the  Ah,  and 
write  (4.2)  as 

2m  2  K 

X  =  ^XhAh  A  ^2  XkAk  (4.6) 

k  1  It  2m  +  l 

where  satisfies  the  condition  (4.5). 

Definition  4  An  LD  code  (4.2)  is  a  (2  K,  d)-Quasi  Orthogonal  Design  ((2  K  d)-QOD)  if  there 
is  a  partition  of  the  set  of  integers  { 1  ....  2  K}  into  2  K/d  sets  S\ , . . . ,  S2k  fd  cflc//  of  size  d, 
such  that 

L  For  each  j  €  $e  and  m  €  S^,  with  £  / 

AjAl  +  AmA)  =  0,  (4.7) 

2.  For  all  f  =  1,...,2  K/d 

J2  AjA\  =  <*,Ia/  (4.8) 

j€Sr 

Note  that  a  {2K,  l)-QOD  is  an  orthogonal  design 

Definition  5  An  LD  code  is  said  to  have  embedded  (2?7i,  rf)-quasi-orthogonality  if  some  subset 
of  the  matrices  Ah  of  size  2m  satisfies  the  conditions  (4.7)  and  (4.8).  For  such  a  code  we  can 
relabel  the  Ah  (and  x kh  and  write  (4.2)  as 

2  m/d  id  2  K 

x=£  £  XhAh  A  ^2  XkAk>  (4.9) 

Ar-(£-l)d-H  k  2m +  1 

where  {/U}jt^j  satisfies  the  (2 m,d)-QOD  conditions  with  partition  S(  =  {(£  —  \)d  A 
1, . . . ,  Id},  for  i  =  1, . . . ,  2 m/d. 
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4.3  Fast  Maximum  Likelihood  Decoding 

Consider  the  received  signal  given  in  (4.3).  The  likelihood  function  of  symbols  x  given  the 
received  signal  y  is  given  by 

p(y\x)  oc  cxp|“2^2  II y  ~  ^H2}  •  (4.10) 

The  maximum  likelihood  estimate  of  x  is: 

x  =  arg  maxp(y|x).  (4.11) 

X 

If  an  LD  code  X  in  (4.2)  is  an  orthogonal  design  then  the  ML  decoding  problem  forrz>,  k  ~ 
1  ...  ,2 K  decouples  into  2K  real  parallel  sub-channels. 

Theorem  4  If  a  linear  dispersion  code  of  the  fonn  (4.3)  has  embedded  orthogonality t  of  order 
m  then  exact  ML  decoding  can  be  implemented  with  complexity'  0(NK~m)  for  square  N- 
QAM  and  with  complexity  0(NK~rn^ *)  for  a  general  constellation 

Remark  5  It  follows  that  in  all  cases  the  complexity'  of  ML  decoding  is  at  most  0(NK~l)  for 
square  N -QAM. 

Proof  Suppose  our  code  X  has  the  form  (4.6).  Write  x  \  -  (ari,  ■  •  ■  ,  £2m)  and  x2  = 
(x2m  +  i,  •  •  •  ,X2k).  and  decompose  (4.3)  as 

y  =  xiTii  +  x2'H2  +  n  (4  12) 

The  likelihood  function  associated  with  (4.12)  is 

p(V |*J, *2)  OC  exp  II y  -  a'l'Hi  -  a:27-£2||2  j- .  (4.13) 

We  have  isolated  x\,  as  our  approach  is  to  maximize  the  likelihood  function  (4.13)  with  respect 
to  x\  given  x2,  and  then  maximize  the  resulting  partially  optimized  likelihood  function  with 
respect  to  x2. 

Since  the  matrices  {zU,  k  —  1, . . . ,  2??i}  satisfy  the  condition  (4.5), 

nin\  +nin\  =  d,  (4.14) 

where  D  —  diag(ai,  •  •  •  ,a2m),  the  likelihood  function  (4.13)  can  be  written  as 

p(y\xux2)  oc  exp  (  -  i^j(y  ~  X2Ti2)(y  -  x2 H2)' 

i  2m  .  2m  (4.15) 

+  4^2  Hrtjli:j(a:2)|2}  n«p{-ftl-i  -ij(*2)|2}, 

j  =  1 

where  _ _ 

i.(x2)  =  (yn\  +  yn\  x2(n2n\  4  n2n\)}  n~\  (4.16) 
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Thus,  given  X2<  the  ML  dceoding  problem  for  x\  decouples  into  m  parallel  sub-decoding 
problems  We  have 


Xi(x2)  =(Q(x1(x2)),- ••  ,Q(im(x2))) 

(4. 17) 

sQ(*i(x2)), 

where  Q  denotes  the  ML  decoder  for  the  sub-decoding  problems.  For  square  Ar-QAM  the  ML 
decoder  Q  is  just  quantization  which  has  complexity  (7(1 ).  For  a  general  constellation  Q  will 
involve  a  search  with  linear  complexity  in  the  size  of  the  constellation.  Substituting  (4  17)  into 
(4. 15)  we  obtain  the  optimization  problem  for  X2'. 

X2  —  arginiii  \\y  -  xx  (x2)'H\  -  xiH 2||2  .  (4.18) 

X2 

which  can  be  substituted  back  into  (4  17)  to  obtain  the  ML  estimate  of  X\ 

Xi-Xi(x2).  (4.19) 

Thus,  (4  18)  and  (4. 19)  provide  an  algorithm  for  obtaining  the  ML  estimate  of  xt,i  —  1, . . .  ,  2/\, 
which  involves  at  most  Nh~7n  evaluations  of  the  right  hand  size  of  (4.18)  for  square  N- 
QAM.  Lor  a  general  constellation  we  need  to  evaluate  (4.18)  Nh  times  and  each  evalu¬ 
ation  involves  2m  parallel  searches  through  N  alternatives  to  give  an  overall  complexity  of 
0(NK~m+l).  m 


Theorem  5  If  an  LD  code  (4.2)  has  embedded  (2m,  d) -quasi -orthogonality,  then  exact  ML 
decoding  can  be  implemented  with  complexity  0(Nh  jor  transmj[[e(]  symbols 

taken  from  a  square  N-QAM  constellation  while  for  a  general  constellation  the  complexity'  is 

0(NK~m*d). 


Proof  Writexf  =  (x«-i>d+i, 
Then  we  ean  rewrite  (4.3)  as 


Tx^i),  i  =  1,  •  •  •  ,2m  andx2m  +  i  =  (x2m+i,  •  •  •  , x2k)- 


2  tn/d 

y  =  XtHt  +  X2m  +  i'Hm  1 1  +  n. 


(4.20) 


The  likelihood  function  associated  with  (4  20)  is 


p(y\{xf}:x2m  +  i)  °c  OXp  - 


2a2 


2  m/d 


y  ^  "  xfHt  a?2m+i7^ 


2m  +  1 


(4.21) 


Now  since  the  matrices  {A^,  k  =  1 ,  •  •  •  ,  2m}  satisfy  the  condition  (4.7)  then,  for  t  ^ 


ntn\,  +  utu\,  =  [  M  1. 


Mj  0, 


(4.22) 
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where  is  a  d  x  d  matrix  and  0^  denotes  ad  x  d  matrix  with  all  elements  are  zero,  and  so 
the  likelihood  function  (4.21)  can  be  written  as 

p(y|{Xf},X2mu) 

a  exp  |  X2m  1  ~  X2m* 1^2m4i)t 

j  2  m/d  _ 

+  ^  Xl(X2m+  l)  (li.e'Hl  +1tt g<(X2m4-l)T| 

t-  1 

2m  fd  j  _ 

X  Y[  exp{(“  ~~  *#(*2m+l))(W«?i{  +  'Ht*H\)(Xt  -  *f(*2rn+l))T}, 

t  1 

(4.23) 


where,  for  (.  =  l,  •  ,  2 m/d, 

X  f  (  X  2  m  +  1  ) 

=  (yttj  +  yWl  -  seamti(Wam+iWj  +  Wam*i?4))  (w.wj  +  W,W*)  . 

(4.24) 


Thus,  given  X2m  +  ii  the  ML  decoding  problem  for  {x^f1  =  1,  •  ,  2 m/d}  decouples  into 

2 m/d  independent  sub-decoding  problems  We  have 

Xp(x2rn  +  \)  —  mill  (x^  X  f  (x  2m  f  1  ))  (  f  4  |  j  (x^-x^(x2m  +  l  ))T.  (4.25) 

Xf  '  ' 


For  a  square  N- QAM  constellation  the  complexity  of  the  above  2  m/d  minimization  is  0(N^d~ 1  )/2) 
using  conditional  optimization,  while  for  a  general  constellation  the  complexity  is  0(Nd). 
Suhstituting  (4. 1 7)  into  (4. 1 5)  we  obtain  the  optimization  problem  for  X2m  *  \ : 


X2m+i  =  argmiri 
*2111+1 


2m/ d 

V  ^  ^  Xf(X2m  +  l)7"^  —  X2frt  l  l7^2;ri  +  l 
*  1 


(426) 


which  can  he  substituted  back  into  (4.25)  to  obtain  the  ML  estimate  of  x< 


X(  =  Xf  (x 2m  + 1 ) ’  t  =  1,  ••  •  ,2 m/d.  (4.27) 

Thus,  (4.25)  and  (4.26)  provide  an  algorithm  for  ohtaimng  the  ML  estimate  of  xt .  i  —  L  ,  2 K 
which  involves  at  most  NK~m  evaluations  of  the  right  hand  size  of  (4.26)  for  square  N- 
QAM,  with  each  evaluation  having  complexity  0(Ar(<i-1)/2),  giving  an  overall  complexity 
of  0( /vK”m  +  (rf_l)/2).  For  a  general  constellation  we  need  to  evaluate  (4.26)  Nh  ~rn  times 
with  each  evaluation  involving  2m /d  parallel  searches  through  Nd  alternatives  giving  an  over¬ 
all  complexity  of  0(NK~m*d).  ■ 

Of  course,  we  can  apply  this  general  theorem  to  reproduce  the  conclusions  of  Chapter  3  for 
multiplexed  orthogonal  [55]  and  quasi -orthogonal  designs  [56].  However,  somewhat  surpris¬ 
ingly,  this  theorem  also  applies  to  the  perfect  space-time  hloek  codes  of  Oggier  et  al  (7]  which 
wc  will  show  below  have  the  structure  of  multiplexed  quasi-orthogonal  designs. 
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4.4  Examples 

4.4.1  The  Golden  Code 

Recall  that  the  Golden  Code  codewords  take  the  form 


X 


1 

7! 


/.Ti  x3\  fr  0\ 

\ix3  xi  J  l  0  nj 


/  X 2  X4 
W.X‘4  X2 


(4  28) 


where  a  =  ( 1  +  ifi)  and  a  =  (1  +  ir).  This  ean  he  written  in  terms  of  a  real  LD  code  as 

8 

X  =  £  skAk  (4.29) 

k  1 


with  ,S*1  =  X\R,S-2  —  X\I,Ss  =  X2R  S4  =  X2/,S5  =  X.3fl,A’6  =  X3I ,  S7  =  X4  R-Sg  =  X4  /, 

where  subsenpts  nR  and  nl  represent  real  and  imaginary  parts  of  symbol  xn,  and 


Ax  = 


A&  = 


a  0 
0  a 

0  a 
ict  0 


»  A2 

,  ^6 


?a  0 


ar  0 


n  1  1  ^3  (  n  I  ’  ^4 

0  ?rv  /  \  0  a// 


0  ia 


n  /  1A7  =  .  .  n 

—  Or  U  /  \  lOtfi  U 


0  ar 


>  As  = 


ictr  0 
0  iafi 

0  iar 
—  a//  0 


(4.30) 


We  ean  interpret  the  Golden  Code  as  a  multiplexed  pair  of  (4,2)-QOD  leading  to  a  fast  ML 
decoder.  To  see  this  consider  the  matrices 

t =(.,“)  ■nj  2=(“l)  ,4-3i> 

and  note  that  ZH  —  I2.  From  T  we  ean  construct  a  pair  of  diagonal  matrices 


B 1  —  (1  +  i) I2  ~  iT 
B2  =  -il2  4-  T 


(4.32) 


which  satisfy 

Tr(BjBl)  =  5  Sjm.  (4.33) 

Note  that  A\  —  B\^  A2  =  B\  22  =  tB\ ,  A3  —  B2 ,  A4  =  B2Z?  —  iB2 ,  ^4a,-4  4  —  ^4* 2,  A:  — 
1,...,4.  It  ean  be  verified  that  the  subsets  {A*}*  j  and  {>1^  =5  satisfy  the  (4,  2)-QOD  con¬ 

ditions  (4.7)  and  (4.8)  with  partition  {1,3},  {2,4}  and  {5,7},  {6,8}  respectively  Therefore 
the  Golden  Code  is  the  multiplexing  of  a  pair  of  (4,2)-QOD  with  multiplexing  matrix  2, 


X  =  Xi  (.Tifl,  X2fl,  Xi/,  X2/)  -F  X2(X3fl,  X4fl,  X3;,  X.i/)2  (4.34) 


where 

4 

Xl(fli,  a2,  a3,  a4)  =  (4.35) 

k  1 


52 


The  received  signal  can  be  written  as 


r  -  x\H\  +  X1H2  +  n, 


where  X\  =  {x\r<  x2/?,xi/,x2/),  x2  =  (x3*,  x4/?,  x3/,  x4/)  and  also 


WiWj  +WiWl 


/a  6  0  0\ 
tf  00 
0  0  a  6 
^0  0  6 


=  w2?4  + 


(4.36) 


(4.37) 


where  a,  6,  c  are  some  constants  Theorem  5  implies  that  the  Golden  Code  has  ML  decoding 
complexity  0(N2y/N)  (see  [57,  52|)  for  square  N  QAM  constellations  using  the  decoding 
algorithm  of  the  last  Section. 


4,4.2  3  x  3  Perfect  STBC 


The  perfect  3  x  3  STBC  transmits  nine  complex  (usually  N-HEX  constellation)  information 
symbols  over  three  time  slots  from  three  antennas.  The  codeword  of  the  3x3  perfect  STBC 
can  be  expressed  as 


2 

X  =  Vb,  +  i 

i  —  0 


/  X3i  fl 
JX3l  + 3 
^?X3i  +  2 


X3l  +  2 

X3t+1 
7X31  +  3 


a*3t  +3\ 

x3i+ 2 
X31  +  1  j 


(4.38) 


where  the  diagonal  matrices  Dt  are 


b  1  =  (i  +  j)/3  +  e 
b2  =  (  i  2 j)h  +  je2 

B3  =  (-1  2  j)l3  +  (1  +  j)S  +  (1  +  j)G2  (4.39) 


with  0  =  diag(0j ,  02, 03 ),  0,  —  2  cos(2'7r/7)  j  =  e2”1  1  and  Bt  satisfy 

Tr  (BjBl)  =  76Jtn.  (4.40) 


Similarly,  let 


Xi  =  Xi/,X4/,X7/), 

X2  =  (X2fi,X5/?,T8/?.X2/tX5/,  Xs/), 

X3  =  (X3K,X6/?,  X9R,  X3/,  X6/,  X9/), 


then  the  3  x  3  perfect  STBC  (4.38)  can  be  written  as  three  (6,3)-QOD  multiplexed  as  follows 


where 


X  =  Xi(x,)  +  X2(i  2)Z  +  X3(i3)22 

6 

X?(.Si,(S2,53,S4,S5,.S6)  = 

k  1 


(4.41) 


(4.42) 
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(4  43) 


Z  = 


/O  1  Q\ 
0  0  1 

vj  0  °y 


and 


A\  =  B\ ,  A2  4-  F$2 1  A3  —  B3 ,  A 4  =  ,  A$  =  i/^2,  =  iBz- 

In  this  case  the  received  signal  is  given  by 

r  =  X\H\  4-  X2H2  4-  X3H3  4-  n 


(4  44) 

(4.45) 


and 

HiH\  +  HiH\  =  n2nl  +  HjHl  =  n^nl  +  n3nl 

/M  j  03  \  (4.46) 

\  O3  M3  J 

Theorem  5  implies  that  the  3  x  3  perfect  STBC  has  ML  decoding  complexity  0(N7)  for  square 
QAM  constellations. 


4.4.3  4x4  Perfect  STBC 


The  4x4  perfect  STBC  transmits  16  complex  (N  QAM  constellation)  information  symbols 
{.t*}]?!  over  four  times  slots  from  four  antennas.  The  codewords  can  be  expressed  as 


x  =  y>t4 


(  1  X4i+2  *T4i  +  3  #4t+4^ 

ix  4if4  ^-4i3-l  ^41+2  ^-4*  +  3 

iX4i  +  3  iX  4ij-4  3*41-1-1  X-4i  +  2 

\iX4i±2  IX  4t+3  ?X4i+4  ^4t  +  l  J 


(4.47) 


where 

Hi  =  (1  -3 i)/4  4-i©2 
B2  =  (1  -3i)B  +  iB3 
lh  =  -H 4  +  (-3  +  4i)0  +  (1  -  t)©3 
B4  =  (-1  +  i)I 4  -  3©  +  Q2  +  ©3 


with  0  =  diag(0i ,  02 1 03. 04),  0i  =  2  eos(2t7r/15),  and 

Tr(BjBl)  =  1 5<5jm 


Write 


X\  =  (.Ti/f,  XtjRy  X9fi,  X\3R,  Xi/,X5/^97,  X\3i)y 
X2  =  (X2R>X6RiX\0R,Xi4R,X2l)XGhX\0I,X\4i)) 
X3  =  (X3R,  X7R,  XUR,  X15R,  X3/,  XJJ,  In/,  X15/). 
X4  —  (.2*4/?,  X%RsX\2R,  3*16/?.  X4I ,  XhJ  ,  X \21 ,  X\Gl) . 


(4.48) 
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The  4x4  perfect  STBC  (4.47)  can  be  written  as  four  (8,  l)-QOD  multiplexed  as  follows 

X=  X,(a:,)  +  X2(a:2)Z  +  X3(*3)Z2  +  X4(x«)Z3,  (4.49) 


where 


and 


8 

Xl(si,52,*‘3»S4,55t«6>57tS8)  =  ^SkAk, 

k-  1 

/o  1  o  o\ 

0  0  10 
0  0  0  1 
0  0  0 J 


(4.50) 


(4.51) 


A i  =  B\ ,  j4'2  —  B'2i  A ^  =  63,  /I4  =  /?4,  A$  =  iB\9A$  =  —  iB,i,A%  =  t/^4. 


The  received  signal  is 


r  =  X1K1  4-  X2H2  4  X3W3  4-  X\H  \  4  n 


(4.52) 


and 


HiH\  +  H\H\  =  n2n\  +  n2n\  =  h3h\  +  n3n\  =  +  h<h\ 


/m4  o,  \ 

(o4  M.J’ 


(4.53) 


Theorem  5  implies  that  the  4x4  perfect  STBC  has  ML  decoding  complexity  ()(A,13\//V)  for 
square  QAM. 


Chapter  5:  Space-Polarization-Time  System 


The  performance  of  MIMO  systems  is  highly  dependent  on  the  availability  of  a  rieh  scattering 
environment  (e.g.  densely  urban  and  indoor)  as  well  as  sufficient  antenna  spacing  to  achieve 
multiplexing  or  diversity  gain.  In  rural  and  remote  environments  the  wide  open  space  and 
flat  or  smooth  undulating  terrain  give  rise  to  LOS  conditions  dominating  between  the  trans¬ 
mitter  and  the  receiver  and  consequently  to  a  loss  of  diversity.  In  order  to  achieve  significant 
multiplexing  or  diversity  gain,  many  wavelengths  spacing  between  antenna  elements  at  the 
base  station,  and  up  to  a  wavelength  at  a  mohile  unit  are  required  [58].  As  the  polarization 
of  transmitted  signals  is  mostly  preserved  by  the  LOS  environment,  the  use  of  dual-polarized 
antennas  (polarization  diversity)  is  a  promising  cost  and  space  effective  alternative,  where  two 
spatially  separated  uni-polarized  antennas  are  replaced  by  a  single  antenna  structure  employing 
orthogonal  polarization.  In  this  chapter  we  investigate  the  performance  of  various  transmission 
schemes  for  coding  onto  a  dual  polarimetric  antenna. 

We  begin  with  a  brief  overview  of  electromagnetic  polarization  and  polarimetic  transmission 
systems  by  following  [59]. 

5.1  Polarization 

Polarization  is  a  property  of  the  electromagnetic  plane  wave  solutions  of  Maxwell’s  equations. 
In  an  infinite  medium  in  which  there  are  no  sources.  Maxwell's  equations  have  so-called  plane 
wave  solutions  of  the  form 

E(x,t)  =  £e'ikx-U1,\  (5.1) 

B(x,t)  =  Be,(fc'x-U"\  (5.2) 

where  t  denotes  time,  x  €  R5  denotes  a  point  in  space  and  £  and  B  g  C3  The  (angular) 
frequency  of  the  plain  wave  is  w  and  k  is  called  the  wave  vector  of  the  plane  wave.  In  order  to 
be  solutions  of  Maxwell's  equations  the  following  relations  must  hold 


k 

uj 

k  =  /tc  — , 

Cr 

(5.3) 

k 

£  =  0, 

(5.4) 

k 

5-0, 

(5.5) 

B 

=  y/Jl€  k  X  £ 

(5.6) 

Equations  (5. 1 )  and  (5  2)  are  a  complex  solution  to  Maxwell's  equations  Physical  solutions  are 
given  by  both  the  real  and  imaginary  components  of  this  solution  The  unit  vector  k  determines 
the  direction  of  propagation  of  the  plane  wave.  Equations  (5.4)  and  (5.5)  imply  that  the  E  and 
B  are  in  the  plane  perpendicular  to  the  direction  of  propagation  k.  Furthermore,  equation  (5.6) 
implies  that  E  and  B  are  perpendicular  to  each  other. 

The  concept  of  polarization  of  a  plane  wave  is  related  to  the  nature  of  the  vector  £.  Suppose 
that  c i  and  €2  6  M3  are  orthogonal  real  unit  vectors  which  are  hoth  orthogonal  to  k ,  the 
direction  of  propagation.  Since  £  is  orthogonal  to  k ,  it  is  a  complex  linear  combination  of  €1 
and  €2 

£  =  E\E\  +  25/2  €2 1  (5.7) 
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where  E\  and  E2  arc  complex  numbers.  If  the  complex  phases  of  E\  and  E2  are  equal,  so  that 
Ei  =  \Ei\e1*  and  E2  =  then  the  real  electric  field  is 


E(x ,  t)  =  \/\E\\2  +  \E2\2  e  cos(fc  •  x  -  uit  4-  <£),  (5.8) 

where  e  is  the  unit  vector 

e  =  (|£,|e,  +  \E2\e2)  / ( %/ 1  jEi  |2  +  |E2|2).  (5.9) 

We  see  that  in  this  case  the  clectrie  field  oscillates  between  ±\/\E\  |2  +  | £7^ | 2  e  always  re- 
maimng  parallel  to  e.  In  this  ease  the  plane  wave  is  said  to  be  linearly  polarized  in  direction 
e  In  situations  where  the  plane  wave  propagates  parallel  to  the  surface  to  the  Harth,  waves  lin¬ 
early  polarized  in  a  direction  perpendicular  to  the  Earth’s  surface  are  called  vertically  polarized 
and  those  linearly  polarized  in  a  direction  parallel  to  the  Earth’s  surface  are  called  horizontally 
polarized. 

Another  important  type  of  polarization  occurs  when  |Ei|  =  IE2I  and  the  phase  of  E\  is 
different  from  that  of  E2  by  ±7r/2  ,  so  that  E2  —  ±iE\.  The  real  electric  field  then  has  the 
form 


E(x ,  t)  =  \E\  |  (ei  cos(fc  •  x  —  ut  +  <j>)  T  €2 sinffc  •  x  —  ut  +  <f>)) .  (5. 10) 

Notice  that  here  E  just  rotates  around  a  circle  of  radius  |  E\  |  with  angular  frequency  te  radi¬ 
ans/sec.  Such  waves  arc  called  circularly  polarized.  The  two  possibilities  T  in  Equation  (5. 10) 
correspond  to  the  E  rotating  in  opposite  directions  around  the  circle.  These  two  possibili¬ 
ties  arc  referred  to  as  being  left-handed  and  right-handed  circularly  polarized  depending  on 
the  convention  chosen.  Finally,  we  note  that  a  plane  wave  which  is  not  linearly  or  circularly 
polarized  is  said  to  he  elliptically  polarized. 


5.2  Polarimetric  Transmission  Systems 

The  vector  potential  for  a  sinusoidally  oscillating  source  (antenna),  that  is,  a  source  consisting 
of  localized  charges  and  current  with  charge  and  current  densities  of  the  form 

,>{x,t)  =  p(x)e'“l  (5.11) 

J(x.t)  =  J(x)c1^  (5  12) 

which  is  small  compared  to  the  wavelength  is 


i{kr-uit) 

A(x,t)  =  -ikp- -  (5.13) 

r 

where  for  free  space  k  =  u/c  and  r  is  the  distance  between  the  transmitter  and  the  receiver. 
Here  p  is  a  possibly  complex  valued  electric  dipole  moment  of  the  source 

p  —  j  x'p(x')dxf.  (5.14) 
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The  corresponding  cleetne  and  magnetic  fields  are 


B  =  V  x  A  =  k2(n  x  p) — 


2  € 
E  =  -V  x  B  =  - 


(^k2(n 


x  p)  x  n  +  (3n(n  p) 


-p)(^ 


(5.15) 

7))' <5161 


where  x  -  rn.  In  the  far  field  or  radiation  zone  (At  1)  this  reduces  to  the  limiting  form 


i(kr- u/1) 

B  =  Ar(n  x  p) -  (5.17) 

r 

ei(kr-^t) 

E  —  k2  ((n  x  p)  x  n) - .  (5.18) 

i' 


The  electric  field  can  be  rewritten  as 


i(fer— U?l) 

E  =  k2  {p  —  (p  *  n)n)  — 

r 


(5.19) 


Notice  that  p  -  (p  n)n  is  the  projection  of  the  dipole  moment  p  onto  the  plane  perpendicular 
to  n.  At  a  point  x  =  rn  in  space  the  electric  and  magnetic  fields  look  locally  like  an  outgoing 
plane  wave  with  polarization  vector  £  —  p  (p  n)n. 

For  antennas  which  are  small  compared  to  a  wavelength,  which  is  often  the  case  in  mobile 
wireless  communication  systems,  we  can  imagine  having  three  feeds  into  the  antenna  config¬ 
ured  in  such  a  way  that  we  have  complete  control  over  the  dipole  moment  p  of  the  antenna.  For 
a  larger  antennas,  or  for  concreteness  in  the  ease  of  small  antennas,  we  consider  so-called  triad 
antennas.  A  triad  antenna  is  composed  of  three  orthogonal  dipoles  oriented  along  euclidean 
directions  as  shown  in  Figure  5.1.  At  the  transmitter  a  triad  antenna  can  generate  an  arbitrary 
oscillating  dipole  moment  and  consequently  an  arbitrary  polarized  electric  field  at  the  receiver, 
subject  only  to  the  constraints  imposed  by  the  physics  of  the  electromagnetic  field,  A  triad 
antenna  at  the  receiver  allows  the  receiver  to  measure  the  electric  field  as  shown  in  Figure  5.2. 

The  use  of  triad  antennas  allows  us  to  think  of  this  MIMO  system  in  an  unusual  way.  Instead 
of  thinking  of  the  individual  component  antennas  and  the  symbols  coded  onto  them  at  the 
transmitter  and  read  off  at  the  receiver,  we  can  think  directly  of  coding  onto  a  physical  dipole 
moment  vector  at  the  transmitter  and  measuring  the  resulting  electric  field  at  the  receiver 
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Figure  5. 1:  Schema  tic  of  a  triad  antenna 


Figure  5.2.  A  triad  antenna  system.  Electric  field  (E)  at  the  receiver  is  perpendicular  to  the 
direction  of  propagation 

A  special  ease  of  the  triad  antenna  is  the  more  usual  dual  polarimctric  antenna.  A  dual- 
polari7ed  antenna  consists  of  two  eo-loeated  polarized  antennas,  typically  horizontal/vertical 
(0°/90°)  or  slanted  (+45°/  —  45°)  as  shown  in  Figure  5.3.  The  signals,  say  Xi  and  £2, 
arc  transmitted  on  the  two  different  polarizations,  and  r  1  and  r2  arc  the  signals  received  on 
the  corresponding  polarization.  Note  that,  although  there  is  one  physical  transmit  and  one 
physical  receive  antenna,  the  underlying  channel  is  a  two  input  two-output  channel  since  each 
polarization  mode  is  treated  as  a  separate  physical  channel. 

90 


Figure  5.3 :  Schematic  of  a  dual -polarized  antenna 


5.3  Space-Polarization-Time  System:  Alamouti  Signaling 

Nabar  ct  al  [58]  have  studied  the  performance  of  a  system  with  one  dual-polarized  transmit 
and  one  dual-polarized  received  antenna  with  Alamouti  signaling  and  spatial  multiplexing. 
Deng  et  al.  [60]  have  extended  these  results  to  the  case  of  two  dual-polarized  transmit  and 
one  dual-polarized  received  antenna  with  Alamouti  signaling  and  proposed  a  particular  hybrid 
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transmission  scheme  In  this  section  we  further  analyze  Deng  et  al's  system  and  investigate 
the  performance  of  two  Alamouti  multiplexing  in  the  presence  of  polarization  diversity  and 
compare  with  the  performance  achieved  using  a  uni -polarized  scheme  and  the  hybrid  scheme  of 
Deng  et  al.  We  analyze  and  predict  the  performance  of  the  systems  with  different  transmission 
schemes  by  applying  the  parameter  A  developed  in  Chapter  2.  The  parameter  A  measures  the 
“ quaternionic  angle "  between  the  two  Alamouti  channels. 


5.4  MIMO  Channel  Model 


Consider  a  MIMO  wireless  communication  system  with  nt  transmit  and  n,  receive  antennas. 
The  channel  is  assumed  to  have  quasi-statie  fading  represented  hy  a  nr  x  i}t  matrix  H  with 
entries  which  are  complex  Gaussian  random  variables.  The  channel  matrix  can  be  decomposed 
into  the  sum  of  an  average  (LOS)  component  and  a  variable  (scattered)  component  as 

h = /,  1''a-h  -  v  i  J  a-  <52»> 


where  £{H}  =  yj  —77  H  is  the  channel  mean.  The  factor  yJjjx  and  are  energy 

normalization  faetors  and  arc  related  to  the  Riecan  /\  factor,  with  K  being  referred  to  as 
the  A'— factor.  Note  that  K  =  0  corresponds  to  the  case  of  a  pure  Rayleigh  fading  channel 
and  K  —  00  corresponds  to  the  case  of  a  pure  LOS  channel  K  =  10  is  used  to  represent 
the  typical  suburhan  environment  of  a  personal  communication  system  (PCS),  with  not  very 
high  huilding  density  and  partial  LOS  [60).  The  elements  of  the  matrix  H  are  zero-mean 
circularly  symmetric  complex  Gaussian  random  vanahles  and  unit  vananee.  Rt  —  R)  2  Rlf  2f 
and  Rr  =  R\  R)  2  are  transmit  and  receive  correlation  matrices  The  LOS  component  H 
is  modelled  as 

H  =  {3a,a]  (5.21) 

where  at  and  ar  are  transmit  and  receive  steering  vectors  and  ft  is  a  complex  gain.  If  the 
transmit  and  receive  antennas  are  uniform  linear  arrays  with  spaeing  d  then 


^  —  ij  f  ~  tnc^cos  l7r(nr  -  1  )rfcos  0i 


ar  =  (l,e 


nrd  cos  6, 


-i7r(n,  —  l)dcos^r  yi 


(5.22) 


are  the  transmit  and  receive  steering  vectors  respectively.  We  take  H  to  be  normalized  such 
that 


Tr(H*H)  =  mil  r, 


(5  23) 


so  that  |/i|  =  1. 


5.5  Dual-Polarized  Channel  Model 


Consider  a  system  with  one  dual-polarized  transmit  and  one  dual-polarized  receive  antenna. 
The  2x2  channel  matrix 


H  = 


(hvv  h  hv 
hvh  /ihh 


(5.24) 
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is  a  polarization  matrix  where  /iw  and  hbb  represent  path  gain  from  vertical  transmit  to  vertical 
receive  antenna  and  from  the  horizontal  transmit  to  the  horizontal  receive  antenna,  respectively. 
Similarly,  hbV  and  hv h  represent  path  gain  from  the  horizontal  transmit  to  the  vertical  receive 
antenna  and  vice  versa.  The  elements  of  the  matrix  H,  which  are  denoted  as  —  v,  h)  are 

zero-mean  circularly  symmetric  complex  Gaussian  random  variables  whose  variances  depend 
on  the  propagation  conditions  and  the  antenna  characteristics  We  use  the  model  given  in  [58], 
that  is 


£{|/iw|2}  =  E{|/)hh|2}  =  1 

£{IM2}  =  £{IM2}  =  7. 


(5.25) 


where  0  <  7  <  1  describes  the  cross  polarization  discrimination  (XPD)  or  separation  of 
orthogonal  polarization  for  the  variable  component  of  the  channel.  7  is  a  composite  of  the 
properties  of  the  antennas  and  the  scattering  environment  (coupling  between  elements  due  to 
scattering).  Good  discrimination  of  orthogonal  polarization  corresponds  to  a  small  value  of  7 
and  vice  versa.  The  elements  of  the  matrix  H,  which  are  denoted  as  (i,j  =  vh),  are  fixed 
complex  numbers  satisfying 


IM2  =  |/ihhl2  =  1,  IM2  =  IM2  =  7/' 


(5.26) 


where  0  <  7/  <  1  is  related  to  the  XPD  for  the  LOS  component.  For  pure  LOS  conditions, 
7/  is  solely  a  function  of  the  antennas'  ability  to  separate  the  orthogonal  polarization. 

It  was  stated  in  [58]  that  the  experimental  data  collected  in  [01]  and  [62]  shows  that  the 
elements  of  H  are,  in  general,  correlated.  The  correlation  coefficients  are  therefore  defined  as 


t  _  £{^wfehv}  _  ft  {fry  h/»hh}  r  _  ft{ftw/ly*h}  _  ftfohvfthh} 

\/7  v/7?  ’  v/7  v/77 


(5.27) 


where  t  is  referred  to  as  the  transmit  correlation  coefficient,  and  r  is  the  receive  correlation 
coefficient  The  expen  mental  data  [58]  also  reveals  that  the  correlation  between  the  diagonal 
element  matrix  /?w  and  h uu  and  the  off-diagonal  elements  /?vh  and  hhv  is  typically  very  small. 
Therefore,  for  simplicity,  we  assume  that 


EfhwhJh}  =  £{hhv/Ch}  =  0. 


(5.28) 


Consider  a  two  user  system  with  two  dual-polarized  transmit  and  one  dual-polarized  receive 
antenna,  where  each  user  employs  the  Alamouti  ,STBC  over  a  dual -polarized  transmit  antenna 
The  equivalent  polarization  matrices  from  the  first  user  to  the  first  and  second  receive  antenna, 
can  be  represented  as 


Hi 


h  w  hbv 

^  vh 


and  H2  = 


/lvh  h-hh 

h-hh  ^  vh 


(5.29) 


Similarly,  the  equivalent  channel  matrices  for  the  second  user  to  the  first  and  second  receive 
antenna  are 

<7hv\  (  fjv  h 

Qx  =  I  .  .  ,  and  g2  =  .  .  1  (5.30) 


9\th  9w 


9bb  Gvb 
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We  consider  three  different  transmission  schemes.  The  first  is  a  uni-polarized  system  as 
described  in  Section  2,  i.e.  four  vertically  polarized  transmit  and  two  vertically  polarized 
receive  antennas,  all  spatially  separated  as  shown  in  Figure  5.4. 
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Figure  5.4:  Scheme  l :  Uni -polarized 


The  LOS  components  of  the  equivalent  channel  matrices  arc  modelled  as 


Ml  1 


Mi\  _ 

,n2 


ll2  H3  \  _ 

"M3  M 2  I 


ftA  fl5  \  _ 

_  .  , 02  = 
M 5  Fa 


Me  F7 
M?  Me 


(5.31) 


where  •  • .  ,M7  are  arbitrarily  chosen  unimodular  complex  numbers  that  represent  the 

different  LOS  paths.  Note  that  in  this  case  7  —  1. 

The  second  scheme  consists  of  two  dual-polarized  transmit  and  one  dual-polarized  receive 
antenna  and  employs  Alamouti  STBC  on  each  pair  of  dual-polarized  transmit  antennas  as- 
shown  in  Figure  5.5. 
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Figure  5.5:  Scheme  2: Dual-polarized 


This  scheme  will  be  referred  to  as  the  dual-polarized  system  and  the  equivalent  LOS  com 
ponent  channels  arc  modelled  as 


«i  = 


1  ^  1 


Vt7  1 


-1  y/Tr 


s/Yiv  /»’  )  \~n' s/vf‘\ 


/i  fi 


(5.32) 
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where  p  is  again  a  unimodular  eomplex  number.  Under  the  assumption  that  the  receiver  is 
in  the  far  field  of  the  transmit  antenna,  then  p  =  oxp(  — A),  where  A  is  the  carrier 
wavelength,  n  is  the  unit  vector  in  the  direction  of  the  receiver  and  A  is  the  position  of  the 
seeond  antenna  relative  to  the  first  antenna.  Furthermore,  since  the  two  orthogonal  components 
of  each  of  the  dual  polarized  antennas  are  eo-loeated,  their  relative  path  difference  is  zero. 

The  third  scheme  is  that  proposed  by  Deng  et  al.  [60].  It  also  consists  of  two  dual-polarized 
transmit  antennas  and  one  dual-polarized  receive  antenna,  but  employs  Alamouti  STBC  on  the 
same  polarization  across  the  two  transmit  antennas  as  shown  in  figure  5.6. 


Vertical  of  1“  Tx 


Horizontal  of  1st  Tx 
Horizontal  of  2m*  Tx 


Figure  5.6:  Scheme  3: Dual-polarized  hybrid 


In  this  ease  the  equivalent  LOS  channel  matrices  are  modelled  as  follows 


H  i- 


( 


i 

-m* 


v/7j  \A7>\  r  _(  s/V  \A7>\  /  1  »\ 

v/w  Su  )  ^-y-sTuv-'  V-'1’  V 


where  p  is  as  defined  above. 


(5.33) 


5.6  Analysis  and  Simulation 

We  now  analyse  the  three  transmission  sehemes  desenbed  above  by  applying  the  parameter  A 
obtained  in  Chapter  2  We  consider  both  the  Rayleigh  fading  channel  and  the  Rieean  channel 
where  A  =  10  In  the  simulation  results  reported  here,  the  data  symbols  were  taken  from 
an  8-PSK  constellation  and  deeoded  with  the  algorithm  given  in  (2.42)  and  (2.49).  For  the 
polarization  sehemes  we  set  the  XPD  parameters  to  7  =  0.4  and  7/  =  0.3  and  for  correlated 
channels,  the  transmit  and  receive  correlation  coefficients,  as  defined  in  (5.27),  were  ehosen  to 
be  t  =  0.5  and  r  =  0.3. 

Figure  5.7  illustrates  how  the  angles  between  the  channel  matrices  |A|2  are  distributed  for 
an  uneorrelated  Rayleigh  fading  channel,  for  the  three  different  transmission  sehemes.  Fig¬ 
ure  5.8  shows  the  symbol  error  rate  ohtained  under  the  same  channel  conditions.  Figure  5.9 
and  Figure  5  10  show  the  corresponding  results  for  the  ease  of  an  uneorrelated  Rieean  channel 
when  A  =  10  11  is  elear  that  the  presence  of  the  LOS  component  eauses  the  means  of  the  |  A| 2 
distributions  to  shift  and  the  distributions  to  heeome  more  concentrated  around  the  mean.  The 
effect  is  quite  pronounced  even  for  relatively  small  values  of  K.  This  behavior  holds  generally 
and  leads  to  the  conclusion  that,  when  present,  the  angle  hetween  LOS  paths  is  the  dominant 
factor  in  the  performance  of  various  sehemes. 
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Figure  5.1 1  shows  the  |A|2  distributions  for  three  different  LOS  paths  for  correlated  chan 
nels.  Figure  5  12  shows  the  corresponding  results  with  the  XPD  parameter  set  to  7/  —  0.1 
rather  than  0,3.  Note  that  each  row  represents  a  different  LOS  path  and  each  column  rep¬ 
resents  different  transmission  schemes.  It  is  evident  that  the  |A|2  distributions  for  the  uni 
polarized  and  dual-polarized  schemes  have  considerable  variability  as  the  relative  phases  in 
the  LOS  components  vary.  However,  what  stands  out  is  the  stability  of  the  |A|2  distribution  for 
the  dual-polarized  hybrid  transmission  scheme  The  distribution,  and  hence  the  performance 
of  this  scheme,  is  insensitive  to  the  relative  phases  of  the  LOS  components  of  the  channel 
We  can  understand  this  behavior  by  considering  the  angles  between  the  LOS  components  of 
the  channels  which  we  denote  by  Alos.  From  (2.33)  and  (5.20),  we  have 


A  LOS  = 


wg 

w 


(5.34) 


For  the  dual-polarized  hybrid  system 


I  Alos  |  = 


‘Wv 

1+7/' 


(5.35) 


and  for  the  dual-polarized  system  we  have 


|Alos|  =  «/cos2<£  +  ^f—rs\n2<t >, 


(1  +7/)2 


(5.36) 


where  <f>  —  Z/x. 

We  immediately  note  from  (5.35),  that  |Alos|  for  the  dual-polarized  hybrid  scheme  depends 
only  on  7/.  This  explains  why  the  dual-polarized  hybrid  scheme  is  robust  to  changes  in  the 
relative  path  lengths  between  the  two  dual  polarized  transmitters,  as  observed  in  Figures  5.1 1 
and  5. 12. 

We  note  that  for  the  dual-polarized  scheme 

,->v  '  <  |A,  os|  <  1-  (5.37) 

1  +  7/ 

We  have  demonstrated  that  the  performance  parameter  A  obtained  in  Chapter  2  can  be  used  to 
analyze  and  predict  the  performance  of  different  transmission  schemes  involving  polarization 
diversity  and  the  Alamouti  codes.  The  results  show  that  the  performance  is  dominated  by  the 
“angle”  between  the  component  Alamouti  channels.  In  designing  transmission  schemes  that 
involve  multiplexing  Alamouti  coded  transmissions,  one  should  aim  for  multiplexing  schemes 
which  keep  the  channel  volume  away  from  zero  (see  Section  2.8  on  channel  volume).  This  is 
the  case  for  the  dual-polarized  hybrid  scheme  analysed  here,  which  has  a  normalised  channel 
volume  dominated  by  the  LOS  channel  volume 

vios  =  (1  -7/)2/(l  +7/)2.  (5.38) 
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Figure  5.7:  Distribution  of\\\2for  three  ciijferent  transmission  schemes  over  an  uncorrelated 
Rayleigh  channel 
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Figure  5.8:  Performance  of  three  different  transmission  schemes  over  an  uncorrelated  Rayleigh 
channel 
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Figure  5. 9:  Distribution  of\X\2for  different  transmission  schemes  over  an  uncorrelated  Rice  an 
channel 
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Figure  5. 1 0:  Performance  of  three  different  transmission  schemes  over  an  uncorrelated  Ricean 
channel 
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figure  5.1 1:  Distribution  of\X\2  for  three  different  LOSs,  7/  —  0.3 
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Figure  5.12:  Distribution  of\X\2for  three  different  LOSs,  77  =  0.1 
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5.7  Polarization-Time  System:  The  Golden  Code 


In  the  last  section,  we  investigated  the  performance  of  three  different  transmission  schemes 
using  Alamouti  STBCs  across  a  pair  of  dual-polarized  antennas  at  the  transmitter.  The  results 
indicate  that  the  performance  is  dominated  by  the  LOS  component  and  that  polarization  diver¬ 
sity  along  with  an  appropriate  transmission  scheme  can  provide  stahility  of  performance  across 
changes  in  environmental  conditions. 

Typical  of  uni-polarized  space-time  systems,  the  performance  of  the  Golden  Code  degrades 
considerably  as  the  LOS  component  increases  due  to  the  loss  of  diversity.  Figure  5.13  shows 
the  simulated  performance  of  4-QAM  and  16-QAM  modulation  in  various  channel  conditions 
(pure  scattering  K  —  0,  Riccan  K  —  10  and  pure  LOS  I\  =  oo  in  uncorrelated  channel). 
These  simulation  results  show  that  the  performance  degrades  as  the  K—  factor  increases.  In 
this  chapter  we  consider  the  effect  of  LOS  on  the  performance  of  the  Golden  Code  and  intro¬ 
duce  a  possible  remedy  w  ith  the  use  of  dual-polarized  antennas.  We  analyse  the  performance 
of  the  Golden  Code  in  terms  of  angle  between  the  channels  corresponding  to  the  two  receivers, 
which  allows  us  to  predict  the  performance  without  necessarily  needing  to  resort  to  simula¬ 
tions  Analysis  and  simulation  results  show  that  with  the  introduction  of  polarization  diversity 
the  performance  of  the  Golden  Code  can  be  made  consistently  good  across  both  rich  scattering 
and  LOS  conditions 


Figure  5. 13:  Performance  of  the  Golden  Code  for  uncorrelated  uni-polarized  channels  from 
pure  scattering  to  pure  LOS 


5.8  The  Effect  of  LOS  on  the  Performance 

Codewords  in  the  Golden  Code  take  the  form 

fa  0\  (  Xi  +  xiT  X3  +  X4 r 

X  - 


0  a)  \i(x  3  T  X4p)  X  i  -f  X  2/i 


(5.39) 
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where  a  =  (1  +  */*)/\/5,  a  =  (1  +  ir)/\/ 5  and  {xj^  ,  G  C  C  Z[z]  are  transmitted  sym¬ 
bols,  where  C  is  a  signal  constellation  taken  to  be  2m  —  QAM  and  m  bits  per  symbol.  The 
parameters  r  and  p  are  the  real  roots  of  the  polynomial  x2  —  x  —  1,  that  is,  the  Golden  ratio 
r  =  (1  +  \/5)/2  and  its  algebraic  conjugate  p  —  —  1/r  =  (1  —  \/5)/2  The  diagonal  matrix 
diag(a  rv)  serves  to  equalize  transmitted  signal  power  across  the  two  transmit  antennas. 

Let  rj  =  (r 1 1 ,  r\2)  and  r2  =  (r2i.  r22)  be  the  two  received  signal  vectors  at  the  first  and 
the  second  receiver  and  the  components  are  the  signals  received  over  two  consecutive  time 
slots.  For  convenience,  in  what  follows  we  will  use  the  rescaled  channel  gains,  h\3  —  o htJ 
and  h.2j  =  ah2j  c  =  (x  1,0:3)  and  s  =  (x2,X4),  so  the  received  signals  can  be  written  as 


where 


r  1  =  c.h  -f  sh  +  ni 
t 2  ~  eg  +  sg  +  n2 


h  — 


9  = 


(  hnT  h2ifJ. 
\ih2ifi  h\iT 


(h\2T  fi22fl 

ih22p  J112T 


(5.40) 


(541) 


where  hij  is  the  channel  gain  from  transmit  antenna  1  to  the  receive  antenna  j. 

Given  that  the  channel  gains  arc  known  at  the  receivers  and  each  symbol  is  transmitted 
with  equal  probability,  optimal  decoding  is  provided  by  the  maximum  a  posterior/maximum 
likelihood  (MAP/ML)  estimate  as  follows.  Let  r  —  (r\,r2)  W  =  (/i, h)T  and  Q  =  ( g,g )7 
and  rewrite  (5.40)  as 

r  —  xA  4-  n,  (5.42) 


where  A  —  (W,  Q)  and  x  —  (c,  5).  The  likelihood  function  of  x  given  ihe  received  signal  r 
is  given  by 

p(r\x)  (x  exp  -j-  -(x  —  x)AA  (x  —  (5  43) 

where 

x  -  rA^AA  )~V  (5.44) 


Taking  the  prior  distribution  of  x  to  be  uniform  on  the  constellation  C ,  we  obtain  the  MAP/ML 
estimate: 

x  —  arginaxp(r|x).  (5.45) 

16  c4 


The  performance  of  (5  45)  in  decoding  x  is  determined  by  the  determinant  of  AAt  or  equiva¬ 
lently  A*  A.  vSincc 


||WH*  la  H'Q  \ 
WGWlh) 


(5.46) 


where 


HW|i2  =  (l  +  r2)|/li.|2  +  (l  +  /i2)|/t2ir, 
II<?II2T-(1  +  T2)|/112|2  +  (1  +  M2)|/122|2, 
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wc  have 

det(A*A)  =  ||*|£|K?||t(l  |A|2)2  (5.48) 


where 

H'Q 

IIW||TII«?IIT' 


(5.49) 


The  parameter  A  is  an  inner  product  of  two  unit  cyclotomic  vectors  and  it  represents  the  differ¬ 
ences  (angle)  between  the  channel  components  H  and  Q.  A  similar  measure  to  this  parameter 
which  represents  the  angle  between  two  Alamouti  users  was  given  in  Section  2.4,  [48]  and  was 
used  to  predict  the  performance  of  various  Alamouti  transmission  schemes  involving  dual- 
polarized  antennas  (see  [40]  for  details).  Good  performance  corresponds  to  small  values  of 
| A| 2  and  vice  versa. 

Consider  the  capacity  of  the  equivalent  channel  matrix  A  which  is  given  by  [9], 


C(A)  =  -/J 


2E 


log  (  det  (  I4  +  A  A 


(5.50) 


where  E9  is  an  average  transmitted  signal  power  and  e*  are  eigenvalues  of  A1  A  with  multi¬ 
plicity  2  and  are 


ei,2 


\ml  WQfr 


2 


1  ± 


\\n\\l  +  \\G\\l{1  |A,) 


(5.51) 


This  implies  that  the  capacity  of  the  channel  depends  on  the  parameter  A.  If  |A|2  =  0  then  we 
have 

e,  =  ||H||2  and  r2  =  ||0||2  (5.52) 


and  the  corresponding  capacity 

C|A|2 :  o(A)  =  2 E 


log 


If  | A | 2  =  I  then  the  channel  is  rank  deficient,  that  is 

Cl  =  ||W||2  +  ||(7||2  and  e2  =  0, 
and  the  corresponding  capacity  is 

q^=i(A)  =  2 E  log  (l  +  ~  (IIHII2  +  ll^ll2 


(5.53) 


(5.54) 


(5.55) 


Consider  the  case  of  a  pure  scattering  channel,  i.c.  K  =  0,  then  the  values  of  |A|2  is  uniformly 
distributed  between  [0  1]  and  can  be  explicitly  written  as 


2  l(M)r|2 

K  0  —  -  *  jj2  ||  m2 

n ||  ||e|| 


(5.56) 
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Figure  5. 14:  Distribution  of\\\2for  uncorrelated  uni-polarized  channels 


where 

(h,g)T  =  (1  4  T2)k\xhi2  +  (14*  /i2)/?2l/i22, 
|  n\  —  (1  +  T2)|/?.i  t  |2  +  (1  4-  /X.2 ) |/i-21 1 2 9 

1^11  =  (1  +  7"2  )  |/l  1  2  |2  4-  (1  4  /£2  )  I  /?  22  I  2  • 


It  is  also  clear  from  (5.56)  that  in  the  pure  scattering  environment,  the  channel  becomes  rank 
deficient,  i.e  |A|x=o  =  1,  only  if  hnh.22  =  ^21^12  or  equivalently  det(H)  =  0.  As  K 
increases  H  dominates  the  behaviour  of  the  system  and  when  K  =  00,  which  corresponds  to 
pure  LOS  channels,  we  have 


/  ]  e%n  cos  or  \ 


H  - 


•  I  7T  COS  Or 

ic  1 


1 


lie 


ynelK,:os0‘-  T 


Q  = 


ITT  COS  Oh  ^,17t(cOS  Oy—COS  Oft)  \ 


ITT  (cos  er  —  cos  On) 


\ific 


re 

a tt(cos  O  f  —  cos  On) 


pe 


e  n 

i tt( cos  0-f  -cos  On  ) 


ITT  COS  0  n 


Substituting  these  values  into  (2.33),  we  obtain 

[Alios  =  1> 


(5.58) 


i.e.  the  channel  is  rank-deficient  hence  resulting  in  the  performance  degradation. 

Figure  5.13  shows  the  results  of  simulated  performance  of  the  Golden  Code  with  4-QAM 
and  16-QAM  modulation  in  various  channel  conditions  (pure  scattering  (K  =  0),  Ricean 
(K  =  10)  and  pure  LOS  {K  =  00)  in  uncorrelated  channels).  Figure  5.14  shows  the  empirical 
distribution  of  |A|2  for  K  =  0,  10  and  100.  Clearly,  as  the  K -factor  increases  the  mean  of 
the  distribution  of  | A[ 2  shifts  toward  |A|2  =  1.  Thus,  as  expected,  the  distribution  of  |A|2 
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indicates  that  the  performance  degrades  as  K  increases.  The  distribution  of  | A | 2  agrees  with 
the  simulation  results. 


5.9  Golden  Code  and  Polarization  Diversity 


Consider  Golden  Code  signaling  using  one  dual-polarized  transmit  and  one  dual-polarized 
receive  antenna  Based  on  the  channel  model  given  by  Nabar  et  al.[58]  as  described  in  Sec¬ 
tion  5.5,  the  LOS  component  can  be  explicitly  written  as 


1 

Vv  \ 

(yJTt  1  \ 

*v/t7 

1 

n  _ 

*  \TTJ 

T  y/T/fl 

1  3/  — 

s/riT  /< 

\jVrin  T  j 

Substituting  this  into  (5.34),  we  hnd 

\  1 2  _  _ (2  -f  r2  +  /J2)27/ _ 

'  lLOS  [(1  +  r»)  +  (1  +  M2)7/][(l  +  r2)7/  +  (1  +  K2)) ' 


(5.59) 


(5.60) 


If  7/  <  1.  then  |A|£0s  <  1.  As  (5.60)  shows,  the  antenna’s  cross-polarization  diserim 
ination  (7/)  will  determine  the  mean  performance  of  the  system.  For  a  pure  LOS  condition, 
the  performance  depends  solely  on  the  antenna's  ability  to  separate  the  orthogonal  polariza¬ 
tion  Therefore  as  the  K—  factor  increases  the  effect  of  polarization  diversity  is  to  reduce  the 
probability  that  the  channel  will  be  rank-defieient,  thus  resulting  in  a  diversity  gain 

We  ean  prediet  the  performance  of  the  Golden  Code  in  the  presence  of  polarization  diversity 
based  on  the  distribution  of  the  parameter  |  A | 2 .  Figure  5.15  shows  the  distribution  of  |A|  2  for 
K  =  0  10,  100  when  dual-polarized  antennas  are  used.  Compared  to  the  ease  of  uni-polarized 
antennas  in  Figure  5. 14,  the  mean  of  the  |A|2  for  K  >  0  is  now  shifted  away  from  the  value  1 
The  simulation  results  in  Figure  5.16  demonstrate  the  benefit  of  dual-polarized  antennas  in  the 
high  K  environments  as  predieted  in  Figure  5.15  Note  that  for  pure  LOS  conditions  (A  =  00) 
the  channel  is  fixed,  therefore  the  probability  of  error  decays  exponentially  with  SNR  and  so 
no  diversity  order  ean  be  defined. 

Figure  5.17  shows  the  distribution  of  |  A | 2  in  a  correlated  environment,  with  t  =  0.5,  r  = 
0.3.  We  ean  see  from  the  distribution  that  for  a  pure  Rayleigh  channel  ( K  =  0),  |A  2  is  no 
longer  uniformly  distributed  but  biased  toward  1.  Hence  the  performance  will  degrade  com¬ 
pared  to  an  un  correlated  environment.  Figure  5. 1 8  demonstrates  that,  for  K  >  0,  the  better  the 
separation  of  orthogonal  polarization  (smaller  7/),  the  better  the  improvement  in  performance. 


5.10  Summary 

We  investigated  the  performance  of  three  different  transmission  schemes  for  using  Alamouti 
STBCs  across  a  pair  of  spatially  separated  dual-polarized  transmit  antennas.  The  results  indi¬ 
cate  that  the  performance  is  dominated  by  the  LOS  component  and  that  polarization  diversity 
along  with  an  appropriate  transmission  scheme  ean  provide  stability  of  performance  across 
ehanges  in  environmental  conditions  We  went  on  to  investigate  the  performance  of  the  Golden 
Code  using  polarization  in  various  environmental  conditions  from  rich  scattering  to  pure  line 
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of  sight.  Our  analysis  leads  to  a  performance  parameter  which  represents  the  angle  between  the 
channels  to  the  two  receivers  and  can  be  used  to  predict  performance.  We  have  analysed  how 
the  performance  of  the  Golden  Code  degrades  as  the  environmental  K -factor  increases.  When 
the  Golden  Code  is  used  to  code  across  orthogonal  polarizations  of  a  dual -polarized  transmit 
antenna,  rather  than  across  spatially  separated  antennas,  we  have  shown  that  polarization  di¬ 
versity  can  lead  to  consistently  good  performance  across  environmental  conditions  from  rich 
scattering  to  pure  line  of  sight  The  simulation  results  agree  with  our  analysis.  We  note  that 
the  approach  is  not  restricted  to  the  Golden  Code  but  can  be  applied  to  most  STBCs. 


Figure  5. 15:  Distribution  of\\\2for  uncorrelated  channels  with  polarized  antennas ,  7/  =  0.4 


Figure  5. 16:  Performance  of  the  Golden  Code  for  uncorrelated  channels  with  polarized  anten¬ 
nas,  7/  =  0.4 
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Figure  5. 17:  Distribution  of  \\\2  for  correlated  c  hannels .  t  =  0  5,  r  =  0 .3,7/  =  0.4 


i*-'2 


Figure  5.18:  Distribution  of\X\2  for  uncorrelated  channels ,  7/  =0.3 
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Chapter  6:  Fully  Polarimetric  MIMO 


In  the  previous  chapter  we  investigated  the  performance  of  a  MIMO  system  using  the  Golden 
Code  with  a  dual-polarized  antenna  system  and  found  that  the  performance  can  be  made  con¬ 
sistently  good  across  environmental  conditions  from  rich  scattering  to  pure  LOS.  However,  we 
assumed  that  the  planes  of  the  dual  polarized  transmit  and  receive  antennas  remained  aligned 
The  simulation  results  for  dual-polarized  systems  show  that  the  performance  degrades  consid 
crably  if  we  allowed  three  dimensional  (3-D)  rotation  of  the  transmit  and/or  receive  antennas 
If  we  arc  to  profitably  use  the  performance  stability  introduced  by  polarimetric  antennas  for 
mobile  terminals  then  we  need  to  address  the  stahility  in  transmission  and/or  reception  under 
the  3-D  rotation  of  the  antennas  in  the  MIMO  system  To  this  end  we  investigate  the  use  of 
triad  antenna  at  the  transmitter  and  receiver.  A  triad  is  composed  of  three  orthogonal  dipoles 
oriented  along  orthogonal  directions  as  described  in  Section  5.2  We  derive  the  channel  model 
and  show  that  the  capacity  is  invariant  under  arbitrary  rotation  of  the  transmitter  and/or  receiver, 
since  the  channel  matrix  is  always  a  3  x  3  matrix  of  rank  2.  We  propose  a  3  x  3  STBC  suitable 
for  this  system  which  achieves  full-rate  and  furthermore  has  a  fast  decoding  algorithm  We 
give  the  simulation  results  for  the  performance  of  the  triad  system  using  the  proposed  3x3 
STBC  with  and  without  antenna  rotation  We  then  compare  the  results  with  a  dual-polarized 
system  transmitting  Golden  Code  codewords. 

6.1  System  Channel  Models 

Consider  a  MIMO  wireless  communication  system  with  three  transmit  and  three  receive  an¬ 
tennas.  We  assume  a  quasi-static  fading  channel,  represented  by  a  3  x  3  matrix  H  with  entries 
which  are  complex  Gaussian  random  variables.  For  each  symbol  interval,  p  =  (pi ,  P3)1 

denotes  a  vector  of  transmitted  complex  symhols  The  received  signal  vector  y  =  (y\ ,  ?/2.  1/3)1 
can  be  expressed  as 


y  =  Hp  +  w, 

where  w  denotes  complex  Gaussian  noise  with  covariance  E{ww *}  —  <J2I$. 


(6.1) 


The  channel  matrix  H  can  be  written  as  the  sum  of  an  average  (LOS)  component  H  and  a 
variable  Rayleigh  fading  component  H 


(6.2) 


where  £{H}  =  y  /\/(l  +  K)  H  is  the  channel  mean.  H  is  a  complex  Gaussian  random 
matrix  with  independent  zero-mean  and  unit  variance  entries,  i.e.  ht]  e  CA^(0, 1).  H  and 
H  are  weighted  hy  the  Rieean  K -factor.  Note  that  h  =  0  corresponds  to  the  case  of  pure 
Rayleigh  channel.  As  K  increases  H  dominates  the  hchavior  of  the  system  and  K  =  00 
corresponds  to  a  pure  LOS  channel.  Rt  =  R)  2 R\  /2t  and  R,  =  /Z,1  2/ZT1/2t  arc  the  transmit 
and  receive  correlation  matrices  respectively 

Consider  a  system  with  a  triad  antenna  at  both  the  transmitter  and  receiver.  A  triad  is  com¬ 
posed  of  three  orthogonal  dipoles  oriented  along  Euclidean  directions  e\ ,  62.  and  as  shown 
in  Figure  5.1  At  the  transmitter  a  triad  antenna  can  generate  an  arbitrary  dipole  moment  p. 
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As  shown  in  Section  5.2  the  far  field  electric  field  generated  by  the  triad  is 


(6.3) 


The  term  h  -  nnt  is  a  3  x  3  projection  matrix  with  rank  2  and  it  constrains  the  oscilla¬ 
tion  direction  of  the  electric  field  to  be  perpendicular  to  its  propagation  direction  as  shown  in 
Figure  5.2.  In  spherical  coordinates,  the  propagation  direction  can  be  represented  as 


/  sin  0  cos  ^ 


n  —  sin  0. sin  0 
^  cos  0  J 


(6.4) 


where  0  <  0  <  7r,  0  <  0  <  27t.  The  use  of  triad  antennas  allows  us  to  think  of  this  Ml  MO 
system  in  an  unusual  way.  Instead  of  thinking  of  the  individual  component  antennas  and  the 
symbols  coded  onto  them  at  the  transmitter  and  read  off  at  the  receiver,  we  can  think  directly 
of  coding  onto  a  physical  dipole  moment  vector  at  the  transmitter  and  measuring  the  resulting 
electric  field  at  the  receiver  with  the  LOS  channel  given  by  the  matrix 


(6.5) 


It  is  difficult  to  build  co-located  antennas  without  having  some  form  of  mutual  coupling.  To 
model  this  coupling  between  the  three  components  of  the  triad  antenna  we  introduce  mixing 
operators  Mt  for  the  transmitter  and  Mr  for  the  receiver  so  that  the  LOS  channel  becomes 
MlHMt  The  mixing  operators  Mt  and  Mr  are  defined  as 


/  1  P\2  pia\ 


pi  2  1  P23 


Ml  -  Mr  = 


(6  6) 


vW 

yp13  p23  1  ) 


where  0  <  ptJ  <  1  is  a  coupling  between  components  i  and  j  and  is  directly  related  to  the 
cross  coupling  between  the  antenna  elements  comprising  the  triad,  and  N  is  a  normalisation 
term  given  by 


N  —  max(l  -f  pi2  4-  p i ;* ,  1  +■  P12  +  p23»  1  T  P?3  T  pli)- 


(6.7) 


Assuming  that  the  receiver  is  at  position  rn  relative  to  the  transmitter  and  denoting  the  three 
orthogonal  axes  of  the  triad  transmitter  and  receiver  (initially  assumed  to  be  aligned)  by  ci, 
e2,  e3,  the  LOS  component  H  then  consists  of  elements  given  by 


(H)ij  =  ejWej.  i,j-  1,2,3. 


(6.8) 


The  elements  of  H  depend  on  the  propagation  conditions  and  the  antenna  characteristics  and 
satisfy  [58], 


£(M2}  =  1.  and 

E{IM2}  =  £{lM2}  =  %  i,i  =  1,2,3, 
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(6.9) 


where  0  <  7  <  1  describes  attenuated  cross  coupling  for  polarization  multiplexing.  In  a  rich 
scattering  environment  7  tends  to  be  close  to  1  (for  ranges  beyond  1.6  km,  7  is  always  close 
to  1)  [58].  Here  we  assume  that  7  =  1,  i.e  the  scattering  changes  the  polarization  states 
randomly,  henec  in  a  pure  scattering  environment,  a  polarimetne  antenna  system  is  equivalent 
to  a  uni  polarized  antenna  system.  The  correlation  matrices  Rt  and  Rr  are  modelled  as 


/I  f  A 

n  r  r\ 

«<  = 

Of  ,  R,  = 

r  1  r  j 

V  '  V 

r  1  j 

where  t  and  r  arc  transmit  and  receive  correlation  coefficients  respectively. 

For  a  dual-polarized  antenna  consisting  of  co-locatcd  vertically  and  horizontally  polarized 
antennas,  consider  a  system  with  one  dual-polarized  transmit  and  one  dual  polarized  receive 
antenna.  The  LOS  channel  matrix  is  given  by 


Hd  =  /’'Hi' 


(6.11) 


where  R  ~ 


and  H  is  as  defined  in  (6.5).  Similarly,  H  j  is  modelled  as  described  for 


triad  antennas  and  the  correlation  matrices  are  defined  as 


Rt  = 


6.2  The  Effect  of  Antenna  Rotation 


Consider  the  case  of  pure  LOS  ( K  —  oo).  If  the  ehanncl  is  known  at  the  receiver,  then  the 
capacity  is  given  by  [2] 


C(H)  =  Fh 


log2  det  ^/:i  +  ^HHf) 


(6  12) 


where  E,  is  the  average  transmitted  signal  power  From  (6.5)  and  (6.8)  it  follows  that  H  H  is 
a  3  x  3  rank  2  matrix  with  equal  eigenvalues,  and  so  the  eapaeity  (6  12)  can  be  written  as 


6(H) 


Fy 


21og2(l  +  ^A) 


(6.13) 


where  A  is  the  eigenvalue  of  H  H 

Assume  that  the  transmit  and  receive  triad  antennas  have  been  rotated  or  equivalently,  the 
transmitted  and  receive  signals  have  been  rotated.  The  rotated  signal  ean  be  represented  as 

V  =  Rp  -  p,  Re !  +  p2Re2  +  f3Re3.  ^ 

y '  =  Sy  =  2/iSei  +  y2Se2  +  2/3 ^>e3, 

where  S  and  R  are  3  x  3  real  orthogonal  rotation  matrices.  The  elements  of  the  corresponding 
channel  matrix  H  become 

(H')0  =  e\ S' H  Re j,  i,j  =  1,2,3.  (6.15) 
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Since  S  and  R  are  orthogonal,  S^TiR  has  the  same  singular  values  as  H ,  i.e.  H  is  still  a  3  x  3 
rank  2  matrix,  therefore  the  capacity  is  invariant  under  the  rotation  of  the  transmit  and  receive 
antennas, 

C(H')  =  C(H).  (6.16) 

This  is  not  the  ease  when  dual-polarized  antennas  are  employed  at  the  transmitter  and  re¬ 
ceiver.  In  that  ease  the  change  in  orientation  (3-0  rotation)  of  the  dual-polarized  transmit  and 
receive  antennas  will  not  preserve  the  capacity  of  the  rotated  channels.  Moreover,  a  transmitter 
or  receiver  can  he  rotated  such  that  the  channel  matrix  becomes  rank  deficient  (in  fact,  rank  1 ). 
For  example,  consider  a  rotated  channel  matrix, 

,  (e\s'HRe i  e\s'HRe2\ 

e\S^'HIie2 J  (6.17) 

If  a  dual-polarized  receive  antenna  is  rotated  such  that  one  of  the  branches  is  parallel  to  the 
direction  of  propagation  n,  i.e. 


Sex  = 

then  (6. 17)  becomes 

H',  = 

Substituting  (6.18)  into  (6  8),  we  obtain  =  0,  and  we  have  a  rank  1  channel  matrix 
However,  if  the  dual-polarized  receive  antenna  is  rotated  in  a  plane,  then  the  capacity  is 
invariant  under  rotation  [63]  That  is,  under  a  rotation  of  the  receive  antenna,  (6.1 1)  becomes 

Hd  -  /jtSfHP.  (6.20) 

If  PS  =  SP  (i.e  P  and  S  commute),  then  the  capacity  is  invariant  under  the  rotations,  and 
we  obtain, 

Hd  =  5*11^  (6.21) 

Note  that  P  and  S  commute  if  and  only  if  S  fixes  the  plane  defined  by  P,  i.e.  if  S  is  a  rotation 
around  the  normal  to  the  plane  defined  by  P. 

The  above  analysis  shows  that  in  a  LOS  environment,  a  system  with  a  triad  antenna  at  both 
ends  of  the  link  can  provide  robustness  (preserve  capacity)  to  any  orientation  of  the  transmit 
and/or  receive  antennas,  which  is  not  the  case  for  dual -polarized  antennas.  We  will  see  that 
this  robustness  extends  to  the  hybrid  LOS  and  rich  scattering  environments  described  by  (6.2). 


and 

Se  2  =  e2, 

(6.18) 

n^He  i 

n]He2\ 

e\'He  i 

e^He-z  ) 

(6.19) 

6.3  The  Full  Rate  STBC  for  a  Triad  Antenna 

In  this  section  we  introduce  a  STBC  which,  in  LOS  environments,  provides  full-rate  for  the 
triad  antenna  system  discussed  above.  We  emphasize  that,  although  our  system  consists  of  one 
physical  transmit  and  one  physical  receive  antenna,  the  underlying  channel  has  three-inputs  and 
three-outputs.  The  proposed  code  is  a  3  x  3  STBC  which  transmits  six  complex  information 
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symbols  {x,}?  j ,  x2  €  Z[i]  over  three  time  slots  from  three  effective  transmit  antennas.  The 
codewords  can  be  expressed  as 


X=/i 


(  x\  x2  x3\ 

jx 3  Xi  X2 


+  B2 


(  X4  X5  X6\ 
j;r6  x.1  x5 


^jx2jx3xi y  yjx5 jx6x4 y 


(6.22) 


where  the  diagonal  matrices  Bi  are 


Bi  -(1  +j)/3  +  e, 
b2  =  (-1  -2j)h  +  je\ 

with  0  =  diag(0i,02,03),0i  =  2cos(2,3r/7),  j  =  e2irl/3.  This  is  a  sub  code  of  the  3x3 
perfect  STBC  introduced  by  Oggicr  et  al  [7]  which  transmits  nine  complex  information  sym 
bols  over  three  time  slots  from  three  transmit  antennas.  However,  this  sub-code  still  provides 
rate  2  and  can  be  decoded  with  essentially  maximum  likelihood  performance  using  an  algo¬ 
rithm  introduced  in  [47]  with  complexity  G(N 3),  where  N  is  the  size  of  the  underlying  signal 
constellation. 


6.4  Simulation  Results 


We  simulated  the  performance  of  the  triad  and  dual-polarized  antenna  systems  allowing  re¬ 
ceiver  rotation.  The  simulation  included  pure  LOS  (K  ~  cx),  Ricean  (A'  =  10)  and  Rayleigh 
(A  —  0)  channels.  The  system  was  initialized  so  that  the  polarization  of  the  receiver  and  the 
transmitter  were  aligned,  1  e.  we  set  n  =  e3.  The  receiver  was  then  suhject  to  a  random  ro¬ 
tation  drawn  from  the  uniform  distribution.  The  counterclockwise  rotation  ahout  an  arhitrary 
unit  vector  (vi,v2,  v3 ),  hy  an  arbitrary  angle  a,  is  given  hy 


S  = 


/  v 1  +  (1  -  v2  )c  V\  t?2(  1  -  r)  -  vj  s  v  1 1^3  ( 1  -  c)  +  i>2s\ 
viv2(\  -  e)  4-  V3S  v\  +  (1  -  v2)c  v2vs(l  -  c)  -  v\s 


^UlV3(l  -  c)  -  V2S  V2V3(\  -  c)  +  V\ s  v2  +  (1  -  u3)c y 


where  c  =  coso  and  s  =  sin  a.  We  assume  the  mutual  coupling  between  the  three  elements 
is  the  same  both  at  the  transmitter  and  the  receiver.  That  is  the  mixing  operators  are  chosen  as 
follows: 


Mt  =  Mr 


J _ 

V  l+V 


/!  P  P\ 
P  1  P 

V  f>  7 


For  a  triad  antenna  system,  the  transmit  codewords  arc  the  3x3  STBC  given  in  (6.22)  and 
the  information  symbols  are  taken  from  a  16-QAM  constellation  As  shown  in  Figure  6.1,  the 
performance  of  the  system  is  stable  under  the  rotation  of  the  receiver  across  LOS,  Riccan  and 
pure  scattering  conditions.  The  Golden  Code  [22]  with  information  symbols  taken  from  the 
16-QAM  constellation  provides  a  baseline  for  evaluating  coded  dual-polarized  system  perfor¬ 
mance.  The  performance  degradation  under  3-D  rotation  shown  in  Figure  6.2  is  consistent  with 
the  analysis  given  in  Section  6.2. 
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Figure  6.  /  Performance  of  the  triad  system  with  and  without  receiver  rotation  for  uncorrelated 
channels,  p  —  0.3 
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Figure  6.2  Performance  of  the  dual-polarized  system  with  and  without  receiver  rotation  for 
uucoi  related  channels,  p  —  0.3 
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Figure  6.3:  Performance  of  the  triad  system  with  and  without  rec  eiver  rotation  for  correlated 
channels,  t  =  r  =  0.5  and  p  =  0.3 
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Figure  6  4'  Performance  of  the  dual-polanzed  system  with  and  without  receiver  rotation  for 
correlated  channels,  t  =  r  =  0.5  and  p  —  0.3 


Figures  6.3  and  6.4  show  the  performance  of  triad  and  dual-polarized  systems  respectively, 
in  a  correlated  environment,  with  t  —  r  -  0.5  and  p  —  0.3.  In  this  case,  the  performance  of 
the  triad  system  is  not  only  robust  to  the  rotation  of  the  antenna  but  is  also  robust  to  the  prop¬ 
agation  environment.  In  comparison,  the  dual-polanzed  system  loses  robustness  to  changing 
propagation  environments. 
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Figure  6.5:  Performance  of  the  triad  system  with  and  without  receiver  rotation  in  correlated 
channels ;  t  —  r  =  0.6  and  p  =  0. 1 


SNR(dB) 


Figure  6.6.  Performance  of  the  dual-polarized  system  w  ith  and  without  receiver  rotation  in 
correlated  channels,  t  =  r  =  0.6  and  p  =  0.1 


If  the  polarimetric  antennas  have  good  cross  polarization  discrimination  (XPD)  (small  value 
of  p),  even  in  highly  correlated  channels  as  shown  in  Figure  6.5  and  Figure  6.6,  where  t  —  r  — 
0  6  and  p  —  0.1,  the  triad  system  is  still  robust  across  all  conditions,  while  the  performance  of 
the  dual- polarized  system  degrades  severely  under  rotation  of  the  antenna, 
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As  shown  in  the  simulation  results,  the  use  of  triad  antenna  at  the  transmitter  and  the  receiver 
ean  provide  resilience  to  arbitrary  rotations  of  the  transmit  and  receive  antennas  as  well  as 
stahility  in  performance  across  propagation  conditions.  In  the  next  section  we  consider  the 
multi-user  downlink  where  both  the  base  station  and  each  user  is  equipped  with  a  triad  antenna, 

6.5  Fully  Polarimetric  Multi-user  Downlink 

In  this  section  we  investigate  multi-user  detection  on  a  fully  polarimetric  downlink.  In  particu¬ 
lar,  we  investigate  a  system  which  uses  a  Code  Division  Multiple  Aeeess  (CDMA)  transmitting 
scheme  in  which  each  user  is  assigned  a  two  dimensional  suhspace  and  the  transmitted  infor¬ 
mation  symbol  for  each  user  is  coded  across  this  subspace.  The  received  signal  for  each  user 
after  despreading  is  equivalent  to  the  decoding  problem  of  a  single  user  with  a  spaee-time  block 
code  system  We  demonstrate  that  the  transmitted  information  symbols  can  he  chosen  to  have 
the  structure  of  a  STBC  with  full  rate,  full  diversity  and  low  complexity  decoding 

6.5.1  Multi-user  Downlink  Model 

Consider  the  downlink  of  an  M  user  MIMO  multi-user  system  where  the  base  station  is 
equipped  with  a  triad  transmit  antenna  and  eaeh  user  has  a  triad  receive  antenna  We  consider 
CDMA  as  the  multiplexing  scheme,  i.e.  eaeh  user  is  assigned  a  three  dimensional  subspaee  of 
CT,  that  is,  a  matrix  Am  6  C3xT,  sueh  that 


AmA.\  =  6mil3,  A/. 


(6.23) 


In  order  to  make  this  work  for  A/  users  we  must  have  /’  >  3A/. 

Let  $m  be  an  information  symbol  matrix  for  user  m  with  elements  being  st1  i  =  1  ...  .9. 
where  sz  e  C, 


/S\  S 2  $3\ 


Sm  =  «4  *5  *6 

^7  Sg  Sg  j 


(6.24) 


In  a  pure  scattering  environment,  in  the  case  that  H  =  H  is  a  full  rank  ehannel  matrix, 
this  provides  3  times  the  rate  that  eould  be  achieved  using  single-lransmit  and  single  receive 
antennas  [64,  65]  However  in  a  LOS  environment,  H  =  H  as  described  in  (6.5)  for  a  fully 
polarimelrie  antenna  systems,  is  a  rank  2  channel  matrix,  and  the  rate  is  only  twiee  that  of  a 
single  antenna  system. 

The  signal  transmitted  from  the  triad  at  the  base  station  is 


M 


m  1 


We  ean  write  the  spreading  matrix  of  the  fth  user  At  as  {a\ 1 ,  ae2l  a^*)1  where  a  sequence 
a[  6  C1  x  1  and  a\a  /  =  <5tJ,  then  the  transmitted  signal  (6.25)  can  be  rewritten  as 
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At  each  transmitter  element,  the  transmitted  signal  is  an  element  of  the  subspacc  of  <C3xT 
assigned  to  the  ?n,h  user.  Therefore,  the  transmitted  signal  for  the  scheme  is  equivalent  to  that 
of  spatial  multiplexing  over  three  lime  slots,  that  is,  the  number  of  time  slots  corresponds  to 
the  dimension  of  the  subspace  allocated  to  the  user. 

We  note  that  Doostnejad  et  al.,  [66]  study  a  space-time  multiplexing  structure  for  Ml  MO 
downlink  channels  in  which  they  modulated  the  information  symhols  by  a  two  dimensional 
spreading  matrix  which  is  reminiscent  of  the  class  of  linear  dispersion  codes  discussed  in  [18]. 

Assuming  that  perfect  CS1  is  known  at  the  receiver  of  each  user,  the  received  signal  of  the 
flh  user  is 

Y£  -  H,X  +  W, 

A  (6.27) 

=  Y.  +  W- 

m  1  ,m/t 

where  W  is  the  complex  Gaussian  noise  matrix  with  entries,  WI;  6  CAf(0,  a2)  and  6 
C3x3  is  the  channel  matrix  of  the  £lh  user. 

Multiplying  (6.27)  by  A\  the  £th  received  signal  is  despread, 

yt  =  YfAl  =  H(Se  +  nf.  (6.28) 

where  the  noise  nt  —  W A\  is  the  projection  of  the  noise  onto  the  user’s  subspaee  and  is  still 
white.  Equation  (6  28)  shows  that,  after  despreading,  the  decoding  problem  of  each  user  is 
equivalent  to  that  of  a  single  user  w  ith  STBC  even  though  the  two  dimensions  are  polarization 
and  the  spreading  suhspace  components,  rather  than  space  and  lime  This  suggests  that  choos¬ 
ing  the  structure  of  the  information  symbol  matrix  S  is  equivalent  to  choosing  an  STBC.  In  the 
following  sections,  we  consider  a  suitable  choice  for  hoth  the  spreading  matrix  and  the  symbol 
matrices  which  are  compatible  with  the  triad  structure  and  admit  a  fast  decoding  algorithm 

6.5.2  Symbols  Matrix  for  Multi-user  with  Triad  System 

As  shown  in  Section  6.5. 1 ,  the  decoding  prohlcm  at  the  receiver  for  each  user  after  despreading 
is  equivalent  to  that  of  decoding  for  a  single  user  w  ho  is  using  an  STBC.  For  the  case  of  a  single 
user  system  with  triad  transmit  and  receive  antennas  we  have  previously  proposed  the  use  of 
a  polarization-lime  code  [67 ]  This  code  has  full  rale  (under  LOS  conditions),  full  diversity 
and  is  fast  deeodahle.  It  is  a  sub-code  of  the  3x3  perfect  STBC  introduced  by  Oggier  et  al. 
[7]  The  proposed  code  transmits  six  complex  information  symhols  {.T!}*  l  (usually  Ar-HEX 
constellation)  where  .t,  is  an  Eiscnstcin  integer,  x,  6  (a  +  u)b\a,b  €  Z,u;  =  e2irx ^3}  over 
three  lime  slots  from  three  effective  transmit  antennas  The  symbols  matrix  for  this  code  has 
the  following  structure 


s  =  YB'+' 

i  o 


j  «3i  +  l 

*’3t  +  2 

S3i  +  3  > 

$3i  +  l 

S3t+2 

,  (6.29) 

(U>S3i4  2 

«3i  +  l  j 

where  the  diagonal  matrices  B1  arc 

B\  =  (1  +  u)/3  +  0, 

B-2  =  (— 1  —  2oj)  I 3  +  U/’B2, 
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withe  =  diag(0! ,02,03),e,  =  2coe(2*7r/7).w  =  e2l"r\ 

For  each  user,  the  despread  procedure  reduces  the  problem  to  an  equivalent  decoding  prob¬ 
lem  of  a  single  user  with  STBC.  As  shown  in  [67]  this  code  (6.29)  provides  rate  2  across 
environmental  conditions  ranging  from  pure  scattering  to  pure  LOS  environment,  as  well  as 
resilience  to  arbitrary  rotations  of  the  transmit  and  receive  antennas.  The  code  can  be  decoded 
with  essentially  maximum  likelihood  performance  using  an  algorithm  introduced  in  [47]  with 
complexity  0(N3). 


6.5.3  Decoding 

As  wc  have  seen  the  signal  received  by  the  (th  user  after  despreading  has  the  form 

y  —  TLS  +  n. 

which  can  be  rewritten  as 


(6.30) 


V  =  (l/p  3/2^3)  =  (Sl,**2,  «3)?ii  4  (54,55,56)^2  + 
where  yt  =  (yli ,  yl2,  ^3)  and 

U\  =  (Wi,GlfCi),  Ui  =  (H2.G2X2), 

where 


Hi  = 


(  bilk  11  6,2/121  6l363|\ 

W’6,3/i31  6,1611  6,2621 


^’6,2621  0*6,3631  6,i  611 


(6.31) 


(6.32) 


(6.33) 


/ 


and  similarly  for  Qt  and  C,  for  the  second  and  third  receive  antenna  respectively.  Let  c  — 
(si ,  s 2,  53)  and  d  =  (54, 55,  se).  Wc  can  write  (631)  as 


y  —  c7i\  4  dTi 2  4  n. 


(6.34) 


The  likelihood  function  associated  with  (6.34)  is 

p(y\c  d )  rx.  exp  (  III/  ~  c7<i  -  d'H2||2j  (6.35) 

Based  on  the  conditional  optimization  described  in  [46],  the  decoding  algorithm  can  be  sum 
marized  as  follows: 

lfdet(‘W2'Wj)  >  det(WiW{) 

11  -  II2 

c  =  argmin  y  -  c7i\  -  d{c)'H2 

c  ii  li 

d  =  Q(d(c))t 

otherwise, 

d  =  argmin  ||y  -  c(d)7i\  -  d'H2 1|2  ,  (6.36) 

d 

c  =  Q(c(d)), 
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where 


d(c)  =  (y  -cW,)wS(«a«Jr\ 
c(d)  =  (y  -  dH2)H\(HiH\)  ' . 


(6.37) 

(638) 


Here  Q(c)  =  (Q(ci  ),  •  •  •  ,  Q(c3)),  where  Q  is  the  quantizer  for  the  HEX  constellation. 

6.5.4  Spreading  Matrix  of  Length  3m 

Finally,  we  need  appropriate  spreading  sequences  for  this  triad  based  multi-user  scheme.  As 
such  we  would  like  to  have  spreading  sequences  with  a  length  which  is  a  power  of  three, 
rather  than  the  usual  power  of  two  Fortunately,  spreading  sequences  of  length  3m,  which  are 
analogous  to  the  Walsh-Hadamard  sequences,  can  be  constructed  as  follows  (See  for  example 
[68]). 

Begin  with  the  Z3n  consisting  of  7?i-tuples  of  elements  from  Z3,  i.e  the  integers  modulo  3. 
For  each  b  €  Z™,  we  construct  a  function  on  Z™  with  value  at  a  6  Z™  given  by 

fb(a)=ujba  (6.39) 

for  all  a  6  Z37,  where  lj  —  exp(2ni/3)  is  a  cube  root  of  unity  and  •  denotes  the  usual  dot 
product  on  Z The  functions  fb  take  values  in  the  set  {l,u;,u;2}.  This  set  of  functions  is 
orthogonal,  i.e. 


( fb  fb ')  =  ^  fb{a)fb'(a ) 
a  ez"' 

=  ^2  u;(h  h  )  a  (6-40) 

aCZ"* 

=  3mSb,b'. 


Choosing  an  ordering  for  the  elements  of  Z”1  we  can  write  these  functions  as  sequences  and 
stack  them  in  a  matrix.  For  example  for  m  =  2,  we  obtain  a  9  x  9  orthogonal  matrix  as  follows 


a 


ih = 


1  1 

1  u/ 

1  u,2 

1  1 
1  u/ 
u,2 


v; 


In  genera]  //m  is  unitary,  i,e 


Urn  Hi  =  /7^//m=3m/3-. 


(6.41) 


(6.42) 


The  spreading  matrices  Au  for  our  scheme,  for  a  maximum  of  3m  users,  consist  of  sets  of 
mutually  exclusive  rows  from  the  matrix  II  m  \  1 . 
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6.5.5  Inhomogeneous  Multi-user  Downlink  M1MO 

The  MIMO  CDMA  scheme  described  in  the  previous  section  can  be  generalized  to  the  inho¬ 
mogeneous  MIMO  downlink  where  the  base  station  is  equipped  with  nt  transmit  antennas  but 
each  user  can  be  equipped  with  a  different  number  of  receive  antennas.  Note  that  the  base 
station  is  equipped  with  a  fixed  number  of  transmit  antennas  no  matter  how  many  users  there 
are.  Hong  et  al  [69]  proposed  a  scheme  which  allows  different  number  of  receive  antennas, 
however,  the  base  station  requires  the  number  of  transmit  antennas  to  be  the  sum  of  all  users' 
receive  antennas  in  order  to  support  such  a  capability. 

We  give  an  example  of  two  transmit  antennas,  um-polarized  (vertically  or  horizontally  po¬ 
larized),  spatially  separated  at  the  base  station  and  eaeh  user  can  be  equipped  with  either  one  or 
two  receive  antennas  The  code  structure  for  the  user  is  chosen  to  be  the  structure  of  Alamouli 
signaling  [14]  or  the  2  x  2  full  rate,  full  diversity  algebraic  STBC  such  as  the  Golden  Code 
[22],  Silver  code  [27]  etc.  or  even  just  spatial  multiplexing.  Let  us  use  the  Golden  Code  for  the 
dual  antenna  users. 

The  user  with  a  single  receive  antenna,  Alamouti  scheme  is  chosen. 


S  = 


(6.43) 


For  the  users  with  two  receive  antennas,  the  Golden  Code  structure  is  chosen, 


(q(si  +  rs2) 
ia(s 3  4  lis*) 


0(43  4  TcS4) 
o(si  4  I1S2) 


(6  44) 


where  a  =  (1  4  iy)/\/Th  a  —  (1  4  ir)/\/ 5,  r  =  (1  4  \/5)/2  and  /1  =  (1  —  \/5)/2. 

For  M  users,  the  spreading  matrix  for  each  user  is  a  two  dimensional  subspace  matrix,  i.e. 
Am  €  C2xT  where  /  =  2 M  is  the  length  of  the  sequence.  The  transmitted  signal  at  the  base 
station  is  given  by 


A/ 


?n  -  1 


+  ST  a2n\ 

\  s?a\n  4  s4a?  )  ‘ 


(6.45) 


If  the  user  has  a  single  antenna  then  their  received  signal  is  given  by 


V  =  (hn,fi2i) 


Si  CL  1  4  S'2  Q>2 


M 


+  On,/*2i)  ^2  S,„Am+n ,  (6.46) 

T7i  1  ,m^f 


-S2d  1  4  Si  a2 

where  y  is  the  1  x  T  received  signal  vector.  After  despreading,  this  becomes 


/  (  S2\ 

y  =  (Aii>/*2i)  )  4  n  ,  (6.47) 

V"62  *7 


where  y  —  y(ai,a2)t  and  nr  ~  n(ai,a2)*.  Equation  (6.47)  is  the  well  known  Alamouti 
decoding  problem. 
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Similarly,  if  the  ith  user  is  equipped  with  2  receive  antennas,  after  despreading  the  decoding 
problem  becomes 

y  =  y(ai,a2y 

An  hn\  /  ct(si  +  rs2)  n(s3  +  rs4)\  ;  (6.48) 

y/112  h22J  ya(s3  +  fis^)  q(si  +hs2)J 

where  y  is  the  2  x  T  received  vector.  This  is  the  decoding  problem  of  the  Golden  Code  for  a 
single  user  which  can  be  decoded  with  an  essentially  ML  decoder  with  quadratic  complexity 
given  in  [46]. 

6.6  Summary 

We  have  quantified  the  performance  gains  that  result  from  the  introduction  of  a  triad  antenna 
at  both  the  transmitter  and  receiver  in  a  MIMO  communication  system.  We  have  shown  that, 
in  a  LOS  environment,  the  capacity  of  the  channel  is  invariant  to  the  rotation  of  the  transmitter 
and/or  receiver.  Simulation  results  show  that  system  performance  is  stable  across  a  full  range 
of  propagation  environments  from  LOS  to  pure  Rayleigh  scattering.  An  advantage  of  the  triad 
system  over  the  baseline  provided  by  dual-polarization  was  shown  to  be  a  resilience  to  arbitrary 
rotations  of  the  transmit  and  receive  antennas.  The  practicality  of  coded  transmission  for  the 
triad  system  was  shown  through  design  of  a  full  rate  3x3  STBC  that  provides  full  rate  in  LOS 
conditions  and  admits  low-complexity  decoding 

We  have  extended  the  triad  system  to  a  multi-user  scenario  in  which  we  considered  a  CDMA 
transmitting  scheme  for  the  multi  user  downlink  in  a  system  where  both  the  base  station  and 
each  user  is  equipped  with  a  triad  antenna.  The  information  symbols  for  each  user  arc  spread 
over  their  assigned  spreading  matrix  of  orthogonal  sequences.  The  received  signal  for  each 
user  after  despreading  is  equivalent  to  the  decoding  problem  of  a  single  user  using  a  3  x  3 
STBC  This  allows  the  information  symbols  for  each  user  to  be  coded  on  an  STBC  which 
achieves  full-rate,  full  diversity  and  has  a  fast  decoding  algorithm. 
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